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MAGNETORESISTANCE AND FIELD DEPENDENCE OF THE 
HALL EFFECT IN INDIUM ANTIMONIDE! 


GASTON FISCHER AND D. Kk. C. MACDONALD 


ABSTRACT 


Magnetoresistance and Hall-effect measurements in InSb are described. This 
semiconductor has charge carriers with sufficiently long mean free paths, J, 
that it is possible, even at room temperature and with available magnetic fields, 
to obtain //r values considerably greater than unity, r being the orbital radius of 
a charge carrier moving in the applied magnetic field. The classical two-band 
theory has been found to account rather well for the results up to the highest 
magnetic fields employed. A review of the underlying assumptions of this theory 
is presented, and simple formulae are derived which allow the concentrations 
and mobilities of both types of carriers to be calculated from the magnetic 
field dependence of the resistivity, py, and of the Hall-constant, Ay. The para- 
meter A = [(Ay—Ao0)/Aol/[(o”—po)/po| provides a useful means to check the 
consistency of the theory and can give some indication of the variation of the 
mobilities with the magnetic field. 


1. THE THEORY OF THE TWO-BAND SEMICONDUCTOR 

We shall describe experiments carried out on InSb, a semiconductor whose 
charge carriers have sufficiently long mean free paths / that it is possible, 
even at room temperature, to achieve //r values greater than unity* (7 is the 
orbital radius of a charge carrier moving in the applied magnetic field 7). 
Indeed the resistance and the Hall effect are highly field-dependent (see 
Figs. 1 and 2) and immediately call to mind the theoretical predictions of 
the two-band model (Wilson 1936). This model has often been proposed to 
account for the behavior of meta/s where the two bands are usually referred 
to as s- and d-bands respectively, being first suggested for the transition 
metals, although in fact it is applied more generally. However, the success 
of this model with metals has been rather limited, for in most metals the 
magnetoresistance is found to be rather linear with 7 and the Hall coefficient 
shows little dependence on // (cf., for example, Kapitza 1929; Wilson 1936, 
1953; Chambers 1956; MacDonald 1957). The discrepancy of experiments 


1Manuscript received January 14, 1958. 
Contribution from the Division of Pure Physics, National Research Council, Ottawa, 


Canada. ; ; 
Issued as N.R.C. No. 4683. Loe 
*It should be noted that the small effective mass, m, of the carriers contributes to the 


large values of //r, remembering that ?/r = —eHr/mc = b(H/c), where r is the collision time 
and 6 the mobility. 


Can. J. Phys. Vol. 36 (1958) 











528 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


on metals with theoretical predictions was recognized already by Sommerfeld 
and Bethe (1933) in their well-known Handbuch article. On the other hand, 
as we shall show, the two-band model appears to account on the whole remark- 
ably well for the conductive behavior of InSb. In view of the very considerable 
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Fic. 1. Resistivity of p-type InSb as a Fic. 2. Hall coefficient of p-type InSb as 
function of the magnetic field. The figures a function of the magnetic field. The figures 
on each curve give the temperature in on each curve give the temperature in 
degrees Kelvin. degrees Kelvin. 


amount of work already carried out by others on InSb and related compounds, 
we do not wish to claim any great originality for our measurements, but 
wish to present a discussion of Wilson’s two-band theory and to show that 
this theory appears reasonably capable of accounting for the galvanomagnetic 
properties of a semiconductor such as InSb, even in magnetic fields for which 
//r is much greater than 1. It thus appears that the basic assumptions of the 
two-band model are quite well fulfilled in InSb, with possibly one exception: 
the assumption that the mobilities of the charge carriers are independent of 
H. The present work does not extend to fields above 8 kOe., but a few measure- 
ments at much higher magnetic fields have been presented in an earlier note 
(Fischer and MacDonald 1957) and the conclusions to be drawn do not differ 
from those of the present paper. The upper limit of 8 kOe. in the present 
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work is still well above the critical field giving //r = 1 for the conduction 
electrons. 
To analyze a system of two non-interacting and tsotropic bands, with arbi- 
trary densities of electrons and holes, let us follow Chambers (1952), who has 
shown that if a; and o2 are the conductivities which would be found for each 
band separately in the field H, and 1; and A» the corresponding separate Hall 
coefficients, then the conductivity of the complete system will be: 


(1.1) 5 = aita2)"+H'or'o2"(ArtA2)" 

' oi to2+Ho102(o11 1 +024 3) 
and its Hall coefficient: 
(1.2) 4 _ O14 Lite’ do+H'or'o2’. AA» (lite) 

~ ay (a1 +02)" )?4-H'e;'92"(Ay+As) . 

If then we assume constant mobilities }; and b: in each band (independent 
of H and energy, so that o, = —meb; and «A, = 1/mec, and correspondingly 
for band 2), we have 


_ (midi nabs) "+ (H/c) )* (n1—n2) “bi “be” 


(1.3) ee ferry ‘c)” (nyb2+-n2b1) bib» 
(1.4) ae 1 yb, *—nobe *+ UT/ Cc) “(m1 —n2) La 
; ae (yb;+n2b2)° + (H/ ‘c)" (ny —n2)"by be" 
Note that if we set m. = 0 or db, = 0 these results reduce to the single-band 
formulae ¢ = —meh; and A, = 1/mec. 


In order to simplify the comparison between theory and experiment we 
shall transform (1.3) and (1.4) so as to express: 


p Or py (the resistivity) in Q-cm., 


An in cm.*/coulomb, 

N12 in cm.—*, 

bie in cm.?/v-sec., 

V5 § in oersteds, 

e in coulombs (—1.60X 10~'° Clb.). 


These transformed formulae then read: 


(1.5) _ _ Lmybitnsbe + (107° H)” (m1b2 +n2b1) by be" 
= Ps ~ @ (nibi-+-nob2) + (10° H)?(m1—n2 ‘i * 


1 mb," '— mobs "+ (107 *H)’ ‘(m1—n» )bi°b 2 


Ax = € (1b, +n2b2) + (10 °H)*(my— nz) bide" 


It is of interest to note that if we define the parameter 


oa = (59/58) 


then on this model 


(1.8) A= (ny — nN») bbs ‘(n\b,? — nebo") 


(1.6) 
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and is thus independent of HZ. This single parameter A gives a very convenient 
means of checking the validity of the theory proposed. Let us briefly review 
the assumptions underlying this treatment. The requirement of non-interaction 
between the bands should be fairly well realized, owing to the low concentrations 
of charge carriers. For the same reason we might expect that the assumption 
of constant mobilities in each band would be acceptable as the electrons will 
all be near the bottom of the conduction band and all holes near the top of 
the valence band. With metals, however, these two conditions, or their equiva- 
lent, may be less justifiable and this may contribute to the inadequacy of the 
two-band model when applied to metals. The tsotropy condition is known to be 
fairly well realized in InSb (cf. Pearson and Tannenbaum 1953; Hrostovski 
and Tannenbaum 1954; Dresselhaus et al. 1955; Herman 1955). We suppose 
therefore that the dependence of A on 7 should tell us something about our 
assumption that the mobilities are independent of H. We shall show in an 
appendix that, with reasonable assumptions, if A should vary owing to varia- 
tions of the mobilities with 7, then we expect that 


(1.9) Acr=0) > Acu=oo): 


2. DERIVATION OF THE CHARGE-CARRIER CONCENTRATIONS AND THEIR 
MOBILITIES FROM THE FUNCTIONS p(H, T) AND An(H, T) 
Starting from the functions described by (1.5) and (1.6) we shall derive 
Ny, No, 6, and by. To do so let us write 


” 1 A+(10°°H)’B 
Sy =e - - —— 
(2.1) p(H, T) = CH (10 2D’ 
i 1 E+(10-°H)°F 
2» = eee ee 
Saad Ag(, T) eG+(10°H)*K 


and we can then set up our problem so as to solve the following system of 


equations: 

(2.3) mb+tnb. = A (= VC = VG), 
(2.4) (mb2+ n2b,)bibo= B, 

(2.5, 2.5’) (n;—n2)2b2b2 = D = K, 
(2.6) 101" — nob" = Ek, 

(2.7) (1 — M2) b;2b.? = F, 


It is seen at once that for the above system to be compatible an additional 
condition must be satisfied among the coefficients 4 to K, namely: 


(2.8) (AB—D)./D+A*F—-ED = 0. 


Iquations (2.5’) and (2.8) express a strong interconnection between the pre- 
dicted behavior of the resistivity, p(J7, 7), and the Hall coefficient, Ay(H, T). 
If the two-band model considered here is to account for the behavior of 
InSb, then these two relations should be satisfied by our experimental results. 
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Assuming that our sample is p-type, as in the present work, we know that 
nN, — Nz is always negative, and the solution of our system may be conveniently 
written as: 


i GA] 


(2.9) ms, = £- ; 

, 2 ») J ’ 

(2.10) _ fees Pieae) 
2a 6 

5+, 

(2.11) b, = aoe 

—_ 5-8 

2.12) b, = SF 

with 

(2.13) a = AB—EVD-D, 

(4.14) B= E+VD, 

and 

a 6 = V/(6'+4a). 


For an n-type specimen ”,—% is always positive, and the corresponding 
formulae are obtained by changing the sign in front of /D wherever it 


appears. 


3. EXPERIMENTAL PROCEDURE ON InSb 
Our sample was a polycrystalline prism cut from a zone-purified rod of 
InSb. It appeared to be p-type and was of the following dimensions: 


length d = 1.25 cm. width w = 0.37 cm. thickness ¢ = 0.095 cm. 


The measurements of the magnetic field were made by an induction method, 
to an accuracy of about 2%. Resistivity and Hall voltage have been measured 
to an absolute accuracy of about 5% with a bridge described by Dauphinee 
and Mooser (1955), the error margin being mainly due to uncertainties in 
the geometry. The relative accuracy is raised to 2%, and for many parameters 
as well as for the calculation of A of (1.7) only this relative error matters. 
To measure the temperature we used a copper resistance thermometer made 
from wire calibrated by Dauphinee and Preston-Thomas (1954). It was 
possible to stabilize and measure the temperature within less than one degree, 
and the range covered extends from about 80° K. to 310° K. Fig. 3 gives a 
sketch of our sample holders. 

Dividing the measured Hall voltage by H yields the Hall coefficient; this 
greatly increases the uncertainty in the low field values, and as a result the 
Hall coefficient is less well known than the resistivity. For this reason more 
experimental data at low magnetic fields would have been helpful; the measure- 
ments were performed at the following field values: H = 0, 0.125, 1.5, 2.0, 3.0, 
4.0, 5.0, 6.0, 7.0, and 8.025 kOe. We have also to mention that below 170° K. 
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we had some contact difficulties and our results are much less reliable below 
that temperature. 


TO TOP OF 
CRYOSTAT 


TO TOP OF 
CRYOSTAT 
AND SPRINGS 


ony 





Fic. 3. Sample holder and container. A, non-inductive heater; B, non-inductive resistance 
thermometer; C, copper container; D, German-silver tubing; E, sample holder (lucite); F, 
sample; G, current lead; p, c, and A stand for potential, current, and Hall electrodes respectively. 


4. DISCUSSION OF THE EXPERIMENTAL RESULTS AND THE 
COMPUTATION OF THE PARAMETERS A TO K 

At first sight one expects the two-band model to account well for the curves 
displayed in Figs. 1 and 2. One feature for which the theory cannot account, 
however, strikes one immediately; this is the minimum value observed in 
Ay which, below 180° K., is found for fields of about 5 kOe. This behavior 
has also been found by Howarth et al. (1957), who have performed very 
extensive measurements of the Hall coefficient as a function of temperature 
and magnetic field, and try to account for their results in a somewhat more 
sophisticated way than we do here. The shape of the function A,g(H, T), 
as given by (2.2), only allows extremal values of Ay versus H for fields equal 
to 0 or ~, and consequently we expect our theory to be inadequate in the 
low-temperature branch. 
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Another discrepancy appears on either a plot of 4, versus H? (cf. Fig. 4), 
or of A versus H (cf. Fig. 10). Theory predicts behavior according to the full 
line of Fig. 4, while our measurements follow rather the dashed curve. Con- 
sidering a typical set of measurements at 242.4° K. we notice a departure 


3500 
“Ay 
3 
°™Zoip N | 


~ 





2000}- 


0 5 -.. 15 20 25 
H?, (kOe)* 


Fic. 4. Hall coefficient at 242.4° K. as a function of the square of the magnetic field. At 
low field strengths the experimental points fall below the theoretical curve (in full). 


from the predicted behavior in the low-field range of 0 to 3 kOe., and it is 
interesting to compare these field values with those necessary to produce an 
1,/r; ratio equal to unity for conduction electrons in InSb: 


(4.1) 1 ‘ry = H.b,.10-8, 


where H is expressed in oersteds and }; in cm.?/v-sec. Assuming for 6; a reason- 
able value of 5X10‘ cm.?/v-sec. (cf. Fig. 8) the field necessary to render 
1,/r; equal to unity is thus: 


(4.2) H = (h/r:)/(108/b;) = 2 kOe. 


We have also seen in Section 1 that the parameter A defined by (1.7) should 
be independent of the magnetic field, unless the mobilities are themselves 
field dependent, and if they are, then we expect (cf. equation (1.9) above): 


(1.9) he > &.. 


It is evident from Fig. 10 that a quite significant change of A takes place for 
fields below about 5 kOe. and that above that value A is reasonably constant. 
One might conclude from this behavior that there is a definite variation of the 
electron mobility when /;/r; varies from zero to around unity, while no further 
change occurs above that value. We should note, however, that the magneto- 
resistance (cf. Fig. 5) seems unaffected by this variation. Furthermore, on the 
same type of plot the Hall-effect measurements of Howarth et al. (1957) do 
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seem to follow quite well the predictions of the simple theory, and we feel 
that more accurate measurements would be required to settle with certainty 
this question of mobility variation with magnetic field. 





0 5 10 15 20 25 
H?, (kOe)? 


Fic. 5. Resistivity versus the square of the magnetic field at 242.4° K. 


We now face the problem of working out the coefficients 1 to A of the 
field-independent model, from measurements which are in fact field dependent. 
One method of approach is to use those measurements in the neighborhood 
of zero magnetic field. This can be done in a simple manner, using the zero 
field slopes of the functions p(/7, T) and Ay(H, T) plotted versus H?. Let us 
denote by ¢ and y, respectively, these slopes. Using the formulae of section 
2 above we then have 


(4.3) ¢= «ee AB—D; 
(4.4) y= acta} SAntth D| - A°F—ED, 
(4.5) v/¢ = —VD, 

and the formulae (2.13) to (2.14) can be rewritten as 

(4.6) a = +E /¢, 

(4.7) B = E-w/¢. 


In this way we can express ™, m2, 61, and 2 entirely by quantities corre- 
sponding to the field 1 = 0 and thus not directly affected by the field depen- 
dence of the mobilities. We have, however, lost the means of checking the 
self-consistency of the measurements by way of the relations (2.5’) and (2.8). 

Before we became aware of the possible field dependence of the mobilities 
we computed the coefficients A to K directly and independently. Two methods 
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were used, a least mean square approach and the fitting of three arbitrary 
points. The various sets of parameters which we obtained in this manner did 
not agree well with one another; for example they were very dependent on 
the choice of the three arbitrary points when the latter method was used. 
In addition it was found that these parameters were not very consistent with 
the relations (2.5’) and (2.8). Several attempts to use them to calculate the 
n; and 6; would often lead to imaginary hole mobilities }.. It now appears 
clear why this may happen. On the other hand, using the parameters valid 
in the immediate proximity of the field H = 0 gives of course a poorer represen- 
tation of the actual behavior of p(/7, 7) and Ay(/7, T) over the whole field 
range. 

When the ratio //r increases beyond about unity it is known that the 
geometry of the sample can play an important part in magnetoresistance and 
Hall effect (Welker 1954; cf. also MacDonald 1957). We have not checked 
this point directly in the present work by using samples of different sizes, but 
it is possible that our rather short specimen exhibits geometry-dependent 
results. This might explain—at least in part—the apparent inconsistency of 
our parameters .1 to K with respect to the equations (2.5’) and (2.8). 


5. COMPUTATION OF THE CHARGE-CARRIER CONCENTRATIONS AND 
THEIR MOBILITIES 

With the parameters determined as described in the preceding section 
and the formulae (2.9) to (2.12), (4.6), (4.7), and (2.15) we obtain the results 
displayed in Figs. 6 to 9. The uncertainty of the parameters in the low and 
high temperature regions, due to the smaller variation with H, reflects itself 
in the greater uncertainty of the results. Nevertheless our results are in 
general agreement with the usual findings for InSb. Even the normally difficult, 
or often ignored, determination of b: turns out satisfactorily, and it is interest- 
ing to note that we find an increase of electron mobility with temperature, 
below 300° K., which appears to be characteristic of p-type InSb (Welker 
1954). 


6. THE BEHAVIOR OF A 


Let us consider only that range of temperatures where reasonably accurate 
values of the parameters of p(/7, 7) and A,g(H, T) can be derived (say from 
180° K. to 280° K.). As we have mentioned, it seems fairly definite, in spite 
of the considerable error margins, that all A-curves of Fig. 10 bend ultimately 
toward smaller absolute values when H/ tends toward zero. Thus we find, from 
our measurements, Ap > A... And this limiting behavior is what one would ex- 
pect if any change of the scattering mechanism takes place when the magnetic 
field increases. No marked change of A occurs above about 5 kOe. for fields 
up to some 40 kOe. The influence of holes, if their scattering cross section is 
field dependent, would appear only at very high fields because of their much 
lower mobility, and the change would presumably be spread over such a wide 
range that it might hardly be noticeable. In fact a value of b: = 1500 cm.?/ 
v-sec. gives 67 kOe. as the field for which /2./rz = 1. Moreover, it appears that 
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Fics. 6, 7, 8, and 9. These graphs show respectively the electron concentration, the hole 
concentration, the electron mobility, and the hole mobility versus temperature. Possible error 
margins are indicated by vertical lines. The dashed curve of Fig. 8, above 220° K., indicates 
results published by Welker (1954). 


significant changes in the energy gap and effective masses occur at such 
high fields (Burstein e¢ a/. 1956; Keyes et al. 1956), and it would then be 
difficult to distinguish between the many causes affecting the behavior of A. 


7. CONCLUSIONS 


1. Broadly speaking the simple two-band model of Wilson accounts well 
for the behavior of InSb. 
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Fic. 10. A = [(Ag—A0)/Ao)/[(oH — 0) /po]. The full curves correspond to the present work. 
The dashed curves show part of the data of a previous note (Fischer and MacDonald 1957). 


2. It is possible, in the frame of this model, to work out the numbers of both 
types of charge carriers and their mobilities from the observed functions 
p(H, T) and Ay(H, 7). In fact either of the pairs p(#7, T) and Ap = Ax(0, 7) 
or po = p(0, T) and Ag(H, T) will suffice. The complete knowledge of both 
functions, however, allows us to check the theory more fully by way of the 
formulae (2.5’) and (2.8). 

3. An analysis of the function A(/J, 7) as defined by (1.7) gives some in- 
dication of the dependence of the charge-carrier mobilities on the magnetic 
field. 

4. In order to carry out fully the analysis outlined under 2 and 3 above, the 
measurements must be performed with great accuracy. Our present measure- 
ments cannot claim that they are accurate enough for definite conclusions 
always to be drawn. 

5. When the computation is carried out with the knowledge of the field 
dependence of only one of the two functions p(H, 7) and A,g(H, T), more 
reliable results will generally be obtained if p(#Z, T) is fully known, for p 
results from a direct measurement, while 4g is derived from the quotient of 


two measured quantities. 


We should like to thank warmly Dr. J. W. Leech for helpful discussions and 
for checking calculations in this paper. 


8. APPENDIX 


We wish to consider, under reasonable simplifying assumptions, the limiting 
relationship of A,, to Ao when some dependence of the charge mobilities on 
magnetic field is recognized. 
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We assume that the separate resistivities of each band, p; and pz, and the 
separate Hall coefficients, 14; and A», remain finite as H — ©. We assume 
further that 

Ply Plo = P24z/P2 and A1,/A1, = A>,/Ax; 
and in particular we write 
Pin/Pio = ¥ (where we assume y > 1) 


and 
A1,./A1 = a (where we assume a < 1). 


lor the over-all resistivity we write 
P-»/ Po rm Q > 1, 


The condition that Q > 1 follows from the assumptions of our model, and is 
very generally observed. Possible exceptions may occur in InSb at very low 
temperatures (Fritzsche and Lark-Horovitz 1955; Frederikse 1957), and in 
some dilute alloys of a metal such as gold, also at very low temperatures 
(cf. e.g. MacDonald 1951, where references are given). 
By some elementary algebra following equations (1.1) and (1.2) in the 
text one readily derives the relation 
) A+ Gay, 


ea) Me Vahl 8) | 


RQ Ie 


It tollows that 
(8.2) Ao > AL 


the equality sign obtaining only if a = y = 1, corresponding to the model 
with mobilities independent of the magnetic field. Reasonable values of a 
and y which might be encountered (cf. Wilson 1953) are a = 0.8 and y & 1.3. 
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A THEORETICAL RATE-AMPLITUDE RELATION IN METEORIC 
FORWARD-SCATTERING! 


C. O. HINES 


ABSTRACT 


The theory of the forward-scattering of radio waves by ionized meteor trails is 
applied to the development of a rate-amplitude relation. This relation expresses 
the anticipated occurrence rate of scattered signals which exceed a chosen 
amplitude level, as a function of that level. It is compared with preliminary 
observational data, and found to be in good agreement both qualitatively and 
quantitatively. Closest agreement is obtained only with an appropriate choice 
of two scaling factors. These provide an abstract of the observations in a form 
which is convenient for further study and interpretation. 


INTRODUCTION 

Radio signals scattered from ionized meteor trails exhibit a wide range of 
peak amplitudes, whose distribution can be determined directly from standard 
amplitude—time recordings, after appropriate analysis. The records are scaled 
at several amplitude levels in turn, and the number of excursions of the signal 
above each chosen level, during a fixed time interval, is noted and compared 
with the others. The ‘rate—amplitude’ relation thus obtained varies appreciably 
with the observational interval and with the radio frequency, but it always 
exhibits a preponderance of weak signals. The relation can be written in the 
form 


(1) Ne « A e's 


where N, is the number of signals which exceed the counting amplitude 4 ,, 

and a isa slowly varying function of A,, lying usually in the range 1 < a < 4. 
It is the purpose of the present paper to determine the theoretical analogue 

of (1), on the basis of simplifying assumptions which can be applied readily 

to varying circumstances. These assumptions lead to the intermediate con- 

clusion that 

(2) N,«As* for A.<Ap, 

N aA for Aj Ae, 


where A, is some transitional amplitude dependent on the position and 
orientation of individual trails. An integration over all observable positions 
and orientations then transforms the abrupt transition given in (2) into a 
smooth rate-amplitude relation comparable with (1). 

The form of the derived relation is found to be in excellent agreement with 
preliminary observational data, and complete congruence can be established 
if the rates and amplitudes are scaled suitably. Variations in the observations 
can then be analyzed as variations of two scaling factors, to provide a useful 
abstract of the available information. 











1Manuscript received December 19, 1957. 
Contribution from the Radio Physics Laboratory, Defence Research Board, Ottawa, 
Canada. The work was performed under project PCC No. D48-95-11-01. 
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The specific development undertaken here is applicable only to forward- 
scatter experiments in which the transmitter and receiver are separated by 
many hundreds of kilometers. The corresponding backscatter analysis would 
follow the same pattern, however, while the general case of arbitrary trans- 
mitter—receiver separation could be treated in a similar but much more tedious 
fashion. 


BASIC ASSUMPTIONS 


It will be assumed that all signals observed in practice are scattered specu- 
larly—that is, that the incident and scattered rays make supplementary 
angles with the axis of the meteor trail in each case. This implies that the only 
trails of interest are those which lie tangent to one or other of a family of 
ellipsoids of revolution, drawn about the transmitter 7 and receiver R as 
common foci. The condition will be slightly altered, however, to simplify the 
problem of geometry: the ‘cylindrical approximation’ will be adopted, whereby 
the trails of interest are assumed to be those which lie tangent to one or other 
of a family of circular cylinders drawn about the 7-R line as common axis. 
Predictions based upon this approximation have been compared with those 
based on the accurate condition in other statistical problems, and the agree- 
ment is generally good (Hines and Pugh 1956; Hines 1958). The approximation 
reduces the requisite labor by an order of magnitude at least, and in the 
present analysis it should be quite acceptable provided the transmitter— 
receiver separation is many hundreds of kilometers. 

Two distinct classes of meteor trail, termed ‘underdense’ and ‘overdense’ 
respectively, are to be considered. The distinction lies initially in the degree of 
coupling which exists between the scattering electrons of the trail, or in the 
degree of free penetration of the trail by the incident wave, but it can be 
related to the line density of electrons along the trail. 

The underdense trails, with the lower ionizations, provide signals in which 


the peak received power is 
(3) Py = PyRKi\G@G7Gpe sina rere (rp+rpe)— (1 —cos’B sin’)! 


(Eshleman and Manning 1954). Here, Pz is the transmitted power; A, is a 
function of the radio frequency alone (given below) ; g is the number of electrons 
per unit length of the trail at the point S where the specular condition is 
fulfilled; rp and rg are the distances 7S and RS respectively; Gr and Gz are the 
corresponding antenna gains, relative to isotropic radiators; @ is the angle 
between the incident (linearly polarized) electric vector at S and the direction 
SR; B is the angle between the axis of the trail and the plane 7'SR; and 2¢ is the 
angle TSR. 

The peak signal power received from overdense trails is given by a similar 
relation, 
(4) Py = PpK2qQ'G7Gp sin?a rp rg" (rp + rg)—(1 —cos?8 sin’p)-! 


(Hines and Forsyth 1957), which differs only by virtue of the A factor and 
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the exponent of g. Although (4) is not as well founded as (3), it provides the 
only practical point of departure for the present analysis.* 

It will be assumed that (3) applies exactly for all g < q7z, and that (4) 
applies exactly for all g > qr, where q7 is a transitional value of g chosen to 
equalize (3) and (4): 

(5) Qr = (K2/K,)*/*. 
This value may be determined from the relations 
(6) Ky = \3yee4/256r4m 2, 

Ke = V9 yp'e/6447 bm? 
in which XQ is the radio wavelength, uo is the permeability of free space (in a 
rationalized system of units), e is the electronic charge, m, is the electronic 
mass, and e is the base of the natural logarithms. From these, gr = 0.75 X 10" 
electrons/meter. This estimate of the transitional g is somewhat less than 
others which are quoted frequently, and which are based on other criteria, 
but any single choice is artificial and therefore somewhat arbitrary. The 
criterion adopted here seems to be the one most suited to the present analysis, 
since it provides a continuous variation in peak received power. The dis- 
continuity implicit in the assumed rate of variation cannot be expected in 
practice, of course. However, backscatter measurements made with high 
gain antennas during a period of shower activity, and therefore concerned 
primarily with trails which were localized in position and orientation, indicate 
that the transition is complete within a range of two-to-one in signal amplitude 
(cf. Kaiser 1953). The observed forward-scatter transition for sporadic meteors 
(random orientations, broad-beamed antennas) extends over an amplitude 
range of ten-to-one or more, so the foregoing approximation for the elementary 
transition is reasonable. 

It will be assumed that the maximum line density of electrons, and the 
distribution of electrons along the trail, are governed by the relations 


(7) max = K3mcost 
and 
(8) ¢ = Codecs 


respectively, where K3 is a constant, m is the mass of the parent meteoric 
particle, and ¢ is the zenith angle of the meteor radiant. Z is a rapidly varying 
function of height and a very slowly varying function of trail orientation, 
derived by Herlofson (1947). Its dependence on orientation will be ignored 
here, and the form 


9 mes) a o( to) [ 
(9) Z= exp Hi 1 3 XP 7 





*Equation (4) is equivalent to equation (18) of the paper cited. A correction to equation 
(17) of that paper should be noted: a bracket should be inserted at the end of the first line, and 
another at the start of the second line (following the multiplication sign). 
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employed by Eshleman (1957) will be adopted. Here, ho is a constant height 
at which maximum ionization is produced, J/ is the local scale height of the 
atmosphere, which will be taken to be 7 km., and h is the height variable. 

It will be assumed that the number of meteoric particles with masses which 
exceed a given mass is inversely proportional to that mass (cf. Kaiser 1953). 
Further, it will be assumed that the radiant distribution of sporadic meteors 
is isotropic—that the meteoric particles are incident on the earth from all 
directions in equal numbers. This assumption is certainly faulty, and it 
precludes any possibility of predicting temporal variations in the rate- 
amplitude relation due to actual anisotropies, but it will serve the immediate 
purpose of this study. 

The last two assumptions can be expressed mathematically with the aid 
of two polar coordinates, 6 and y. These are used to locate radiant positions on 
that half of the celestial sphere which is visible from the mid-point of the 
forward-scatter path: 6 measures the angle between the direction of the radiant 
point and the direction of the 7—# axis, while y measures the angle between 
the vertical plane through the 7-& axis and the inclined plane in which @ is 
measured. The isotropy condition implies that the number of meteor radiants 
which lie in the element d@ dy is proportional to the solid angle sin 6 dé dy. 
When this is combined with the projection factor 


(10) cos ¢ = sin @ cosy 


(neglecting earth curvature), it implies that the number of trails that pierce a 
horizontal element of area dx dy, with orientations in the range d@ dy, is 
proportional to sin*@ cos y dé dy dx dy. (The angle 6 now measures the angle 
between the axis of the trail and the 7-2 axis, but it has the same value as 
before.) Finally, when combined with the assumed mass distribution, these 
relations imply that 


(11) dt*N = Kym, sin’6 cos y dé dy dx dy, 


where Ay is a constant and d‘N is the number of trails which pierce the 
horizontal area dx dy per unit time, with orientations in the range d6 dy, 
produced by meteoric particles whose masses exceed some value m,. 

Of greater immediate interest is the corresponding number of trails whose 
ionizations at some specified height exceed some specific value, g,. From (7), 
(8), and (10), this would be the number whose masses exceed 


(12) m, = q-K3'Z— csc 6 sec 
and from (11) this number is given by 
(13) d‘*N = K3K,4Zq,“ sin*6 cos*y dé dy dx dy. 


Absolute estimates of rates will require the evaluation of the product K3K,, 
and this may be accomplished by integrating (13) over the limits0 < 6 < a, 
—1/2 <  < 7/2, to get a total incidence rate on the element dx dy. The 
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integral should be evaluated at height 4o, where Z = 1, so that the maximum 
ionization densities will be involved: 


(14) (2N)nan, = (22/3)KsKaqc dx dy. 


It has been estimated that a maximum ionization of 2X10" electrons/meter 
corresponds to a fifth magnitude meteor (Millman and McKinley 1956), and 
that meteors of fifth magnitude or brighter are incident at an average rate 
which corresponds to 2X 10~* hr.-! km.~? (Watson 1941). The estimate 


(15) K3K,4 = 2X10" electrons/meter/hr./km.*/steradian 
may then be derived from (14). 


DISTRIBUTION OF POTENTIALLY OBSERVABLE TRAILS 

Potentially observable trails are taken to be those which satisfy the con- 
dition set out in the cylindrical approximation. For any given orientation, 
such trails are confined to a relatively small region of space, which will now be 
located. For this purpose, a Cartesian coordinate system may be introduced 
whose origin is at the mid-point of the 7—R axis, whose positive x axis lies 
along the 7—R axis in the direction 6 = 0, and whose positive z axis is directed 
vertically upwards. The positive y axis points in the direction 6 = 7/2 = —y, 
in accordance with ‘right-handed coordinate’ conventions. 

The ho height provides a convenient reference level, since the greatest part 
of the total ionization in meteor trails is confined to a narrow range of heights 
near it. This level may be specified by the surface z = 2, where 2 = 20(x, y) 
if a spherical earth is assumed and 2) = constant if a plane-earth approximation 
is adopted. (In the latter case, however, the approximation is improved by 
taking account of the depth of the mid-point of the 7—R axis below the 
actual surface of the spherical earth, when converting an assumed hp into the 
corresponding 20; 20 exceeds fy by this depth.) For purposes of discussion, a 
constant 29 will be assumed and trails will be located initially by the x and y 
coordinates of the point where they intersect the z = 20 level. 

In the cylindrical approximation, the point which corresponds to the specular 
point S lies in the inclined plane y tan y = 2 (cf. Hines 1955, with a change 
of coordinates). A trail which pierces the 2 level at «, y then has its effective 
point of specularity at a level given by 


(16) z= ycosysin ¥+2 sin’y 


(see Fig. 1). With the aid of this relation, trails may be located by their z 
coordinate at the specular point rather than by their y coordinate at the 2 
level, and the conversion 


(17) dy = csc ¥ sec y dz 


may be introduced. Equation (13) then transforms into 


(18) d4*N = K3K,4Zq-"' sin*6 cot y dé dy dx dz, 
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which gives the occurrence rate of trails whose orientations lie in the range 
dé dy, whose x coordinate at z = 2 lies in the range dx, whose z coordinate 






AXIS OF METEOR TRAIL 


Z, LEVEL 






A=(z, cot p, 2,) 


* B=(y,z,) 
.C=(y cos*w+z, cos w sin y, 
\ Y y cos w sin y+ z,sin*y) 


T-R AXIS 


Fic. 1. Projection*on the y, z plane, showing the geometry of meteor scattering when the 
cylindrical approximation is used to locate the specular point S and its projection C. 


at the point of specularity lies in the range dz, and whose ionization density 
at that point exceeds q,. Z is to be evaluated at the point of specularity, and 
it may be rewritten for this purpose as an explicit function of z: 


g 20-3 l es) 
€ = — ex ee — ey ini 
(19) Z 4 exp( li ) 1 1 esol i 


from (9), by setting Ao —h = zo—z. 

Equation (18) expresses the distribution of potentially observable trails 
in terms of the parameters 0, y, x, 2, and q,. (The last of these does not appear 
in differential form, since the specification that g exceeds q, already contains 
an implicit integration over the distribution of ionizations.) It is now necessary 
to determine the fraction of these trails that will in fact be observable, and 
so to determine the distribution of observable trails. 


DISTRIBUTION OF OBSERVABLE TRAILS 


Observable trails are taken to be those which are detected by means of 
signals which exceed some counting level. The analysis may be developed in 
terms of a numerical amplitude, A, related to received power P by 


(20) A = (P/P*)}, 


where P* is some standard power. The counting level may then be specified by 
an amplitude level 4... 

When rewriting P, and Po in terms of amplitude, as will now be done, two 
simplifications may be introduced. The angle 6 in (3) and (4) may be replaced 
by 6, since these two are equal when the cylindrical approximation is employed 
to locate the specular point. Again, it is permissible to evaluate both ¢ and 
the combination 


9 


(21) g? = G7Gpr sin’a rp rg "(rp +reR) 


at y tan Y = z = 2, and to ignore further variations of these quantities with 
y and zg, since any such variations are minor compared with the dominant 
variation of g with s (and with y, along y tan y = 32). 
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With these substitutions, the peak amplitude received from an underdense 
trail may be expressed as 


(22) Ay = (P7/P*)} Ky3gg(1—cos’ sin’¢)- 
and that received from an overdense trail as 

(23) Ao = (Pr/P*)* K2q'/4g(1 —cos’0 sin’¢)-?. 
The transitional amplitude 

(24) Ap = (P7r/P*)* Ky-°K2/%g(1 —cos’6 sin?¢)-? 
is obtained when g = qr. 

Observable trails, as already indicated, are those which provide amplitudes 
that exceed A,. If A, < Ar, a trail must have at the specular point a g which 
exceeds 
(25) dy = A(Pr/P*) 7K Ho —cos’6 sin’¢)? 


if it is to be observable, by virtue of (22). The occurrence rate of potentially 
observable trails that have such q’s is given by the distribution 


(26) d4N, = K3K,ZA, \(P7/P*)Kpg(1 —cos’6 sin’) sin *0 cot ¥ do dy dx dz 


by virtue of (18). On the other hand, if 4, >A7, a trail must have a q which 
exceeds at the specular point 


(27) go = A,*(P7/P*)?K2g-4(1 —cos’6 sin’)? 
if it is to be observable, by virtue of (23), and then the relation 
(28) din, = K3:KyZA—-*(P7/P*)?K2g4(1 —cos’6 sin’g)~ sin*é cot y dé dy dx dz 


gives the corresponding distribution of observable trails. 

From (26) and (28), it is apparent that the simple rate—amplitude relation 
described by (2) results from the assumptions which have been made, provided 
that attention is confined to trails characterized by 6, y, x, z parameters in an 
elementary range dé dy dx dz. The simple relation would continue to hold, 
even when applied to the whole array of trails, if the transitional amplitude 
were independent of these parameters. But such is not the case; Az = 
17(0, ¥, x) as given by (24). The general relation can then be obtained only 
by integrating (26) and (28) over ranges of 0, ¥, and x which vary with A,. A 
height-integration must also be performed, although the relevant range in 
this case is effectively independent of A,. 


TOTAL NUMBER OF OBSERVABLE TRAILS 


The height-integration may be completed first, as it is the least complicated. 
Z is the only factor in (26) or (28) which is still taken to be varying with 
height, and it may therefore be integrated separately: 


(29) f Z dz = 9H/A, 
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2 


where 2; = 29—// In 3 is the lower limit of trail formation (Z = 0). Although 
3 = © is taken to be the upper limit of trail formation, the greater part of the 
integral is obtained from a region of about 2/7 in depth, and it is on this 
basis that height variations of @ and g may be ignored. 

The 6-integration can be completed explicitly, so it will be undertaken next. 
kor this purpose, it is convenient to introduce the quantity 


(30) Ur x sec’ [1 —(P7/P*) Ky" K24/8g24 =], 


which is a function of the pair [y, x] through ¢ and g. Inspection will show that 
there exists some region of [y, x] variation for which u7? > 1, if only by virtue 
of vanishing antenna gain or distance attenuation in the g factor. This region, 
which will be termed Region I, would actually embrace the whole of the 
available range [—2/2 < yy < 2/2, — ~< x < @] if A, were sufficiently large. 
This corresponds to the physical fact that only overdense trails can be detected 
above some particular counting level. Further inspection will show that, 
for lower A,’s, a Region II will exist in which 0 < u7* < 1, and that there may 
even exist a Region III in which u7? < 0. The form of the 6-integral differs in 
each of these regions, and it is for this reason that the distinction is made. 
The differential form (28) applies whenever A, > A 7, or, by virtue of (24), 
whenever cos’@ is less than #7”. This is the case for all 6’s, when [y, x] lies in 
Region I, so in that region the @-integration of (28) will extend over the entire 
range 0 < @< 7m. In Region II, on the other hand, the condition obtains only 
for a limited range of 6’s, and the integration of (28) extends only over that 
range. The remaining range of 6’s is included by the integration of (26), which 
is applicable whenever A, < A,z and so whenever cos’?@ > u,*. Finally, in 
Region III, this last condition is obtained for all 6’s, and the integration of 
(26) extends over the entire range 0 < 6 < =. 

The actual integration may be performed most conveniently after a change 
of variable to 


(31) u = cos 8@. 


The range 0 < @ < aw then corresponds to a double traverse of the range 
1 > u > O for purposes of integration. The integrals take on standard forms: 


20 
(32) Vd) =2 | (u°—1)(1—u°b?) "du 
v1 


_ 140° i. 1 
~ 263 i-b/ 8 
in Region I, 
0 
(33) Vo(II) = 2 j (u°—1)(1—u"b’) "du 
Jur 


_ 140° o( ite22) _ 1-8 ur 
= l1—urb/ 1—u,'b’ 5°’ 
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(34) V,.(11) 


ur 
af (u?—1)(1—u7b") du 
v1 


Uy(1—uy b*) es 25°—1 arc sin(u7b) 


in Region IT, and 
0 

(35) V(Il) = af (u>—1)(1—u"b’) 4du 
1 


wo SEY 8 
= BP + pe are sin b 
in Region III. Each of the V’s is a function of y and x through 


(36) b =sing 


and, in Region II, through m7. 

The integration over y and x cannot be performed explicitly, and recourse 
must be made to numerical computations at this point. Formally, the final 
integrals are 


(37) N.(I) = K3Ki(9H/4)A 0‘ (P2/P*)’K2 f i V(I) g‘ cot y dy dx, 
I 

(38) N,(Il) = K3Ky(9H/4)A¢*(Pr/P*)*K?? Sf Vo(II) g‘ cot y dy dx 
It 


+K3Ki(9H/4)A CP y/P*)'K f V,(I1) g cot v dy dx, 


II 


(39) N.(IIL) = K3Ky(9H/4)A 5 (P7r/P*)'K} SJ V(III) g cot y dy dx, 

Ill 
where the region of integration is indicated by I, II, or III. The total occurrence 
rate of observable trails is then 


(40) N. = NAD+NCID+N (11D), 


and its variation with A, gives the theoretical rate-amplitude relation which 
has been sought. 


COMPUTED RATE-AMPLITUDE RELATION 

The computations required for the evaluation of (87)-(89) have been com- 
pleted for one specific case which is realized closely by the experimental 
circuits operated by this laboratory. A transmitter-receiver separation of 
1000 km. was assumed, and the 2 surface was chosen to lie 100 km. above a 
spherical earth. The sin’a factor in g? was altered to a cos*’y factor to take 
account of the selective reception of polarized fields; y is the angle between the 
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polarization vector incident on the trail and the component of polarization 
accepted by the receiving antenna from the direction of the trail. The antennas 
employed experimentally are five-element yagis mounted horizontally, so it was 
assumed in the computations that the radiated and accepted components of 
polarization were horizontal. Again, the antennas are directed along the great 
circle path between the terminals, and mounted at a height intended to produce 
maximum gain over the mid-point of the path at the meteor level. The sche- 
matic antenna pattern employed in earlier work (Hines 1955) was therefore 
assumed in the computations, as a good approximation; a maximum free- 
space gain factor of 10 was increased to 40 by an assumed perfect ground 
reflection. The absolute values employed for H, K,, Kz, and K3K4 were those 
which have been quoted above. 

Examination of (30) and (37)—(39) will reveal that NV, is a function of 4,, 
\, Pz, and P* only in the combination A A~?2P7? P*}, and this product might 
well be used as the independent variable in presenting the computed rate- 
amplitude relation. To illustrate representative values, however, it will be 
assumed that \ = 7.5 meters (frequency = 40 Mc./s.), Pr = 100 watts = the 
nominal radiated power in the experimental work, and P* = 2X10~- watts. 
The last choice permits A, to be interpreted as a microvolt measurement 
across a 50-ohm receiver input resistance, since this is the standard unit 
used in presenting the experimental results. 

The rate—-amplitude relation obtained in these circumstances is presented 
as a continuous curve in Fig. 2. The amplitude scale alone would have to be 
altered in this figure if other \’s or P’s were required, the rate scale alone 
would be altered if /7 or K3;K,4 were changed, while the whole relation would 
have to be recalculated if the path length, meteor height, polarization factor, 
or antenna pattern were changed radically. 

Logarithmic scales are employed in Fig. 2, in order to exhibit clearly the 
limiting behavior at large and small A,’s. As A, increases from transitional 
values, the simple relation N, « A,~‘ is approached. In fact, this relation 
obtains exactly when .1, exceeds some value A p given by 


(41) Ap = (Pp/ P*)1K 6K 2/3 [g sec >] max 


where [g sec ¢]max is the maximum value attained by g sec @ under variations 
of y and x. This amplitude, 4p, is the maximum amplitude attainable in a 
signal scattered from an underdense trail, when the trail has optimum orien- 
tation, position, and ionization (x = 0 = y—7/2 in the present case, 6 = 0, 
and g = qr). When A, > Ap, Region I exists alone and the integral contained 
in (37) is independent of A,; hence, VN, « A,~4. The value of N, corresponding 


to A, = Apis 
(42) Np = K3K,(9H/4)Ki7"K."[g sec b]ncx V(1) g‘ cot y dy dx, 


which, incidentally, is independent of \and P 7. The point P, at which A, = Ap 
and N, = Np, is indicated in Fig. 2. In the circumstances assumed there, 
Ap = 7.9 pv. and Np = 0.022 hr.-!. 

As A, decreases, the simple form N, « A,~! is approached asymptotically 
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but never attained exactly. The asymptotes to the rate-amplitude curve in a 
logarithmic plot are shown by broken lines in Fig. 2. They have slopes of —1 
and —4, and they intersect at a point Q which lies near the full curve in the 
region of greatest curvature. 
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Fic. 2. Theoretical rate-amplitude relation, computed for the specific parameters quoted 
in the text. 


As an abstract of an observed rate—amplitude variation, it appears that the 
position of the Q point would provide the most satisfactory single item of 
information. Not only does it lie in the vicinity of the major transition, but it 
is the only unique point that could be expected to emerge from an observational 
plot. Its coordinates are Ag = 1.3 wv., Ng = 31 hr.-', in the circumstances 
implicit in Fig. 2, although this value for Ag would be altered by the factor 
(\/7.5)*” for observations at wavelength \ (meters). 


COMPARISON OF THEORY AND EXPERIMENT 


Several experimental VHF forward-scatter circuits have been operated by 
this laboratory over the past year or two, under the direction of Dr. E. L. 
Vogan. Dr. Vogan has kindly provided some preliminary rate-amplitude data 
obtained from his program, for comparison with the theory at this stage. 

The data were obtained during the period Aug. 24-30, 1955, over the path 
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from Greenwood, N.S., to Ottawa, Ont., a distance of about 860 km. For 
purposes of analysis, each day was divided into four 6-hour intervals, and 
the data obtained from corresponding intervals on successive days were com- 
bined. The observed occurrence rate was corrected in the manner described 
by Vogan and Campbell (1957), to take account of the possible masking of 
one meteor signal by another. Three frequencies were employed: 49.98, 39.22, 
and 32.22 Mc./s. (A = 6.0, 7.6, and 9.3 m., respectively). 

The experimental data are displayed by vertical bars in Fig. 3. The length 
of each bar indicates the uncertainty in the measurement due to statistical 
fluctuations in the occurrence of meteors, and to an uncertainty in the correction 
factor. 


AUG. 24-30,1955 
A-49.98 Mc/s 
B-39.22 » 
C-32.22 « 


LOCAL -TIME LOCAL -TIME LOCAL-TIME LOCAL-TIME 
03-09 09-15 15-21 21-03 


NUMBER OF SIGNALS PER HOUR EXCEEDING Ac 





S5:5:.0 AA 2 S525 50 WANA 2 Ga5US CF ol IF 2: SoS ise 
COUNTING LEVEL Ac (MICROVOLTS) 


Fic, 3. Observed rates at various amplitudes, shown by vertical bars, with segments of 
the theoretical rate—amplitude relation shown in the positions of best fit. 


A comparison of the scales in Fig. 3 with those in Fig. 2 makes evident 
immediately the general agreement between experiment and theory, insofar 
as the location of the main transition region is concerned. Segments of the 
theoretical curve have been superimposed on the experimental plots in Fig. 3, 
in the positions of best fit; the theoretical form of the variation is obviously 
confirmed. 

Once the position of best fit of the curve has been established, the quanti- 
tative comparison of theory and experiment is most simply carried out by 
noting the location of the Q point of Fig. 2 in the coordinates of Fig. 3, and by 
comparing this location with that predicted theoretically. This is done in 
Table I. 

It is interesting to speculate on the cause of the variations which occur in the 
columns of this table. The basic point must be made, however, that the 
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TABLE I 
Period of 49.98 Mc./s. 39.22 Mc./s. 32.22 Mc./s. 
observation - -— - - -— 
(hours E.S.T.) Ag (uv.) Ne (hr.~) Ag (uv.) Ne (hr.~) Ag (uv.) Ne (hr) 
03-09 1.2 32 1.4 8U 1.3 170 
09-15 1.0 23 3.3 64 LO 150 
15-21 1.0 18 1.2 52 Bis 68 
21-03 Et 37 1.3 96 2 190 
8 


Predicted 0.93 31 3 31 1.8 31 


theoretical predictions result from an assumed isotropic incidence pattern of 
meteors, whereas the actual incidence pattern is known to be anisotropic. 
Anisotropies could be included in the theory by introducing an appropriate 
weighting factor, W(6, y), and they would lead to an altered rate-amplitude 
relation. This relation would vary with time, as the incidence pattern itself is 
varied. The P point would shift somewhat, and the Q point would shift even 
more. Accordingly, it would be premature to undertake any detailed examina- 
tion of the variations of -1g and Ng until the anisotropies have been taken 
into account. 

However, if Fig. 2 and Table I are accepted at face value, the following 
points may be made: Variations in the observed A g at a given frequency must 
result from variations in [g sec @]max, from (41). Such variations might result 
from the diurnal changes which are known to occur in the heights of formation 
of trails. The differences between the mean observed Ag’s and the corre- 
sponding predicted .1 g's could result from errors in the assumed antenna gains 
or radiated powers. Variations in the observed Ng, at a given frequency, 
almost certainly result primarily from diurnal variations in the absolute 
incidence rate. The frequency-variation of the observed N's is not so readily 
explained, however. Within the limits set by the basic assumptions, it can 
result only from a variation in the antenna patterns employed, as is evident 
from (42). This is a plausible explanation, particularly since the antenna gains 
enter (42) in the rapidly-varying form G7°G,’, but it may not be wholly satis- 
factory. As already indicated, the discrepancy cannot be assessed adequately 
until the anisotropic incidence patterns are incorporated in the theory. 


FURTHER APPLICATIONS OF THE BASIC THEORY 

Even at this stage, two further applications of the basic theory warrant 
consideration. They concern, respectively, the detection rate of overdense 
trails, and the spatial distribution of all observable trails. 

The foregoing analysis was concerned with the detection rate of all trails, 
both overdense and underdense, but the detection rate of overdense trails 
alone can be computed by similar means. Obviously, whenever 4, > 47 the 
contribution d*N,* to the observable overdense rate will equal the contribution 
d‘*N,. to the total observable rate, since no underdense trails contribute to the 
latter: d‘N.* = d1N, as given by (28). When 4A, < A7, however, the differen- 
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tial d*N, given in (26) represents both overdense and underdense contributions. 
It must be replaced by 


(43) d4N* = K3KyZK 24K," sin*6 cot w dé dy dx dz, 


which results from (18) when q, is replaced by qr. This is the differential which 
must be integrated over the range u? > uz? in Regions II and III, to obtain a 
replacement for the second term in (38) and for the whole of (39). 

The duration of an overdense signal is proportional to the value of g at the 
specular point. As g decreases towards g7, the duration decreases to a fraction 
of a second and the signal form is not markedly different from the actual form 
of an underdense signal. For observational purposes, an overdense signal may 
not be recognized as such unless the pertinent g exceeds g7 by an appreciable 
factor. To permit ready comparison with observation, then, the theory may 
be applied to determine the detection rate of signals produced by a q equal 
to ngr, say, where n > 1. In this case, the differential (28) would apply for 
A, > n'/44 7, while the differential (43), reduced by the factor m, would apply 
when A, < n!/4Ap. 

The detection rates of trails with g > nqr have been calculated for n = 1, 
4, and 16, under the conditions described in the preceding section. These rates 
have been expressed as percentages of the total detection rate, and they are 
illustrated in this form in Fig. 4. A portion of the rate-amplitude curve has 
been superimposed, for direct comparison of the amplitudes at which the 


\ 
x 100 % 
xX 
XN 
\ * 
\ 
50 
° 
0.1 10 10 


COUNTING LEVEL Ac 


Fic. 4. Percentages of observable trails having g > mgr at the point of specularity, for 
n = 1, 4, and 16, shown as functions of the counting amplitude, with a portion of the rate- 
amplitude relation indicated bv a broken curve for comparison. 


different variations occur. The overdense fraction of observable trails has not 
been determined experimentally with any precision, but rough counts (prob- 
ably comparable to n = 4 or 5) have indicated something of the order 20% 
to 40% in the region of most rapid transition. Such results are evidently 
consistent with the theoretical prediction. 

The second application to be considered here concerns the spatial distri- 





HINES: METEORIC FORWARD-SCATTERING 553 


bution of observable trails. This problem has been treated extensively in the 
past (Eshleman and Manning 1954; Pugh 1956; Hines and Pugh 1956) for 
isotropic incidence, but always with the assumption that all trails could be 
treated as if underdense. That restriction can now be removed. 

The analysis requires only a recasting of (37)—(40). In equations (37)-(39), 
y should be treated as a function of y, through the relation y tan Y = 2 which 
locates observable trails. If these equations are added as in (40), before in- 
tegration, then the form 


(44) N, = V(x, y) dx dy 


will be obtained; (x, y) gives the required distribution function. 

This function varies with A, in a complicated manner, making it impossible 
to give a single contour chart of even the relative distribution W(x, y)/ Ama 
The same applies to the variation with antenna pattern. The relative distri- 
bution has been plotted in Fig. 5, under the conditions previously assumed, 
for three cases of interest. These arise when A, is very large, intermediate, 


© 
S 
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Fic. 5. Plan views of the transmission path, showing the spatial distribution of observable 
trails for three counting amplitudes. The charts indicate relative numbers only, with maxima 
at the point centers and successive contours at 75%, 50%, and 25% of these maxima. 


and very small, respectively. More specifically, the first applies exactly for 
A, > Ap, the last applies asymptotically as A, is decreased below the main 
transition region, while the intermediate case applies during the main tran- 
sition (at A, = 1, in terms of the amplitude scale employed in Fig. 2). The 
results for A, very small are equivalent to those that would be obtained from 
the earlier analysis (Hines and Pugh 1956) in combination with the present 


antenna pattern. 
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CONCLUSION 

The theoretical rate—amplitude relation which has been derived here is in 
good agreement with observational data, both in the form of the variation and 
in the range of absolute values involved. It can be applied profitably to such 
data, to provide an abstract of the observations for further study and interpre- 
tation. 

A quantitative comparison with some preliminary data has revealed a 
variation of numbers with frequency which would not have been anticipated. 
The discrepancy may arise from inadequate approximations to the actual 
antenna patterns, or it may result from erroneous basic assumptions. It merits 
further study. 

Insofar as there is general agreement between theory and experiment, the 
relation (4) for the peak signal scattered by an overdense trail appears to be 
confirmed statistically and its value in the present type of analysis is estab- 
lished. This value is illustrated further by the development of new contour 
charts, which give the spatial distribution of observable trails in systems of 


low sensitivity. 
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RADIOACTIVE ISOTOPES OF RHENIUM AND OSMIUM 
FORMED BY THE BOMBARDMENT OF RHENIUM WITH 
PROTONS! 


J. S. Foster, J. W. HILBorn,? AND L. YAFFE 


ABSTRACT 


Targets of normal rhenium, and of rhenium enriched in each stable isotope, 
were bombarded by protons in the energy range 10 to 80 Mev. These produced 
osmium and rhenium isotopes of masses 185, 183, 182, 181, and 180. After 
chemical separations, the gamma rays and conversion electrons were examined 
for half-lives, and to determine certain conversion coefficients. Independent 
isomers of Os'*3 with 10 and 15.4 hour half-lives were found. Gamma rays with 
energies 168.1, 850, and 1097 kev. were assigned to the 10 hour isomer, and 
gamma rays with energies 114.6, 168.7, and 382.4 kev. were assigned to the 
15.4 hour isomer. In addition, the half-lives of Os!*, Os'®?, and Re'®* were ob- 
served to be 93.4 days, 21.0 hours, and 67.6 days respectively. Activities with 
half-lives of 19 hours, 23 minutes, 18 minutes, and 2 minutes were tentatively 
assigned to Re™, Os!'*!, Re, and Os! respectively 


INTRODUCTION 
This investigation has removed some of the uncertainties in the mass assign- 
ments, half-lives, and decay schemes for several neutron-deficient nuclides of 
rhenium and osmium. 
Active osmium isotopes were produced by (p, xm) reactions in targets of 
natural rhenium and of rhenium enriched in each of the stable isotopes. These 
were bombarded at suitable energies in the McGill cyclotron. 


EXPERIMENTAL METHODS 

(1) Instruments 

The analyses and measurements have been made with a gamma-ray spectro- 
meter and a Slatis-Siegbahn beta-ray spectrometer. 

(a) Gamma-Ray Spectrometer 

A 1} in. by 1 in. cylindrical Nal(Tl) crystal and Dumont 6292 photomulti- 
plier were used with the necessary electronic equipment to make an automatic- 
scanning, automatic-recording, single-channel pulse-height spectrometer. The 
gate was generally kept at 2% of full scale and the resolution of the Cs!* 
G61 kev. gamma ray was 8%. Well-known gamma rays were used for the 
energy calibration and the energy was found to be proportional to pulse height. 
The 511 kev. annihilation radiation from Na*® was used as a reference to 
correct for possible channel-width variations. A convenient variety of scanning 
rates was available. The maximum energies measured were in the neighbor- 
hood of 1 Mev. 


‘Manuscript received in original form December 13, 1957, and in revised form February 7, 
1958. 

A joint contribution from the Radiation Laboratory, Department of Physics, and Radio- 
chemistry Laboratory, Department of Chemistry, McGill University, Montreal, Quebec. 

2Present address: Atomic Energy of Canada Limited, Chalk River, Ontario. 
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(b) Beta-Ray Spectrometer (Slatis and Siegbahn 1949) 

This type of spectrometer has the distinct advantage of being effectively 
screened from external magnetic fields by its heavy cylindrical casing of 
Armco iron. Thus no compensation was necessary for the earth’s field, or for 
the field of the cyclotron magnet. 

The magnetic field within the spectrometer was generated by a symmetrical 
system of internal coils excited by a suitable range of currents. Thus electrons 
of given energy from a point source near one end of the axis may be caused 
to spiral around to a central ring focus, and to continue on to a small detector 
near the other end of the axis. The anthracene detector was 32 in. from the 
source and was screened from direct gamma rays by 1 foot of lead. This crystal 
was a polished button 0.200 in. in diameter and 0.030 in. thick. Light pulses 
from the crystal were carried through the heavy iron pole piece by a tapered 
lucite rod 83 in. long and 2 in. in diameter at the photomultiplier tube. 

Sources were mounted on a 5 pg./cm.? VYNS film (Pate and Yaffe 1955) 
which was spread over a } in. hole in a thin lucite disk. Sources were deposited 
either as a solution or as a slurry, the final diameter being in. or less. 

The instrument had an automatic scanning mechanism which continuously 
varied the magnet current. At the same time, the counting rate was measured 
with a counting-rate meter, and the resulting spectrum continuously recorded 
on a strip chart. Calibration was carried out with the well-known Cs!" and 
Th-F lines. Observations on other known lines from 500 to 10,000 gauss-cm. 
indicated a probable error in momentum of +0.2%, provided thin sources 
were used. 

While the resolution (1.8%) was satisfactory, the transmission (0.3%) in- 
dicated that the shape of the magnetic field could be improved. Subsequent 
modifications increased the transmission to 1.4% with corresponding resolu- 
tion 2.1% (Epp 1955). 


(2) Cyclotron Target Materials 

The rhenium bombarded was (a) normal rhenium metal (Re!—37.07%, 
Re'*’—62.93%) ‘‘spec-pure’’ grade obtained from Johnson, Matthey and 
Mallory, (6) rhenium metal enriched* in Re'® containing 85.38% Re’ and 
14.62% Re!*’, and (c) rhenium metal* enriched in Re'!®? containing 98.22% 
Re!®? and 1.78% Re®, 


(3) Target Mounting 

About 100 mg. rhenium metal was placed in a thin-walled tube (2S aluminum 
of wall thickness 0.0015 in.). The tube fitted conveniently into a C-shaped 
target holder which was bolted to the end of the water-cooled cyclotron 


probe. 


(4) Chemical Separation 
After the irradiation, the aluminum tube was split longitudinally and the 


*We wish to express our thanks to the Stable Isotopes Research and Production Division 
of Oak Ridge National Laboratory, Oak Ridge, Tenn., for providing the enriched rhenium 
metal. 
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irradiated material transferred to a distilling flask which contained 1 mg. of 
osmium carrier as osmic acid. The sample was dissolved in nitric acid and— 
with slight heating in a current of air—osmium tetroxide was distilled into 
ice-cold 6 V sodium hydroxide. Osmium sulphide was then precipitated using 
hydrogen sulphide. The precipitate was washed and the slurry used as the 
source material. 

This procedure was found to produce osmium free from rhenium and other 
contaminating activities. However, a careful check showed that only about 
75% of the osmium was removed in this manner. For some unknown reason, 
the remainder did not oxidize readily to the octavalent state nor exchange 
with the osmium carrier. 

To produce “‘line’’ sources for the 180° beta-ray spectrograph (see later), 
the separated osmium sulphide was dissolved in a minimum of concentrated 
hydrochloric acid plus hydrogen peroxide, and the resulting solution rendered 
just alkaline with sodium hydroxide. A 1 inch length of 0.008 in. aluminum 
wire, which had been thoroughly cleaned with carbon tetrachloride, was 
fastened to the end of the spectrograph holder and suspended for 1 hour in 
the osmium solution. Under these conditions osmium metal deposited electro- 
chemically on the wire. 

EXPERIMENTAL RESULTS AND DISCUSSION 
(a) Os! 

The table of isotopes (Hollander, Perlman, and Seaborg 1953) lists Os!* 
with a ‘‘B”’ designation, i.e., element certain, mass assignment probable. This 
nuclide is characterized by 646 kev., 872 kev., and 879 kev. gamma rays 
as the main transitions after A capture, and has a reported half-life of 93.6 
days (Johns, Nablo, and King 1957). 

Targets of normal rhenium, and enriched Re, were separately bombarded 
with 10 Mev. protons and the gamma-ray transitions of the separated osmium 
examined with the gamma-ray spectrometer. The ratios of the observed 
gamma-ray intensities corresponded with the isotopic abundance of Re'™® 
in the targets. Since the activity was shown chemically to be due to an osmium 
isotope, and the proton energy was 10 Mev., the reaction must be (), 7). 
The activity was thus identified as Os!. 

The assignment was further verified by bombardments of normal and 
enriched Re'? at 20 and 25 Mev. The ratios of the gamma-ray intensities now 
corresponded to the isotopic abundance of Re'’’. At these proton energies 
the more probable nuclear reactions were (p, 2”) and (p, 3”). The former, 
however, would have produced stable Os'™ from Re! and stable Os'* from 
Re'?, The latter would have produced relatively short-lived Os'* from Re! 
(see later) and Os'® from Re'’’. It was concluded that the observed gamma 
rays at 645 kev. and 880 kev. should be assigned to Os". 

The half-life was measured by observing the decay of the 645 kev. peak over 
a period of 300 days. In order to compensate for instrumental drifts the 
511 kev. peak from the annihilation radiation of Na® was used as a reference. 
The half-life of Os'™ was found to be 93.4+ 1.8 days by a least squares analysis, 
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based on a half-life of 2.6 years for Na”, which confirms the value of Johns, 
Nablo, and King. 

The K-conversion coefficient of the 645 kev. gamma ray was measured 
using the beta-ray spectrometer. Figure 1 (upper part) shows the K and L 
conversion peaks of the 645 kev. gamma ray at 573 kev. and 632 kev. The 
K-conversion line of the 660 kev. Cs'37 gamma ray was then scanned in the 
same fashion (Fig. 1, lower part). A comparison of the same two sources was 
then made using the gamma-ray spectrometer to determine the intensity ratio 


CONVERSION LINE SPECTRUM 
93 DAY os’ AND 68 DAY RE'® 












137 
CESIUM STANDARD 


COUNTING RATE 





ELECTRON MOMENTUM 


Fic. 1. Conversion line spectrum from osmium fraction of normal rhenium target exposed 
to 25 Mev. protons. The spectrometer scan was made 1 week after chemical separation to 
allow short-lived products to decay. Hence the source contained Os!*, and Re'®, the decay 
product of 10 hour and 15 hour Os'**. The lower spectrum is from a Cs? standard used in 
measuring the conversion coefficient of the 650 kev. gamma ray from Os", 


of the respective gamma-ray peaks. Assuming the A-conversion coefficient 
of the Cs!*7 line to be 0.087 (Wagonner 1951; Azuma 1954) the A-conversion 
coefficient of the 650 kev. Os!® gamma ray was determined to be 0.0084 + 
0.0015. This indicates an £2 transition (Rose 1955). The K/L ratio from the 
beta-ray spectrometer records was between 4 and 5, and this is consistent 
with the above classification. Johns, Nablo, and King (1957) came to similar 


conclusions. 
(b) Os'83* 
As a product of the bombardment of rhenium with 25 Mev. protons, Jones- 


*These data were presented at the Chemical Institute of Canada Symposium on Radio 
chemistry, September 7-10, 1955, and published in the preprints in preliminary form 








FOSTER ET AL.: RADIOACTIVE ISOTOPES 559 
Stover (1950) reported a 12 hour osmium activity decaying by A-capture. The 
present study showed that targets of rhenium similarly bombarded yielded 
gamma-ray activities with two distinct half-lives. One set of gamma-ray 
transitions decayed with a half-life of 10 hours, the other with a half-life of 


15.4 hours. Bombardments with enriched isotopes showed that both activities 
‘87. The nuclides were chemically 


resulted from reactions on Re'* rather than Re 
identified as osmium. 

These were assigned to Os!’ since the (p, 3”) reaction is the most probable 
at a proton energy of 25 Mev. (Bell and Skarsgard 1956). A (p, 2”) reaction 
would produce stable Os'*. This assignment was later confirmed by Epp 
(1955), who determined the shape of the vield curves and compared them 
with the curves of Bell and Skarsgard. Figure 2 is a reproduction of Epp’s 
curve for the 15 hour activity. His curve for the 10 hour activity is similar 
These isomers fit into the ‘islands of isomerism’ formed by known long- 
lived isomers as systematized by Goldhaber and Hill (1952). 
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Fic. 2. Excitation function for the 380 kev. gamma ray of 15 hour Os'* produced from 
Re, The dotted curves are excitation functions for a bismuth target, normalized to the 
peak of the rhenium curve. 

Figure 3 shows a | hour scan from the sodium iodide gamma spectrometer. 
The 73 kev. peak is due to the A X-ray of rhenium, greatly attenuated by a 
cadmium absorber. The three most prominent peaks at 120, 162, and 380 kev. 
all decayed with a 15.40.83 hour half-life. The three weaker lines at 486, 
846, and 1110 kev. all decayed with a 10+1 hour halt-life. The 15.4 hour 

value was the average of two 100 hour observations of the 380 kev. gamma 
ray. The 10 hour value was the average of two 70 hour observations of the 

1100 kev. gamma ray. Figure 4 shows a beta-ray spectrometer trace of the 

conversion lines of these gamma rays. The relatively thick source contributed 

to the poor resolution. 
Normal rhenium metal powder was then bombarded tor 3 hours at 30 Mev 


This source was exposed in the spectrograph for periods of up to 80 hours 
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It would be expected to contain Os'®, Os'®, and the Re! and Re'** decay 


products, as well as Os!**. 
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GAMMA RAY ENERGY 


Fic. 3. Gamma spectrum from osmium fraction of Re! enriched target exposed to 25 Mev. 
protons. The lines at 120, 162, and 380 kev. all decay with a 15.4+0.3 hour half-life; the weaker 
lines at 486, 846, and 1110 kev. all decay with a 10+1 hour half-life. 
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Fic. 4. Conversion line spectrum from osmium fraction of normal rhenium target exposed 
to 25 Mev. protons. The lines are identified in Table I. Note, however, that the weak lines at 
89 and 145 kev. belong to Re'8, and the line at 573 kev. belongs to Os'®, Two additional 
10 hour lines not included in Fig. 3 were observed at 1025 and 1080 kev. 


The internal conversion spectrum was studied in detail (Epp 1955) using 
a 180° beta-ray spectrograph with 0.2% resolution. The spectrograph was 
accurately calibrated using a thorium (B+C+C”) source. The spectrograph 
was adjusted to record a useful energy range from 20 kev. to 320 kev. It 
required a strong line source. 

In all, 34 lines of varying intensity were observed between 42.9 and 310.7 
kev. Of these, 14 lines have been assigned to five gamma rays of Os'!® as 
indicated in Table I. Lines appearing at the following energies are as yet 
unassigned: 54.7, 58.5, 60.7, 64.4, 66.6, 75.4, 79.3, 81.6, 87.3, 88.5, 89.8, 92.6, 
108.5, 125.8, 138.7, 142.5, 150.3, 160.2, 168.6, and 170.2 kev. The error in 
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measuring energy using the 180° spectrograph was +0.2%. Gamma-ray 
energies were calculated by adding the binding energy for rhenium, 71.7 kev. 
(Hill, Church, and Mihelich 1952), to the energy of the conversion electrons. 

In addition five more lines were observed with the Slatis-Siegbahn spectro- 
meter at energies greater than 311 kev. and assigned to three additional 
gamma rays of Os!**. Gamma-ray energies were calculated as above. The 
agreement between the two columns headed AE confirmed the assignment of 
the various conversion lines to the electron shells shown in the table. 

Precise half-life measurements were made with the gamma-ray spectrometer 
only. The half-lives of conversion lines were determined accurately enough 
to assign them to either the 10 or the 15 hour activity. In the spectrograph 
the “‘half-plate’”’ technique was a valuable aid in making such an assignment. 
Half of the plate was exposed for 15 hours while the other half was shielded. 
The second half was then exposed for a time long compared with the half- 
life while the first half was shielded. The 15 hour electron lines were equally 
dense on the two halves of the plate, while the 10 hour lines were somewhat 
more dense on the first half. Intensity ratios between the conversion lines 
were estimated from the photographic plate and are only approximate. 

The 114.6, 145.5, and 168.1 kev. gamma rays have measured properties 
consistent with 71 radiation. 


(c) Os? 

Beta-ray spectrometer measurements of normal rhenium, which had been 
bombarded at 60 Mev., resulted in a complex line spectrum. One pair of lines, 
a K-L pair at 437 and 495 kev., was very clearly resolved and gave a good half- 
life measurement. A chemical separation showed that they are osmium lines 
and hence converted in rhenium. The corresponding gamma-ray energy was 
found to be 509 kev. 

Observation of the 437 kev. peak for 135 hours gave a half-life of 21.9+0.1 
hours when analyzed by the method of least squares. This is very likely the 
24 hour activity assigned by Jones-Stover (1950) to Os!'**. Evidence for a 
12.7 hour rhenium daughter (Dybvig and Pool 1950; Wilkinson and Hicks 
1950) was indicated by conversion lines at 87.8 kev. and 1.11 Mev., both of 
which grew to a maximum in approximately 22 hours. It can be shown that 
for a 21.9 hour parent and 12.7 hour daughter the maximum should occur at 
23.7 hours. The subsequent decay could not be followed long enough for a 
good half-life measurement. However, observation of the daughter activity 
over a period of 110 hours indicated that it was controlled by the parent half- 
life of 20 to 24 hours, which is consistent with the expected value, 21.9 hours. 


(d) Re 

Wilkinson and Hicks (1950) observed a rhenium activity with a half-life 
of approximately 240 days resulting both from alpha-particle bombardment 
of tantalum and from proton bombardment of tungsten. Turner and Morgan 
(1951) believed that what the former had observed was really a mixture of 
Re!8 and 250 day Re'’® which yielded a half-life longer than their own estimate 
of 120-140 days. Lindner (1951), on the other hand, suggested that the’ 250 
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day radiation observed by Turner and Morgan was not associated with Re'*® 
and that the half-life of Re'®’ was either 150 days or greater than 5 years. 
Jones-Stover (1950) mentioned 120 days as the half-life, but did not specify 
whether this was a measurement or a restatement of some previous value. 
A complete investigation of the decay scheme was very recently published 
by Thulin e¢ al. (1956). They cite a value of 70+10 days for the half-life, 
obtained by Strominger. 

A rhenium target was bombarded at 25 Mev. and the osmium separated 
chemically after irradiation. The separated osmium should have contained 
mainly Os'** and Os'®. The 10 and 15 hour Os! activities were allowed to 
decay for several weeks so that the sample contained the long-lived Re'™® 
daughter and Os'®. Figure 5 shows the characteristic lines of the gamma-ray 
activities obtained and the upper part of Fig. 1 shows the corresponding 
electron conversion lines. 







2 GAMMA SPECTRUM 
93 DAY OSMiuM!®> 
68 DAY RHENIUM'®? 





COUNTING RATE 


GAMMA RAY ENERGY 


Fic. 5. Gamma spectrum from osmium fraction of normal rhenium target exposed to 
25 Mev. protons. The spectrometer scan was made 1 week after chemical separation to allow 
short-lived products to decay. Hence the source contained Os", and Re'*’, the decay product 
of 10 hour and 15 hour Os", 


In Fig. 5 the lines at 645 and 880 kev. arise from Os'**. There is an X-ray 
peak at about 76 kev., and peaks at 116 and 169 kev. which are those of 
Re'*’, The half-life of the 169 kev. peak was measured with the gamma-ray 
spectrometer over a period of 300 days. The half-life was found, by least 
squares, to be 67.6+2.5 days. This value agrees well with that recently found 
by Strominger as reported by Thulin et al. (1956). 


(e) Re!®! and Os"! 

When this work was begun no data on these nuclides had been reported. 
Recently Gallagher, Sweeney, and Rasmussen (1957) reported a 20-hour 
electron-capturing nuclide which they assigned to Re'*!. This nuclide emits 
a gamma ray of energy 365.5 kev. and was produced by the reaction Ta™® 
(a, 4n)Re'®!. 

In our work we have produced a 19 hour nuclide emitting a gamma ray 
of 367 kev. energy by the bombardment of Re'® with 60 and 80 Mev. protons 
and by the bombardment of tungsten with 25 and 35 Mev. protons. The most 
probable reactions indicate the production of Re'*! although Re'*®® cannot be 
excluded. A rapid osmium separation showed the decay of a 23 minute nuclide 
with conversion lines at 93 and 101 kev. which we have assigned to Os''. 
Growth of the 294 kev. and 352 kev. conversion lines of the 19 hour Re™ 
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daughter was observed with the maximum occurring at the calculated time 
(approximately 150 minutes after separation). 


(f) Re!®° and Os'8° 

Bombardment of rhenium with 80 Mev. protons vielded an 18 minute 
rhenium activity with conversion lines at 154 and 209 kev. A rapid chemical 
separation, and comparison of separated and unseparated fractions, showed 
that this must have grown in from an osmium isotope with half-life less than 
2 minutes. 

The same two conversion lines were observed following a bombardment of 
tungsten with 35 Mev. protons, but were not observed following a bombard- 
ment at 25 Mev. The evidence indicates Re®® as the most likely assignment 
for the 18 minute activity. 

Recently Haldar and Wiig (1957) suggested a 20 hour positron emitter as 
Re'S°. It appears that further experiments are required to assign masses 180 
and 181 with certainty. 
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A NEW BAND SYSTEM OF THE P, MOLECULE ANALOGOUS TO 
THE LYMAN-BIRGE-HOPFIELD BANDS OF N.! 


A. E. DouGLAS AND K. SURYANARAYANA RAO? 


ABSTRACT 


Five bands of a new band system of P2 have been photographed at high disper- 
sion and analyzed. The upper state of the system is a 'I, state and lies lower than 
any previously known excited singlet state. The lower state of the new system is 
the ground state of P2 and the analysis of the new bands has given improved con- 
stants for this state. The new system appears to be the analogue of the Lyman— 
Birge—Hopfield bands of Ne. The electron configuration of the low excited states 
of P2 and of related molecules is discussed. 


INTRODUCTION 

The absorption and emission spectrum of the diatomic phosphorus molecule 
has been studied quite extensively but up to the present only three band 
systems involving four states have been discovered. Two of the band systems 
have been investigated at low dispersion by Dressler (1955) and little is 
known of the upper states of these systems. The third system of P2 is an 
extensive one which has been investigated both in absorption and in emission. 
In emission, this system gives many strong bands covering the whole spectral 
range from about 3500 A to the vacuum ultraviolet. Rotational analysis of 
many of these bands has been carried out by Herzberg (1932), Ashley (1933), 
Herzberg, Herzberg, and Milne (1940), and Marais (1946). The band system 
arises from a !2,+—!Z,+ transition in which the lower state is the ground 
state of the molecule. Since all the strongest emission bands involve only 
relatively high vibrational levels of the ground state and since the absorption 
spectrum is difficult to obtain, the lowest vibrational level of the ground state 
that has been investigated is that with v’’ = 5. Up till the present work, the 
1y,+ upper state of this system was the lowest excited electronic state that 
had been found either in emission or in absorption and it has been assumed 
that this 'Y,* state is the lowest excited singlet state of Pe. 

The present work concerns a newly discovered band system of P2 which 
arises from a 'II—' transition. The lower state of the new band system is the 
ground state of P, and the 'II state is a previously unknown state which lies 
lower than the '2,,*+ state. The constants of both states have been determined 
and since this band system makes it possible, for the first time, to analyze a 
band involving the v’’ = 0 level of the ground state, the constants for this 
state have been somewhat improved. 


EXPERIMENTAL 
The spectrum was excited in a pyrex discharge tube about 12 mm. in 
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diameter through which helium containing a trace of phosphorus vapor was 
pumped. In the first experiments the discharge was excited in an electrodeless 
discharge by a 2000 Mc./sec. oscillator but later experiments showed that 
the same spectrum could be excited equally well by the usual type of positive 
column discharge between electrodes which were supplied by a 4000-volt 
60-cycle transformer. It is essential to allow enly a very small quantity of 
phosphorus vapor to enter the tube; otherwise the bands of the 'Z,+—'Z,+ 
system which occur in the same region become much stronger than the new 
bands. The bands were photographed with a 10 meter concave grating at a 
dispersion of 0.32 A/mm. 

Five bands belonging to the new system have been found. The band heads 
are at 2852, 2903, 2970, 3039, and 3112 A, with the 2970 A band being the 
strongest. It is difficult to obtain stronger spectra, which would probably 
show more bands, since the whole region is covered by a large number of 
bands of the '2—'Y system and it is only the stronger bands of the new system 


which show above these. 
RESULTS 

The wave numbers of the lines of the bands are given in Table I. The five 
bands are the 1-0, 0-0, 0-1, 0-2, and 0-3 bands. The 0-1 band is shown in 
Fig. 1. Since each of the bands has one state in common with at least one other 
band, the analysis of the bands presented no difficulty. The analysis shows 
that the bands are 'II—'Z bands of the P2 molecule with the lower state being 
the ground state of the molecule. Successive lines in each branch should show 
a 3:1 intensity ratio since the nuclear spin of phosphorus is one-half. In Table 
I only the strong lines are listed and in fact only a very few of the weaker set 
of lines can be seen at all. In the 0-1 band, which is the strongest of the group, 
the positions of the strongest of the weak set of lines are overlapped by lines of 
the 7-23 band of the '¥,*+—'S,*+ system. Thus the new bands are too weak to 
show that a 3:1 intensity alternation exists but the analysis does show quite 
definitely that all the observed lines must be given odd J numbers. 

The origins and B values of the bands are given in Table II. Bo’ was found 
from the mean values of the combination differences A,F’ (J) of the 0-0, 
0-1, 0-2, and 0-3 bands. The A doubling of the 'II state is very small and of 





TABLE II 
BAND ORIGINS AND B VALUES OF THE OBSERVED BANDS 
Band Yo B' B" 
1-0 35047 . 24 0.27267 0.30285 
0-0 34434 . 30 0.27436 0.30285 
0-1 33659 . 06 0.27436 0.30138 
0-2 32889 . 59 0.27436 0.29988 
0-3 32125.89 0.27436 0.29837 


the order of the experimental error. A rough estimate is that the A doubling 
constant g = 0.00003 cm~'. The B” values for the 0-0, 0-1, 0-2, and 0-3 
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bands were found from By’ and the values of (B’—B’’) determined from the 
P and R branches of the bands. For the 1-0 band, which is rather weak, the 
upper state B value was determined from By’ and the (B,’— Bo’) value de- 
termined from the differences between the Q branch lines in the 0-0 and the 
1-0 band. The constants determined here have been combined with those 
derived by Herzberg (1932), Ashley (1933), Herzberg, Herzberg, and Milne 
(1940), and Marais (1946) to give the constants of P» listed in Table III. 
DISCUSSION 

The fact that the new bands are indeed a band system of P2 is readily estab- 
lished. The values of Bo’’, a’, and AG” (3) found from the new system are 
0.30285, 0.00148, and 775.24 cm~, which may be compared with the values 
0.30256, 0.00142, and 774.82 cm~ given by Herzberg (1950) for the ground 
state of Po. When one considers that all of the values given by Herzberg have 
been obtained by extrapolating from higher vibrational levels of the ground 
state, the agreement is remarkably good. Also the nature of the light source 
and the fact that only lines with odd J values are observed confirms the con- 
clusion that the bands are due to P». 

Since the lower state of the new system is the ground state of P»2, it must be 
a !2,* state (see Herzberg 1950, p. 343). The bands of the new system are 
‘1I—'Z bands and therefore the upper state of the new system must be a ‘Il 
state. From the observed bands it is not possible to determine whether the 
upper state is ‘II, or 'II,. While it is true that a 'II,—'Z, transition is forbidden 
for electric-dipole radiation it may still occur through magnetic-dipole and 
electric-quadrupole radiation. 

The new band system arises from the lowest of the known excited states of 
P, and occurs in a spectral region which is very convenient to observe. In spite 
of this, bands of the new system have not been observed in the absorption 
experiments of Jakowlewa (1931) or Herzberg, Herzberg, and Milne (1940). 
It therefore appears certain that the bands arise from a transition with a 
low transition probability. From an emission source, it is difficult to draw 
any definite conclusions regarding the transition probability but it is worth 
noting that diluting an active gas with a large amount of an inert gas as has 
been done here is often a very effective way of causing the emission of ‘‘for- 
bidden”’ transitions. The experimental facts are therefore quite consistent 
with the assumption that the bands result from a 'II,—'2, transition. 

The low-lying electron configurations of P, and the resulting singlet states 
are 

KKLL(o,3s)? (0,35)? (my3p)* (0,3)? f 
KKLL(oa,3s)? (043s)? (ay3p)* (0,3p) (2,3p) ‘I, 
KKLL(o,3s)* (oy3s)? (au3p)? (o93p)? (4,3p) 'Zyt 1B, 1A, 


Here it can be seen that the only low-lying 'II state is a ‘II, state. Thus it 
seems most probable that the observed bands arise from a ‘II,——'S, transition. 

The bands observed here appear to be analogous to the Lyman-—Birge- 
Hopfield bands of No, which Herzberg (1946) showed most probably in- 
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volved a 'II,—'Z, transition. This assumption has recently been proved ex- 
perimentally by Wilkinson and Mulliken (1957), who observed the O and S 
branches. Probably the only reason the band system described here was not 
discovered many years ago is the fact that it is quite difficult to obtain a long 
absorption path of diatomic phosphorus. Also the weakness of the emission 
bands prevented the observation of these bands in the past and at present 
prevents the observation of O and S branches which, in the manner of Wilkin- 
son and Mulliken, would, if observed, give definite proof of the nature of the 
upper state. 

It is interesting to note that the existence of the ‘II state at an energy of 
about 34,000 cm~! might have been predicted from the known spectra of the 
isoelectronic molecules AICI and SiS which have their first excited singlet 
states at 38,254 and 32,450 cm~! respectively. These two molecules may be 
expected to have electronic structures similar to P2 but their 'II—'Y bands 
are not forbidden since the g or u symmetry is no longer present in the hetero- 
nuclear molecules. Thus the new band system arising from the 'IT, state of P»2 
at 34,500 cm! seems to fit very well into this isoelectronic group. 

Since prior to this work the lowest excited state of Pz which was known was 
the 'Z, state, it had been assumed that the 'Z, state was lower than the 'II,. 
Thus insofar as the simple molecular orbital picture is valid it was assumed 
that the ground state of P: should be written KALL(¢,3s)? (¢,3s)? (0,3p)? 
(1,3p)* and was not closely analogous to the ground state of N» for which the 
configuration is KA(o,2s)? (o,2s)? (m,2p)* (¢,2p)?. The present work shows 


that the 'X* state of the... (o,3p)? (4,3p)* (7,3p) configuration is higher than 
the 'II, state of the ...(¢,3p) (7,3p)* (7,8p) configuration. It appears more 
reasonable therefore to write the ground state of Py, as... (a,3p)* (o,3p)*. The 


ground-state configurations of Nz and P» are therefore exact analogues. If this 
is so, one may reasonably believe that the ground states of P2+ and Siz are 
*y, and !X, rather than *IT, and *Z, as has been assumed (Herzberg 1950; 


—g « 


Douglas 1955; and Narasimham 1957). The present experimental data on 
P+ and Sis are not sufficient to determine the symmetry of the ground state. 
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SPIN-ORBIT INTERACTION IN NUCLEL 
B. P. NiGAM? AND M. K. SUNDARESAN® 


ABSTRACT 


Brueckner’s theory of nuclear structure has been used to estimate the spin 
orbit interaction of a single nucleon interacting with the core arising from the 
tensor force of the two-body interaction. Using Yamaguchi's separable two-body 
interaction, the doublet splitting for Ca*! and O" has been calculated (i) in the 
first Born approximation and (ii) with an exact solution of the reaction matrix 
in which the exclusion principle has been taken into account. It is found that 
the values of the splittings calculated exactly are smaller than the experimental 
values. 


INTRODUCTION 

Calculations of the spin-orbit interaction of a single nucleon interacting 
with a core of nucleons have been carried out previously. Two main approaches 
to this problem can be distinguished. One of these is to renounce any possi- 
bility of describing the nucleon-nucleon interaction by means of a potential 
and to obtain the spin-orbit coupling part of the two-nucleon reaction matrix 
(involved in the definition of the single-particle potential) directly through 
recourse to experimental data and to use this to estimate the spin-orbit 
interaction in nuclei. A thorough investigation using this approach is due to 
Bell and Skyrme (1956), who show that the observed doublet splittings in 
O" and N® can be obtained from the experimental two-body scattering data 
up to about 100 Mev. The other approach is to retain the description of the 
nucleon-nucleon interaction in terms of a potential and to estimate the con- 
tribution to the splitting arising from the tensor force in the potential or 
from any spin-orbit coupling term that may be present in the nucleon-nucleon 
interaction itself. This latter approach has been the basis of the work of Hughes 
and Lecouteur (1950), Keilson (1951), Elliott and Lane (1954), Blin-Stoyle 
(1955), and Kisslinger (1956). When a spin-orbit interaction is included in 
the nucleon-nucleon interaction itself, the first approximation in which the 
reaction matrix of the two-particle system is replaced by the potential itself 
is finite and the splitting can be obtained relatively easily. The second- and 
higher-order corrections to the reaction matrix are assumed small because 
of the inhibition of their effects by the exclusion principle. However, in the 
case of the tensor force, the first approximation to the reaction matrix leads 
to a vanishing result for the splitting because a sum over the spins of the 
core particles is involved. Hence one has to evaluate its contribution in the 
second order, including the effects of the exclusion principle for the nucleons 
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in the core, which is a lot more difficult to calculate accurately. Kisslinger has 
evaluated the contribution to the spin-orbit interaction arising from a two- 
body tensor force with exponential radial dependence. However, he has not 
included the effects of the exclusion principle, as according to him this should 
not cause much error because the main contribution to this interaction should 
come from nucleons with momenta near the Fermi maximum. He obtained 
results which are of the order of magnitude of the observed splittings but 
somewhat smaller. Jancovici (1957) has carried out similar calculations and 
finds however that the calculated second-order contributions are considerably 
smaller than the experimental values. Thus there is considerable uncertainty 
about the actual magnitude of the contribution to doublet splitting arising 
from a two-body tensor force owing to difficulties of correctly evaluating the 
reaction matrix. 

It is the object of this paper to estimate as accurately as possible the 
spin-orbit interaction arising from the tensor force of Yamaguchi’s separable 
two-body interaction within the framework of Brueckner’s theory of nuclear 
structure (Brueckner 1955). The advantage of the separable potential is that 
the equation for the reaction matrix can be easily solved exactly including 
the effects of the exclusion principle and also without it. A comparison can 
be made of the relative contributions from the first-order approximation 
(which in this case leads to non-vanishing results for the splitting) and from 
the higher-order terms of the reaction matrix. Further, the effect of the 
exclusion principle can be easily seen. In the first part of the paper we evaluate 
the contribution to the splitting of the level of a nucleon at the top of the 
kermi distribution from the first-order approximation to the reaction matrix. 
Since there is no repulsive core in this potential, the first approximation of 
replacing the reaction matrix by the potential matrix itself should be good 
(Bethe 1956). This is checked in the latter part of the paper by using the 
exact solution to the reaction matrix including effects of the exclusion 
principle, and the results are compared with the first-order results. The 
theory has been applied to the / = 1 and / = 3 levels of Ca* and / = 2 level 
of O17. It is assumed that the central part of the two-nucleon potential leads 
(in the effective-mass approximation) to an effective mass M*/M = 0.6. This 
value is used in the calculations below. 


FIRST-ORDER BORN APPROXIMATION 


We use the Yamaguchi separable potential given by 
ene » 1 
(1) (ki V\k’) = —57 g(k)g(k’), 


where k and k’ are the momenta of the nucleons in the center-of-mass system 
and 


g(k) = CR) + Jy Stk) TOR), 


(2) S(k) = (91 Kites 


C(k) = (k°+6°)' and 7(k) = —tk’(k°+y")~. 


(6; . 6), 
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8, y, and ¢ are constants which are determined from a fit to the low-energy 
neutron-proton scattering data and the properties of the deuteron. The 
interaction strength parameter X is calculated from the relation 


ff f.5e 
(3) a a ar rt 2 1675 ’ 


where a is the scattering length. 

Since only the tensor part of the interaction can be expected to give rise 
to a spin-orbit interaction, the interaction energy responsible for the dis- 
placement of a level becomes in momentum space (Kisslinger 1956) 


AE = (Wo, VW) 
(4) 


=> J Wo" (Ki) W-*(Ke) Ti(k, K;k’, K’) po(ki’) ¥e(ky’ )d°hi dhe dh dhe’, 


where Wo and y, are the wave functions of the extra core nucleon and the 
core nucleon respectively, and 
9 dX T(k) T(R’) 


tu pe (K-) Gn Ikxk’])5(K-K’), 


(5) Ji(k,K;k’,K’) = — 
where a summation over the spins of the core particles 2 has been carried 
out. This interaction arises only if the two nucleons considered are in the 
triplet state. Later we shall include a statistical factor 3} to take care of the 
singlet states. Here K = k,+ke, k = }(k,—k2), and k’ = 3(k,’—k,’). We 
now integrate over k»’, introducing the Fourier transforms of the core func- 
tions ¥,*(kK2) and y,(k»’). We also introduce the variables R = }(r;+1r2), 
r)= r:—f2, and following Kisslinger (1956) approximatet >> .W.*(ri)¥.-(re) by 
p(R)e-*”, where a = 4.61/Ro?, Ro = 1.18X10- A! cm. Then we obtain 
9it 

6) a = 2 f y.9(e,)p(R) explitk! —k, R)] 

X T2(K1, ki’) yo(ky’) d°ky d’ky’ d°R, 


where 
(7) J2(Ki, ki’) = far d*k exp(—ar’?+ik, .r) exp(—2ik . r) 


Boe (k,k+p) (61, [kxp]) | 
x yexpl— (P -)] rat py ty E 


with p = k,’—k;. From equation (6) we see that the contribution to AE 
will be small for large values of p = k,’—k, owing to rapid oscillations of 
the wave functions. Thus we expand the quantity inside the curly bracket of 
equation (7) in powers of » about » = 0 and retain only the leading term 
(61, [kX p])k?/(k?+?)!, which involves only the first derivative of Yo*(R)yo(R) ; 
the higher powers of p involve higher derivatives of this quantity. Then the 


tThis approximation has been discussed in detail by Kisslinger (1956) in the appendix of 
his paper and depends mainly on the fact that correlé ations between nucleons in a Fermi gas 
are limited to distances which satisfy kr|ri:—rz| « 1. 
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integration over the angular parts of k and r gives, after the variables 
= k/ky, y = kir have been introduced, 


(8) I2(k,, ky’) = — {sip} Es) I (Ri) 


with 


al or ‘ 
a 2 siti: ee ee d {dio »). ad fa 2x») 
I(ki) = 8x J dx y' F(x) J dy 2y exp( p2? )z ‘ d(2xy) aay 
and 
5 2 4 
F(x) = ae ‘) a(*) | 
1 av 


The x and y integrations were carried out numerically using Simpson’s one- 
third rule. 

Now we are interested in the splitting of the level of the extra core nucleon 
and hence we consider a nucleon at the top of the Fermi distribution. Then 


Wo(k,) and Wo(ki’) can be written as 
(9) vo(k:) = 8(ki—Re) V2(R1); —— Wo(Ki") = 6(Ry’— ke) VU(RY’), 


where FY; *(K1) and Vi(k,’) are the usual spherical harmonics involving the 
angular parts of k; and k,’ specified by the unit vectors ky and ky respectively. 
Substituting these wave functions in (6) and making use of the approximate 
expression (8) for J.(k,, ki’), we can easily integrate over k; and k,’ and we 
get, after a slight rearrangement of constant factors, 


ae t° r 41 [ie 1 dp 
- 2 Si? ) 1 | tcl Je heR) pgp Rak, 


where the upper line holds for 7 = /+} and the lower line for 7 = /—}. The 
value of \/y* can be obtained from equation (3) by using the values of 8, 
y, and ¢ given by Yamaguchi and the triplet neutron—proton scattering length. 
We get 

(11) N/y® = 1.274/27°. 


The value of J(R-) for Ca* found by numerical integration was 
T(ky) = 0.614. 


Using these values in equation (10), the expression for doublet splitting 
becomes 


dp 
2141 r jt “(ee -R) R am 


’ : ii. 
(12) A(AE) = 4 (5.561 X 10 *om.*) ; es 

| ji (ReR) R°dR 

e/i) 
where the factor ? is included in equation (12) to take account of singlet 
state collisions. For O'’, the numerical factor in equation (12) turns out to 
be ?(8.781 X 10-**) cm.* To evaluate the actual magnitude of splittings, we 
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have to make some assumptions about the density distribution in the nucleus. 
For this we take 


p(R) = p for0 << R < Ro—ahy,, 
, Si a) ia 
(13) + RAI R. for Ro—aa, < R< Ro, 
=| for R > Ro, 





where A, is the Compton wavelength of the pion divided by 27. We have 
to impose the condition fp(R)d*R = A, from which we can fix the value 
of p. We calculate splittings for two values of a in equation (14), namely, 
fora = landa = 1.5. These values are given in the third and fourth columns 
of Table I. 


EXACT SOLUTION 


The reaction matrix 7 satisfies the integral equation 


= ; * 3. (k’|V\p)(p|7\k) 
i KT |k) = (k’|Vik)+P f a’p SOIR PIE) 
(14) (k’| ( +P Y, db eB) /M* 
where ky has been written as the lower limit of the integral to take account 
of the exclusion principle. Here /* is the effective mass, which we shall 
take in the following to be 0.6 . Since (k’|V|k) is separable, the solution 
of equation (14) is (in triplet states) 


dX ___g(k’)g(k) 


(15) (k'ITIK) = —a5 TCM /M) I)’ 
where 

| c ) +7" (p) 
(16) Jk -Pf dep 


and P signifies that the principal value of this integral is taken. The evaluation 
of J(k) is straightforward and gives 





= \ xy as 3S as \ 
(17) ge JR) = 28g ig py the Bi) § 
a cnt hii 
+13 7B"! — 2 8 716 Fa(x, X1, V1)» 
where 
1 \ 
fi(x, Bi) = —3(62x4+1) Bix In} at! 4 (g x°—1) tan™ Bi( 
__ hh 1 
2(1+8,°) (87x?+1) ’ 
and 


fe(x, %1,. 71) = ~ 57 ae 4 In | po +t 3(Axy'+Brui2+Cx?+Dx1), 
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mh 
ll 


—tan*y1, ba An —3+tan- i), 





EE oR ) a B 
c=H( Gy tF)) = 6a4 78 24 er) 8B 


_ _ ker _ kr Y 
a ky ’ Bi B ’ Y1 y ’ ¥3 8’ 
and 
ee 
(18) x= tyr 5 


The calculation of the doublet splitting proceeds in exactly the same 
manner as in the first order, except that now we use the 7 matrix, equation 
(15), instead of the potential matrix V. The only change this introduces is 
that in equation (7) for J2(ki, ki’), an extra factor, [1+A(M*/M)J(k)|", 
appears under the sign of integration. After carrying out the integrations over 
the angular parts of k and r, we obtain (introducing y = kyr and x = k/k») 


(19) I2(ky,k,’) = — Lax} ki) I'(k1), 


where 


: . 
I'(k) = sx’ [ dx y°F'(x) dy 2y* exp( — 2 7’) 
0 70 ky 
ape _d ees 
dy\ y / d(2xy)\ 2xy /’ 


(20) k 5 k 2 1 
F'(x) = wh s) [1+-(#) | pak ney} 
mu av 


The double integral in equation (19) was evaluated numerically, using Simp- 
son’s rule. The factor [1+A(*/M)J(k)]“ needs careful consideration. For 
x = 1, ie. k = ky, it is zero, but very close to x = 1, on either side of this 
value, it has a singularity and the principal value of the integral is calculated 
at the point below ky. The contribution due to the singularity is estimated 
to be about nine per cent of the total value of the double integral. We obtain 


(21) I2(ky, ky’) = —lsxpp Ey) [1.080] (for Ca*’), 
(22) = _(62XP} Ks) 10.837] (for 0"), 


for a nucleon at the top of the Fermi distribution. Thus we have that for 
Ca", using the wave functions given in equation (9), the expression corre- 
sponding to equation (10) 


; ot” | _ 1] | (f dp 
2. _= 
(23a) AE 5 ME (3 [1.080] 1 (4-1) ji *(kpR) Yee aR dk, 


s 3] eo d 
2: = 9.79 Mey | : | . 5 
(230) 798X107 °°cm —141)1. Ji "(ke RK) R apie 
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and for O!’, the factor in the square bracket of equation (23a) is [0.837], so 
that the corresponding numerical factor in equation (230) is 7.595 X 10-4 cm.‘. 

If we neglect the effects of the exclusion principle, that is put ky = 0 in 
equation (14), the value for the square bracket in equation (23a) for Ca* 
is [1.011]. 


CONCLUSIONS 

The above calculations were performed mainly with a view to seeing the 
effect of the different approximations made by different people in obtaining 
the spin-orbit interaction. We see that the results for the double integral 
over r and k can be compared in different cases. We discuss the case for 
Ca‘!; similar considerations hold for O'!7. When the effects of the exclusion 
principle are taken into account and the exact solution of the reaction matrix 
is used we find that the result for the splitting is larger than the first-order 
result by 57% of the first order. Of this, we find that a contribution somewhat 
less than 10% is due to the singularity alone, arising from the vanishing of 
1+A(.M*/M)J(k) in the denominator. We find further that when the exclusion 
principle is not taken into account, no singularity occurs and the result for 
the splitting is only slightly different from that in the case with the exclusion 
principle. This result may be interpreted as meaning that the exclusion 
principle makes only a little difference in the calculations because we are 
interested in nucleons at about the top of the Fermi distribution. We have 
evaluated the double integrals numerically by using Simpson’s rule checked 
by a graphical method and we feel that the calculations warrant the general 
conclusions drawn above. 

Using the density distribution given in equation (13) with a linear fall to 
zero in a region aA, of the nuclear radius (with @ = 1.0 and 1.5) we obtain 
a table of values for the splitting energies (Table 1), where the calculations 
of this section lead to the figures in the fifth and sixth columns. 


TABLE I! 
SPIN-ORBIT SPLITTING IN MEV. 








Calculated 
First order Exact 
Observed a=1.0 a=1.5 a=1.0 a=1.5 
Ca® (J = 1) 0.5 0.03 0.09 0.05 0.14 
Ca“ (/ = 3) >2 0.29 0.24 0.51 0.42 
OW (] = 2) 5.0 0.22 0.17 0.45 0.35 





We see that the results of our calculations indicate smaller values for the 
spin-orbit splittings than are indicated experimentally. To account for the 
difference, one may perhaps have to invoke the presence of spin-orbit coupling 
in the two-body interaction itself. 
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A STUDY OF NUCLEON FORCES WITH REPULSIVE CORES 


IV. CALCULATIONS OF THE SHAPE-DEPENDENT PARAMETERS 
IN NEUTRON-PROTON SCATTERING! 


R. L. Preston? anp M. A. PRESTON 


ABSTRACT 


The magnitudes of the coefficients P and Q in the phase-shift expansion 
kcotd = —1/a+ }rok?—rePki+r°Qk®+ ... have been investigated in the case 
of repulsive-core NP interactions. For square, Gaussian, exponential, and Yukawa 
shapes, which fit experimental values of the scattering length and the effective 
range, the shape-dependent parameter P has been computed as a function of r-, 
in both singlet and triplet states. Also, for square and Yukawa shapes, which 
fit the scattering length and effective range, Q has been computed in the singlet 
state. 


A. INTRODUCTION 


It is well known that the S-wave neutron—proton scattering phase shift 6 
can be expanded as 


(A.1) kcot 6 = —1/at3rok?— Proevki+Qroek®+ .... 


A complete discussion of the derivation of this formula, together with references 
to the original papers, may now be found in an article by Hulthén and Sugawara 
(1957) in the Encyclopedia of Physics. Some of the results of the present work 
have also been quoted there. For a discussion of the low-energy properties of 
a two-nucleon interaction, it is desirable to know with some accuracy the 
value of the effective range ro. Values of the effective range deduced from 
experimental scattering cross sections, the binding energy of the deuteron, and 
photodisintegration cross sections depend upon the values assumed for the 
higher-order terms in the phase-shift expansion and, in particular, on those 
assumed for P and Q. The values of P for potentials of various shapes without 
repulsive cores have been calculated by Blatt and Jackson (1949). Because of 
our interest in repulsive-core potentials, we have investigated the effect of 
the introduction of a repulsive core on the magnitude of P and Q. The potential 
shapes used are not “cut off’’ at the core radius, r,, as is the case in previous 
reports (Shapiro and Preston 1956), but are displaced from the origin. Thus, 
in considering the Yukawa potential, we are considering a potential which 
has an attractive infinity just at r,, where the infinite repulsive potential 
starts. Since the potentials deduced from meson theories would be expected 
to cut off with a finite, though possibly large, value at the core, the potentials 
considered here may be said to cover the range of reasonable shapes. 
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P, has been calculated at various values of r, for square, Gaussian, expo- 
nential, and Yukawa shapes, which give a singlet scattering length and an 
effective range a, = —23.68X10-" cm. and ro; = 2.6X10-"% cm. P, has been 
calculated for the four shapes in agreement with a triplet scattering length 
and an effective range a, = 5.28X10-" cm. and ro, = 1.56X10—-* cm. 

B. METHOD OF CALCULATION AND RESULTS 
The parameters ro, P, and Q associated with a potential 
V = —Vof [(r—r,.)/rol, ar > 7, 
(B.1) of [¢ ai 0] c 
V =; ei 
are functions of the parameters ro’, P’, and Q’ associated with the potential 
(B.2) V’ = —Vof (r/ro), r > 0. 


V is the same potential as V’ except that it is displaced from the origin a 
distance r,. Since the radial wave function u(r) vanishes at r = r,, the effect 
of displacing the potential is to displace the associated wave function « 
through the same distance, i.e. 


(B.3) u(r—r,) = u'(r). 
Since asymptotically u(r) ~ sin(kr+6), we see that 
(B.4) 5’ = 6—kr.. 


Expressions for the unprimed parameters in terms of the primed, and the 
I I 

primed in terms of the unprimed, then follow by equating coefficients of 
equal powers of k? in the expansions of the two identities 


(B.5a) kcoté = kcot(é’+kr,), 
(B.5d) k cot 6’ = k cot(6—kr,). 


Il 


(a) Evaluation of P 

For the shapes considered P’(a’ro’) is given by Blatt and Jackson (1949). 
a is the reciprocal scattering length, 1/a. Hence, by giving a and ro suitable 
values and using the first two terms in identity (B.5)) to calculate a'r,’ at 
desired values of 7., P was evaluated by equating the coefficients of &* in 
identity (B.5a) for the four shapes in both singlet and triplet states. The 
values are tabulated in Table I and the results for the four shapes are shown 
in Figs. 1 and 2. 

As an independent check, P, has also been calculated directly for a parti- 
cular square well and a particular exponential well. If u(r) is expanded about 
zero energy, 


(B.6) u = u(r) +Rui(r)+..., 
and if its asymptotic form v is also expanded, 


(B.7) , » = vo(r) +kv,(r)+..., 





PRESTON AND PRESTON: NUCLEON FORCES. IV 581 


OF CORE RADIUS 


Core radius, 


107-8 cm. Square 


FOR THOSE 
SHAPES LISTED, WHICH GIVE fro AN 


Singlet parameter P, 


TABLE 


SINGLET AND TRIPLET SHAPE-DEPENDENT PARAMETER P AS A FUNCTION 
“DISPLACED’’ POTENTIALS OF THE 


DaVA 


Gaussian 





0.0000 — .03068 — .O1764 
0. 2000 — .03438 — .02664 
0.43691 — .03721 — .03362 
0.82743 — .03923 — .03876 
1.0000 — .03950 .03940 
1.2345 — .03951 — .08951 
Triplet parameter P, 
0 — .039 - 018 
0.1 — .0419 — .028 
0.2 0443 — .0351 
0.3 — .0462 — .0404 
0.4 — 04781 — .0443 
0.6 — .Q5007 — .0492 
0.7 — .05107 — .05073 
0.8 — .05205 — 05197 
0.9 — .04958 — .O4958 
0 
.0 
.03 
-02 
Y 
P, 01 
0 
-.01 
-.02 
*.03 
*.04 
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CORE RADIUS (107'cm.) 


Fic. 1. The singlet parameter P, plotted 
against core radius r. for several potential 
shapes. 
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AS QUOTED IN SECTION A 





xponential Yukawa 


009813 05445 
— .O010385 01654 
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— .03745 — .03601 
— .03921 — .03887 
— .03951 — .03951 
0275 1355 
0037 0785 
— .O144 0326 
— .0273 0004 
— .0366 — .0208 
— .0473 — .0438 
— .05007 — .0488 
— .05184 — .05154 
— .04958 — .04958 
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.08 
-06 
Y 
-04 
.02 
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E 
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Fic. 2. The triplet parameter P; plotted 
against core radius r- for several potential 
shapes. 
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then exact expressions for P and Q (Jackson and Blatt 1950) are 


(B.8) rP = J (ses ower, 
0 
(B.9) ro QO = J (v,"—uy)dr. 
0 
Here uo and 2 vanish at the origin and satisfy 
/ i duo, yy _ 
(B.10) oe ae +Vuo = 0, 
idm) _h 
(B.11) a +Vu, = vt 
and 
(B.12) v= 1-4, 
a 
I lr 
(B.13) 1 ==7r(ro—r) +5 —. 
2 6a 


For the square well, P was calculated by means of exact integration. Numerical 
integrations of (B.11) using Milne’s formulas as in Shapiro and Preston (1956) 
were carried out to obtain P for the exponential shape. For these two cases 
the particular well parameters and the values of P obtained are listed in 
Table II and are seen to be in agreement with the corresponding values in 


Table I. 


TABLE II 
SINGLET PARAMETERS P, AND Qs FOR TWO PARTICULAR POTENTIALS 























at ex 

Shape Vo, Mev. 10-8 cm. 107% cm. P, Qs 
Square 149.35 1.64 0.8274 —.0392 00364 
Exponential 295.2 0.4400 0.4575 —.02571 .00187 


(b) Evaluation of Q 

When coefficients of k® in (B.5a) were equated, Q was obtained in terms 
of the primed parameters. However, in this case it was necessary to calculate 
Q’(a’ro’) for the desired shapes. For the square and Yukawa shapes in the 
singlet state an approximation for Q’(a’ro’) was obtained by evaluating the 
first two terms of a Taylor series expansion of Q’ 

, , , , dQ’ | 

(B.14) O'(a'ro) = O'0)+-—— (a’ro’). 


8(a’r0’) | amo 
For the two cases the approximations obtained are: 
(B.15) Square: Q,’ = 0.001713+0.0020(a’ry’), 
(B.16) Yukawa: Q,’ = 0.0113+0.113(a’r)’). 
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For the square well Q was also evaluated by exact integrations at five 
values of r,. To obtain the Q’ expansion (B.16) for the Yukawa shape, it 
was necessary to carry out numerical integrations of differential equations 
for vo—Uo, Vi— 41, OU/ Oa, OU,/da, at a = 0. The variation of a is produced 
by altering V, thus producing different mu) and u;. The values of Q, obtained 
for the square well and for the Yukawa are tabulated in Table III. 


TABLE III 


SINGLET PARAMETER Q, FOR SQUARE AND YUKAWA SHAPES FOR 
POTENTIALS WHICH GIVE ds AND fos VALUES QUOTED IN SECTION A 


Qs, square well 





ley a = Os, 

10-8 cm. (a)Fromintegration (6b) From use ofQ,’ Yukawa 
0 OO 14, — .0012 
0.04516 OOL76 0016; 

0.1 00196 — .QUO14 
0.2 . 0024; — .OO14 
0.3 00275 — .0009 
0. 4369 003815 0031, 0001 
0.5 .0032, 0006 
0.6 00354 0016 
0.63195 00348 00364 

0.8274 00364 0036s 0029 
1 0035 
1.2200 00374 

1.2345 .0037 5 .0037 





Q, was also calculated for a particular exponential, and the calculated value 
is listed with the potential in Table II. The results for Q, are compared in 


Fig. 3. 


.004 
.003 
.002 


.001 


-,001 


-,002 
0 0.2 0.4 0.6 0.8 1.0 1.2 


CORE RADIUS (107'=cm.) 


Fic. 3. The singlet parameter Q, as a function of core radius for square and Yukawa 
attractive potentials, with one point for an exponential shape. 


In order that ro’ remain positive, it is necessary that 7, be less than a 
certain limiting value, which for the singlet state is 1.23410—%cm. We 
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have not considered values of r, greater than this. It will be observed that, 
as r, increases to its limiting value, both P and Q tend to a limit independent 
of the shape of the attractive potential. 
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THE FREEZING POINTS OF HIGH PURITY METALS AS PRECISION 
TEMPERATURE STANDARDS 


lil. THERMAL ANALYSES ON EIGHT GRADES OF ZINC WITH PURITIES 
GREATER THAN 99.99+ %! 


E. H. McLAREN 


ABSTRACT 


Extensive thermal analyses have been carried out on eight grades of high 
purity zinc with the suppliers’ analyzed impurity contents in the range <0.5 to 
<70 p.p.m. The liquidus points have been intercompared with a precision of 
about 0.0002° C. and the alloy melting ranges have been examined following 
different rates of freezing. A melting range parameter may provide a more 
sensitive index of the purity than the liquidus temperature since the shape of the 
melting curve is strongly influenced by segregation of impurities due to coring on 
freezing. Two grades of Russian high purity zinc have liquidus points and 
melting ranges so close to those of New Jersey S.P. zinc that it is difficult to dis- 
criminate among them using precision resistance thermometry. This shows that 
a reported discrepancy of about 0.1° C. in the freezing points of zinc between the 
Russian and Canadian laboratories is not due to different impurity contents of 
the zincs that were used. 


INTRODUCTION 

Precision thermal analyses on eight samples of high purity zinc procured 
from four countries (U.S., U.S.S.R., U.IK., and Canada) have been carried 
out for two reasons: to investigate a difference of about 0.1° C. in the freezing 
temperatures of zinc reported by the D.I. Mendeleev Institute of Metrology, 
Leningrad (1955), and by McLaren (1955); and to demonstrate the usefulness 
of precision resistance thermometry in the detection of soluble impurities in 
high purity metals. 

The Russian and Canadian values of the freezing points of high purity zinc 
were 419.59+0.038° C. (Int.) and 419.50+0.01° C. (Int.) respectively. The 
Canadian value was based on the freezing point of the Super Pure grade zinc 
prepared by the New Jersey Zinc Company, Palmerton, Pa., U.S.A., and was 
the average of the ‘freezing point temperatures’ given by eight standard 
thermometers, seven with National Bureau of Standards (Washington) cali- 
brations and the other with a National Physical Laboratory (Teddington) 
calibration. It seems probable that this discrepancy was cue to difficulties in 
the realization of the Boiling Point of Sulphur on the International Tem- 
perature Scale. However, an interchange of high purity zincs was arranged 
with the Russian laboratory so that an intercomparison of freezing tempera- 
tures could determine conclusively whether or not the reported large difference 
was due to impurity contents of the Russian and American zincs. Samples of 
high purity zinc from British and Canadian sources were included in the 
comparison at this laboratory. 
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McLaren (1957a, 6) has shown that precision thermal analysis is capable 
of revealing alloy structures on the thermal curves of zinc, cadmium, and tin 
of purities greater than 99.99+%. The investigation of eight grades of high 
purity zinc provides a broader basis for assessing the usefulness of resistance 
thermometry in the selection of high purity metals and affords a direct com- 
parison between the thermal effects of impurities and the chemical analyses 
of the eight alloys. 


Material 

Table I lists the origin, grade, and form of the zinc metal as received and 
Table II shows the chemical analyses supplied by the sources. Where possible, 
estimates were made of the slopes of the liquidus and solidus at the zine-rich 
end of the phase diagrams for the impurities listed in Table II in order to 
calculate the changes in the liquidus and solidus temperatures from the 
freezing temperature of the pure element and the melting ranges that would 
be expected for the eight samples; these estimates are shown in Tables III 


and IV. 


TABLE | 


IDENTIFICATION OF THE ZINC SAMPLES 





Sample Origin 


R.1 D.1. Mendeleev Institute of 


Metrology, Leningrad 
- 2 ’ 


R.2 US:S.R: 
S.4 New Jersey Zinc Company, 
5.7 Palmerton, Pa., U.S.A. 


Cl ‘The Consolidated Mining and 

C.2 Smelting Company, Trail, 
B.C., Canada 

B.1 Imperial Smelting Corpora- 

B.2 tion Limited, London, U.K. 


*See Cyr (1927), Truesdale and Edn 
tSee Twyman (1951). 


Grade 
Highest Purity 
High Purity 


Super Pure* 
Chemically Pure* 


High Purity No. 678 
Tadanac 
Spectrographically Puret 
Crown Special 


1unds (1939). 


TABLE II 


Form 


Rod, 8 mm. diam. 
Rod, 8 mm. diam. 


Crystalline aggregate 
Rod, 8 mm. diam. 


Cast slab 
Cast slab 


Rod, 7 mm. diam. 
Ball, 2 in. diam. 


IMPURITY CONTENTS OF THE ZINC SAMPLES SUPPLIED BY THE SOURCES 


Sample Fe Cd 
R.1 
R.2 <i2 <0.1 
S.4 <0.5 
ad ' 0.5 
CA 2 0.1 
C.2 1O 10 
B.1 
B.2 <30 20 


Wt. (p.p.m.) 
Cu Pb Sn 


Better than R.2 


0.01 <0.2 <0.02 
<2 

0.5 l 0.5 
0.5 <2 
10 20 
<l 
<10 <10 


- Total 
As (p.p.m. ) 


<0.5 
<0.5 
<2 
0.04 6 
<4 
50 
<1 
<70 








si RP ORR ELI 





i 
E 
y 
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TABLE III 


SLOPE OF THE SOLIDUS AND LIQUIDUS AT THE ZINC-RICH END OF THE PHASE DIAGRAMS FOR 
PARTICULAR IMPURITIES 


Solid Cc CC wt... “C./wt.%. 





Wt. % solubility AT(Zn- slope slope 
Impurity eutectic limit (wt. ) eutectic) liquidus solidus References 
Fe 0.0122 0.0015 +0.05 —4.1 —33 Edmunds (1944) 


Truesdale et al. (1%36) 
Anderson (1953) 


Cd 82.5 2:0 +154 —1.9 —77 Metals Handbook (1948) 
Cu* 1.8 2:7 -4.5 +2.5 +1.7 $3 7 ss 
Pb 0.9 +1.7 —1.9 ss bs 
Sn 91 +220 —2.4 : - 





*Peritectic at 2.7% Cu, liquidus at 1.8 Cu. 


TABLE IV 


ESTIMATED CHANGE IN THE SOLIDUS AND LIQUIDUS TEMPERATURES FROM THOSE OF THE 
PURE ELEMENT, AND THE ASSOCIATED RANGE OF MELTING TEMPERATURES FOR THE SAMPLES 
LISTED IN TABLE II 


Lowering of Lowering of Estimated 


Sample liquidus (° C.) solidus* (° C.) melting range (° C.) 
R.1 <0. 00014 <0.0014 <0.0018 
R.2 <0. 00014 <0.0014 <0.0013 
S.4 <0. Q0048 <0. 0089 <0.0084 
a8 0.0019 0.0170 0.015 
Sal <0.0011 0.0073 <0.0062 
C2 0.0073 0.11 0.10 
B.1 <0. 00019 
<0.016 <0.25 <0.23 


B.2 


*In these estimates it is assumed that the solid solubility of Pb and Sn in Zn is zero (see 
Rodda (1938) and Peirce (1923)) and that the temperature at which the bulk of these im- 
purities dissolves into liquid alloy is primarily determined by the solidus temperature of the 
solid solution containing the soluble impurities in these very dilute alloys. 


EXPERIMENTAL 

The apparatus and techniques used in the thermal analysis have been des- 
cribed previously (McLaren 1957a, 6). For the present work dry nitrogen 
atmospheres were used over the melts. Very fast freezes were made by ex- 
tracting the pyrex tubes holding the crucibles of molten zinc right out of the 
zinc furnace and allowing solidification to take place either in room air or 
in another furnace nominally held at 250° C. 

The analysis was divided into two parts; first, a survey of all samples was 
carried out to obtain freezing and melting curves under a variety of thermal 
conditions using a shortened measurement technique, and secondly a careful 
intercomparison of the plateau freezing temperatures was made using the 
complete precision techniques outlined by McLaren (1957a, 6). 


Thermal Survey 

Resistance thermometry for this survey was based on NR balances at 2 ma. 
current, single daily measurements at the triple point of water, and approxi- 
mate pressure corrections over the melts. The relative errors along a particular 
curve due to this shortened technique are less than 0.0001° C. but the relative 
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difference from curve to curve could possibly be as high as 0.0005° C. For 
survey work of this type where interest lies primarily on the form of the 
thermal curves the 2NR technique allows many measurements to be made 
since there are no undue dead times as the resistance thermometer comes to 
equilibrium at different currents. For the very fast freezes the measuring 
technique was further shortened to 2N balances, enabling several measure- 
ments to be made in a minute. Times for resistance balances were recorded 
to the nearest second. 

Considerable attention has been given to the rounding at the beginning and 
end of arrests; Hume-Rothery et al. (1952) indicate that arrests may be 
sharpened by decreasing the thermal lag of the temperature detector and by 
using slow heating rates. The Meyers resistance thermometer used in this work 
has such a rapid response time that thermal lags due to this detector im- 
mersed in the 1-kg. samples are negligible, except possibly during very fast 
freezes. The rounding at the end of the melting curves may be due to con- 
vection currents in the crucible when the bulk of the sample is liquid and only 
a small fraction of the charge is still solid; this rounding in the few minutes 
prior to the liquidus break is similar from sample to sample. The best criterion 
that good thermal conditions exist for the transformations is the close agree- 
ment between the temperature of the liquidus break on the melting curve and 
the plateau temperature (i.e. liquidus point, see McLaren (1957)) for the 
same sample on either normal or induced slow freezes. 

Freezing and remelting curves were obtained for several series of trans- 
formations which were designed to inhibit or enhance segregation of chemical 
impurities (coring) by control of cooling rates on the freeze and to detect these 
effects by examination of the remelting contour.* For each series of trans- 
formations the furnaces were adjusted to give similar thermal treatment to 
all samples; melting curves were taken immediately after completion of 
freezes to minimize annealing effects in the ingots. Table V gives a summary of 
transformation conditions for the thermal curves included in this paper. Series 
2(b) and 3 were run twice; essentially the same information was obtained both 











times. 
TABLE V 
SUMMARY OF TRANSFORMATION CONDITIONS FOR THE THERMAL CURVES 
Duration 
Series Designation (minutes) Sample Figure 
1. Freeze Fast normal 35-40 All 2 
2(a). Freeze Very fast normal 3 R.1 1 
2(b). Freeze Very fast normal 5-6 All 3 
3. Melt Slow after freezes of Fig. 3 70-100. All 4 
4. Freeze Fast induced 100) 
Slow induced C.1—150 > All a 
Slow induced C.2 -180§ 
5. Melt Slow after freezes of Fig. 5 70-100 All 6 











*Melting contour is used here to mean the profile of the direct time-temperature melting 
curve. 
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Results 

Fig. 1 shows representative freezing curves in sample R.1. On fast induced'and 
fast normal freezes, curves I and II respectively, the zinc solidifies uniformly 
in closed cylindrical shells, providing an excellent environment and shield for 
the resistance thermometer. Curves III and IV were obtained at very fast 
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Fic. 1. Representative freezing curves on sample R.1: I, fast induced; II, fast normal; 
11] and IV, very fast normal. 


freezing rates where the rapidity and lack of uniformity in freezing would; be 
expected to provide extremely unfavorable conditions for resistance ther- 
mometry; hence it is surprising that the arrests noted on these curves are 
within 0.002° C. of the plateau freezing temperatures that are obtained on 
slow induced freezes using the complete precision technique. Undoubtedly the 
very large amount of latent heat associated with 1 kg. of zinc, combined with 
the high purity of this sample, enables these fairly reliable temperature 
measurements to be made even when the ingot solidifies in 3 minutes. 

Fast normal freezing curves on all the samples, obtained by monitoring 
temperatures in the samples as they solidified in the furnaces after the power 
supplies were turned off, are shown in Fig. 2. From this figure it is evident that 
the samples can be usefully divided into three groups on the basis of maximum 
temperatures attained during recalescence; group I includes R.1, R.2, and 
S.4; group II, C.1, B.1, and 5.7; and group III, C.2 and B.2. This figure shows 
freezing occurring over a range of temperatures and distinctive differences are 
noticed among the samples; but as liquidus arrests are screened out by super- 
cooling (about 0.04° C. for these samples) it is difficult to obtain useful in- 
formation on the relative ranges of alloy transformation temperatures from 
these curves. However, the extreme sensitivity and simplicity of using the 
resistance thermometer to sort out high purity metals is demonstrated by 


these curves. 
The ideal range of alloy transformation temperatures is the difference in 
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the temperatures at which melting or freezing begins and ends in a completely 
homogenized alloy; or alternatively the difference in the temperatures of the 
solidus and liquidus for the alloy on the equilibrium diagram. In practice very 
great obstacles exist in realizing ideal melting ranges because it is difficult, if 
not impossible, to attain complete homogenization of a sample; this is a con- 
sequence of the differences in the solubilities of the impurities in the solid and 
liquid alloys and of the low rate of diffusion of the solute in the solid. Earlier 
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Fic. 2. Fast normal freezes on all samples; plotting begins during recalescence. 
Fic. 3. Very fast freezes on all samples. 
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work on New Jersey zincs has shown that the best alloy ranges (i.e. approxi- 
mating textbook thermal curves on simple binary alloys) are obtained from 
melting curves determined immediately after very fast freezes. The rapid 
freezing gives rise to smaller grains and the impurities must be more evenly 
distributed throughout the ingot than after slow freezing, which gives very 
large highly cored grains. The impurity concentrations near the grain boun- 
daries in these cored grains are many times greater than the equilibrium 
concentrations for the ingot. These non-equilibrium concentrations of dis- 
solved impurities cause differences between the true solidus temperatures of 
the alloys and the temperatures at which melting begins, but at the same time 
provide a mechanism for detecting the impurities and their segregation in 
high purity metals. 

In order to distribute the impurities throughout the ingot as evenly as 
possible, in preparation for melting range determinations, the very fast freezes 
of Fig. 3 were obtained. Surprisingly these curves deserve attention in their 
own right because they show that even under these severe conditions the 
resistance thermometer is capable of distinguishing between samples of 
different quality. Fig. 4 shows the melting curves that were obtained im- 
mediately after these freezes; the alloy ranges varied from about 0.001° C. 
for the purest to a few hundredths of a degree for the least pure samples. 

In order to accentuate segregation of the dissolved impurities the induced 
freezes shown in Fig. 5 were made and melting curves following these freezes 
are given in Fig. 6. 

A comparison of Figs. 4 and 6 shows that although group I metals have 
similar alloy melting ranges following the two types of freezing, distinct 
differences are noted for the samples of group II and III. In Fig. 6 segregation 
structures have been developed on the melting curves of samples C.1 and B.1 
and now very substantial fractions of all the group II metals melt well inside 
a range of 0.01° C. For group III samples two observations are made: the C.2 
curve does not have the step near the liquidus break that appears in Fig. 4 and 
the general melting contours for both C.2 and B.2 lie on steeper slopes than 
those obtained in Fig. 4. Melting takes place over greater ranges and only 
the temperatures of the liquidus breaks remain unchanged. 

A more detailed analysis of these melting curves allows comparative 
measurements to be made of the melting ranges of these samples. Measure- 
ments of three melting parameters were made: conventional melting ranges 
from direct melting curves, and the widths at half-maximum and at the 1% 
level from inverse curves. 

Curve I of Fig. 7 reproduces on a larger scale the melting curve of S.4 shown 
in Fig. 4 and shows how the melting ranges of the alloys were estimated from 
the direct melting curves. 

Curve II of Fig. 7 shows the inverse melting curve of I. This was derived 
from the times for successive equal temperature intervals and from the total 
melting time; the percentage of the total melting time in a given temperature 
interval is plotted against the average temperature of the interval. For the 
slow heating rates used in these experiments the furnace block temperature 
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Fic. 4. Melting curves after the very fast freezes of Fig. 3. 
Fic. 5. Induced freezes on all samples. 
Fic. 6. Melting curves after the induced freezes of Fig. 5. 
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remains nearly constant during the transformation; therefore the percentage 
time in a given temperature interval is a fair approximation of the percentage 
of the sample that is melted during the interval. The inverse curve does not 
give as detailed information as the direct curve but it has the advantage’of 
enabling the essential melting range information to be plotted in a smaller 
space. On curve II of Fig. 7 an estimate of the width at half-maximum ‘W’ 
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Fic. 8. Inverse melting curves of samples R.1, R.2, and S.4: I, after very fast freezes; 
and II, after induced freezes. 
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is shown; both W and the conventional melting range shown in I are approxi- 
mately 0.001° C. 

Figs. 8 and 9 show examples of linear and logarithmic inverse plots from 
which estimates of the widths at half-maximum and the widths at the 1% 
level of the melting ranges of all the samples were made; Table VI summarizes 
the measurements of the three melting parameters. 
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Fic. 9 Inverse melting curves of samples C.1, B.1, and 5.7: 1, after very fast freezes; 
and U1, after induced freezes. 


TABLE VI 


ALLOY MELTING RANGES AFTER VERY FAST FREEZES (1) AND INDUCED FREEZES (I1) 


From inverse curves (° C.) 


= Total 

From direct curves (° C.) Width at } max. Width at 1% level Estimated impurity 

- range trom content 

Sample I Il I II I II Table IV (p.p.m.) 

eo. 

R.1 0.001 0.001 0.0011 0.0011 0.0045 0.0041 <0.0013 <0O.5 

R.2 0.001 0.001 0.0010 0.0007 0.0029 0.0032 <0.0018 <0.5 
S.4 0.001 0.001 0.0008 0.0008 0.0053 0.0050 <0.00384 <2 
5.7 0.005 0.007 0.0014 0.0014 0.015 0.015 0.015 6 
C3 0.01 0.03, 0.005* 0.0033 0.0009 0.019 0.0386 <0. 0082 <4 
C2 0 .02-0.038 0.10 0.031 6.045 0.09 0.15 0.10 5) 
B.1 0.01 0.038, 0.005* 0.0027 0.0018 0.020 0.039 <1 
cae 


B.2 0.01-0.02 0.04 0.003 0.004 0.020 0.019 <0.23 <70 


*Melting range of the high purity cores in these samples after induced freezing. 


Plateau Temperatures 

The second part of the investigation involved a careful comparison of the 
plateau freezing temperatures developed on slow induced freezes of S.4 with 
the plateau temperatures on similar freezes on the seven other alloys. Two 
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furnaces were employed in order that sample S.4 and the sample under com- 
parison could be allowed to undergo slow induced freezes simultaneously, 
enabling the measurements of the plateau temperatures to be made inside 1 
hour on both samples. The measuring sequence was as follows: 


Br:, Bz, F.P.. Z2me (S.4), Zine CX. Be, Ba. FP. 


p.m. Zine (%),. Zine (S.4), Be., Ba. T.P. 
Br.— Bridge Ratio; Bz. Triple Point of Water. 


a.m. 


Bridge Zero; T.P. 


Measurements were made with the mid-point of the sensing coil of a single 
Meyers thermometer immersed about 15 cm. in the zine samples. A single 
triple-point cell was used throughout the study: one ice mantle was prepared 
and maintained during the work on R.1 and R.2 and another mantle was 
used for the comparisons on the other five samples. 

Fig. 10 gives results of comparing the plateau temperatures of R.1 and R.2 
with those obtained for S.4 and involves 40 slow induced freezes. Fig. 11 gives 
similar results on the comparisons of C.1, B.1, and 5.7 with S.4, where five 
determinations on each sample were made. Only two precision runs were made 
on each of C.2 and B.2 because the precise location of their liquidus points 
does not appear important in temperature-standard considerations. The inter- 
comparisons are summarized in Table VII; and in addition the final column 
in the table shows the estimated lowerings of the liquidus temperatures from 
the pure element derived from the chemical analyses of the samples and the 
known phase diagrams. 


TABLE VII 


INTERCOMPARISON OF PLATEAU TEMPERATURES 


Plateau Liquidus 


Reproducibility 
temp. rel. depression 


Ratio Rzpn/Rrp 


(equivalent 10—#° C.) 


Mean No. of 108 - to (° C.) from 
date Sample freezes Mean Std. Dev. Std. Dev. Spread S.4(°C.) Table 1V 
18.7.57 S.4, R.l 20 2. 56806227 +60 +1.7 6.7 
R.1 10 6194 +48 +1.4 1.5 0.0001 <0.00014 
S.4 10 6259 +55 +1 ti $.7 0.0 <0.00048 
25.7 .57 S.4, R.2 20 2 SH806359 $ti5 +1.9 5.5 
R.2 10 6374 +60 +1.7 4 +0. 00009 <0. 00011 
S.4 10 63438 +70 +2 .0 §.5 
21.7.5% SAR, R2 10 2. 56806293 +90 +2.t6 9.0 
10° 
Mean Dev. Mean De 
10.8.57 S.4 ) 2. 56806120 + 32 +0.9 3.7 
Cl 5 5998 + 43 +1.2 +.6 0.00035 <0.0011 
14.857 S.4 5 2. 56806155 t 45 +1.6 6.0 
B.1 5 5966 + 19 +0.7 2.¢ 0.00054 = <0.00019 
18.8.57 tf 5 2.56806150 + 46 +1. 1.3 
§.7 5 5582 + 19 +0.7 2.4 0.0016 0.0019 
22.8.57 S4 2 2.56806158 + 58 +1.7 3.3 
eZ 2 3262 +122 +34 7.0 0.0083 0.0075 
24.8 .57 S.4* 2 2. 56806142 + 44 $1.3 2.5 
B.2 2 + 44 +1.3 2.5 —0.0102 <0.016 


3578 


*The average Rz,/Rtp for S.4 for the 39 slow induced freezes was 2.56806224 with a standard 


deviation equivalent to +0.0002,° C. and a spread of 0.0011° C.; these numbers reflect the 
stability of the whole measuring system including thermometer, triple-point cell, and bridge. 





596 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


58 


~0.001°C 


2 
~444KK blu 


65.444xxX b.u. 
RzNn 





48 
3 3 
2 = @ op 28 @ = 
a x 24 + : i Sie igo 8 aiOne 8 ‘iz 7 26 x a 
as ~0,0002°C, hina ~O.0tae c ee 
+ 22 a 24 *. 
g s 
* é 
. 4 + = i oer = 
ao | 8 en 
— 
<5 ~o.oorc*® , ° BON a a 6 ~0.001°C 
z9 . @ © o 8 | 
eu | ' ‘. 






— o 

(i acess ciatiicanaiaiaaiiammea emma aa ain 
as ~0.0002°C 

® ® “-— 





Rs-4-Rx 
pohms 
°o 








3 
3 
2 
N v 
= 3 
o 
° 
3 8 
2 30 - 8 8 
a x ' e°e GS 6 
ex ~ 0.0003 °C ee 
a) | eae 8 ® 
= °° 8 
pee + = 
a, 62 Te Amie, ooo is oe 
2 = 
& & Qo—-O—_ 8 0~9—0~5~4 
Qo 
=? ~0.002°C 
xn 9 
x” 56 | 


Rs-4-Rx 
9 
pohms 





CONSECUTIVE INDUCED FREEZES 


Fic. 10. Intercomparison of plateau temperatures of samples R.1 and R.2 with S.4: Rrp 
and Rzn, zero-current fully corrected resistances of the thermometer at the triple point of 
water and zinc points respectively. 

Fic. 11. Intercomparison of plateau temperatures of samples C.1, B.1, and 5.7 with S.4. 
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DISCUSSION 

The plateau freezing temperatures developed on induced freezes of high 
purity zinc now appear unquestionably to be the liquidus points of the alloys. 

The three melting parameters in conjunction with varied freezing techniques 
provide a means of gauging impurity effects from melting curves. Thus samples 
C.1 and B.1 under all phases of the thermal analyses are strikingly similar 
and one could suggest that they have nearly the same impurity contents. The 
segregation of a specific impurity may account for the minor arrests on the 
melting curves of these alloys following induced freezing and for the steps on 
the melting curves of C.2 after very fast freezes. It is hoped that refinements 
in experimental techniques will allow more detailed measurements to be made 
on the melting contours on a series of special alloys and so allow us to obtain 
information on the presence of specific impurities. Parallel metallographic 
inspection would undoubtedly provide valuable correlations with both thermal 
and chemical analysis. 

As expected, the samples with the highest liquidus points have the smallest 
melting ranges (effects of copper impurity contents are overweighed by the 
other impurities); in this regard the width at the 1% level appears to be a 
more sensitive criterion of purity than the liquidus point of the alloy. 

On the basis of both plateau temperature and alloy melting range the 
Russian zines are certainly as pure as the New Jersey Super Pure Zine and 
the source of the difference in the reported zinc freezing points referred to 
earlier must lie in other aspects of precision resistance thermometry. The 
sample S.4 was first used in January, 1956, and has been used occasionally for 
thermometer calibrations since that time, being removed from the furnace 
after each calibration; the Russian samples were loaded in July, 1957. The 
agreement in the plateau temperatures of these samples is further evidence 
of the high stability of the liquidus point of high purity zinc as a precision 
temperature standard. 

The advantage of thermal analysis for the detection of soluble impurities in 
high purity zinc is that reasonably accurate relative numbers may be identified 
with liquidus points and alloy melting ranges from sample to sample. It is not 
necessary to depend on estimates of upper limits of these physical quantities 
derived from uncertain impurity analyses and even more doubtful knowledge of 
the phase diagrams at the zinc-rich end. The 1-kg. ingots required for precision 
resistance thermometry allow the liquidus-point determination to be made 
on fairly representative samples, which help to eliminate sampling errors 
frequently found in quantitative spectrochemical analysis. 

The disadvantage of precision thermal analysis is that particular impurities 
cannot be immediately identified. 


CONCLUSIONS 


No measurable difference was found between the liquidus points of two 
samples of high purity zinc from the D.I. Mendeleev Institute of Metrology, 
U.S.S.R., and the liquidus point of New Jersey S.P. zinc previously examined 


in this laboratory. 
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New information on the usefulness of precision thermal analysis for metal- 
lurgical inspection of high purity metals has been presented. 
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ON THE FRICTION OF INFLATED RUBBER TIRES ON ICE! 


C. D. NIVEN 


ABSTRACT 


The friction on ice of some small inflated rubber tires was measured on a 
turntable in a cold room. When rolling-friction force was plotted against load, the 
relation was either linear or slightly curved away from the load axis; such curva- 
ture implies that Thirion’s Law does not hold for rolling friction. On the other 
hand when sliding-friction force was plotted against load the curvature was 
toward the load axis as would be expected if Thirion’s Law applied. The co- 
efficient of friction can go as low as 0.01 or even lower for a hard-pumped tire 
when the temperature is near 0° C., but at —1°C. rolling friction on dry ice 
is quite appreciable. The results refer only to measurements at very slow speed. 


INTRODUCTION 

A considerable amount of work (Gough ef a/. 1956; Hadekel 1952; Hample 
1955; Milwitzky et al. 1955; Pike 1949; Taborex 1957) has been done on the 
friction and performance of rubber tires on surfaces other than ice. Férster 
(1956), working in Germany, has included in the various road conditions he 
tested the ‘‘glazed ice’’ condition. Laboratory work has been undertaken on 
small flat samples of rubber on smooth surfaces; the formula |/u = a+dp, where 
a and } are constants and ? is pressure, has been proposed by Thirion (1946) to 
describe rubber friction; this law has been verified over a 10,000-fold load 
range by Denny (1953). Conant, Dum, and Cox (1949) have used flat rubber 
samples on ice and have concluded that the base polymer and the type and 
amount of carbon black were important; Pfalzner (1950) on the other hand 
has stressed the hardness of the rubber. 

Little or no laboratory work has been published, however, up to this date, 
on inflated rubber tire samples on ice. This problem is quite exceptional from 
the theoretical standpoint because it involves perhaps the two most exceptional 
of all common substances in the field of solid friction—with the added compli- 
cation that the soft substance has a gaseous substratum. The work described 
below is an approach to this unique friction problem and has been undertaken 
with a view to getting a better understanding of the ice friction mechanism 
rather than to improving rubber formulations. Many of the observations in 
the cold room confirm the statements and beliefs of winter drivers. One of these 
statements which is of special interest claims that “heavy vehicles grip an icy 
road better than light ones until they start sliding and then they slide worse’. 
The measurements point to the reason why drivers have been led to this 


conclusion. 


METHODS AND RESULTS 


The turntable described in a previous communication (Niven 1954) was 

‘Manuscript received November 19, 1957. . 20-8 ; 
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used for the work. The straight loading arm had, however, to be changed to an 
arm with an arch above the track (see Fig. 1) in order to accommodate wheels 
up to 8 in. in diameter. Furthermore, variation in wheel diameter necessitated 
means to raise and lower the universal at the center of the turntable, where the 
arm was attached. 

In the measurement of the limiting rolling friction the table was gently 
moved by hand, while a torque was applied to the test wheel resting on the 
smooth surface of ice. The torque which just caused the wheel to slip gave a 
measure of the rolling-friction force, provided that the radius of the wheel was 
known; it was necessary then to ascertain the deflection of the tire. Some 
measurements were also made on the sliding friction with the test wheel 
locked. 

The variation in cold-room temperature was about +0.5° C. Since the 
main interest centers around temperatures at or just below 0° C., very little 
work was done at temperatures much below —3° C. For measurements near 
the freezing point, experience showed that it was wise to keep the cold-room 
temperature below 0° C. in order to avoid the danger of melting the surface 
and of thereby increasing the friction. This temperature has been shown on the 
graphs in Figs. 4 and 5 as —1° to 0° C. In Fig. 6, however, the temperature is 
shown as 0° C. because it was much nearer to the melting point of ice, on one 
occasion moisture on the track actually being noted. 

The samples used for the work consisted of three model aeroplane tires 
6 in., 4 in., and 2 in. in diameter and of two small standard tires 8 in. in dia- 
meter referred to below as A and B. A wasseveral years old; B was new from the 
factory. The model tires were tubeless and were not reinforced with canvas; 
the small standard tires had an inner tube and a canvas-reinforced cover and 
could be inflated to 30 p.s.i. The model tires could not be inflated to hard 
pressures and their valve design provided no means for measuring air pressure. 
When fully inflated they were only as hard as a very softly pumped bicycle tire 
and could easily be pressed in with the finger and thumb. The deflection and 
contact areas under load, shown in Figs. 2 and 3, give an idea of their condition 
under test. Their widths were ca. 2$$ in., 13$ in., and 3$ in. respectively. The 
tread was carefully buffed off. When this was done on the 4 in. tire, there 
appeared to be a ridge on the inside so that it was slightly different from 
the other two samples. The rubber in the model tires was a compound made 
up from natural rubber. The hardness of the inflated tire measured, on a Shore 
Durometer Type A, 26 for the 6 in. and 4 in. models and 33 for the 2 in. model. 
These figures take into account the cushioning effect of inflation. 

The two standard tires had a width of about 2} in. They are known in the 
trade as 2.50-4 Industrial-type tires (in the trade tires are described by the 
size of the cross section rather than by the diameter). The tread, which 
incidentally was not buffed off, consisted of four thick ridges about } in. wide, 
spaced about } in. apart; it was composed, in the case of the A tire, of ‘‘hot 
GRS butadiene-styrene synthetic with a black reinforcement of channel 
black and high modulus furnace black’’, and, in the case of the B sample, of 
“synthetic rubber of the styrene—butadiene type and a black of one of the 
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Loading arm showing arch and adjustable universal joint. 
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Fic. 2. Deflection vs. load. O 6 in. model, @ 4 in. model, A 2 in. model, O 8 in. synthetic 
A at zero inflation pressure, @ 8 in. synthetics A and B at 30 p.s.i. inflation pressure. 
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Fic. 3. Apparent contact area vs. load. 0 6 in. model, ™ 4 in. model, A 2 in. model, 
O 8 in. synthetic A zero inflation pressure, © 8 in. synthetic A 30 p.s.i. inflation pressure, 
Q@8 in. synthetic B 30 p.s.i. inflation pressure. 
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new furnace types’. Sample A had a very strong wall and was compressed 
little even under a load of 44 kg. with zero inflation pressure; sample B was 
not nearly so strong. The hardness of sample A was 69 and of B 62 on the Shore 
Durometer. The cushioning effect of inflation did not arise, as these reinforced 
samples were so strong. 

Apparent contact areas were approximately determined by pressing the 
tire against a rubber-stamp ink-pad and then letting it rest on a piece of white 
paper under load. 

Obviously a great many variables were present—temperature, loading, 
diameter of tire, kind of tire, quality of rubber, ice conditions. Over and above 
these variables there were two kinds of friction to deal with, depending on 
whether the wheel was rolling or sliding. It was soon found that the friction 
of a rolling wheel was much more easily defined in words than measured in 
practice, because under certain conditions when the torque was applied there 
was a slow creep, as though the wheel were turning in a very viscous fluid. 
The experimenter has to decide what torque really corresponds to the limiting 
friction value. The writer took the torque corresponding to a very slow creep. 
Under other conditions, when the wheel was subjected to torque, it might 
suddenly ‘‘let go”’, just at the limiting value, and this value would then depend 
on how gently the turntable was rocked. The writer took the value corre- 
sponding to a very gentle movement. 

After the wheel had slipped by application of too large a torque and repetition 
of the measurement was being made at lower torque, far lower loads on the 
weight pan would cause slipping than those which had failed to cause slipping 
only a few moments earlier. It was found, however, that by rocking the turn- 
table without application of torque for a little while—say 20 seconds— the 
wheel would again take its original grip on the ice. The higher torque value was 
assumed to be the more correct one. 

In the measurement of sliding friction of a material like rubber, which has 
usually a very appreciable friction, high turntable speeds cannot be used for 
dry ice tests—at any rate if the sample runs in the same track all the time— 
because water accumulates on the track and the results have little meaning. 
The sliding-friction values reported below refer, accordingly, to very slow 
or ‘‘just-stopping”’ speeds. It is important to emphasize this so that the reader 
will not apply the values reported to the practical problem of braking at speeds 
of 30 to 100 m.p.h. 

The low-pressure tires could be deflected with load until they were almost 
“running on the rim’’. Such deflections were not feasible for the standard 
tires with the equipment available; they were designed to carry about 130 kg. 

When friction force was plotted against load, the points for a series of 
measurements at different loads lay as a rule surprisingly close to a smooth 
curve or straight line in spite of the troubles in assessing limiting torque, but 
there was some lack of repetition due to the ice. That this was to blame there 
was no doubt, because two series of points would give similar curves, but all the 
points of one series would be either lower or higher than all the points of another 
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series. Instances of this can be seen in Fig. 6. The trouble lay not in the method 
of measuring the friction but in the conditions at the interface. 

The ice surfaces were prepared by machining with a wide-edged steel tool 
without employing any special buffing technique afterwards. No attempt was 
made to separate the data on freshly machined ice from those on ice which 
had stood overnight because this would not have cleared up the trouble. 
Buffing the surface by having the 6 in. model tire slide on the revolving 
turntable at high speed, until the surface became as smooth as glass, could 
actually cause an increase in friction (except near 0° C. when the buffed surface 
had if anything the lower value). Presumably buffing permitted closer contact 
and increased adhesion. Decrease in friction could sometimes be related to the 
possible presence of ice crystals—so to speak the ‘‘turnings’’—which had not 
been thoroughly enough removed after machining. 
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Fic. 4. Friction force vs. load curves for model tires. A circle signifies rolling friction, a 
triangle sliding friction, and a cross on top of a symbol water on the ice. 
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A true assessment of ice friction values cannot be reached by using one type 
of specially prepared surface. One must accept the fact that various interface 
conditions are possibilities. Research might conceivably succeed in attaching 
coefficients to specified conditions but the fact must be faced, nonetheless, 
that under practical conditions variation in friction coefficient occurs for 
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Fic. 5. Friction force vs. load curves for 8 in. (diameter) synthetic standard tires. A 
circle signifies rolling friction, a triangle sliding friction, and a cross on top of a symbol water 
on the ice. 
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apparently very little reason. Fortunately these variations in interface condi- 
tions do not generally affect the type of friction-force-versus-load curve 
obtained at a particular temperature. 

Results on the model tires are shown in Fig. 4 and on the standard tires in 
Fig. 5. Fig. 6 shows some special curves relating to the A tire at zero inflation 
pressure at a temperature very near the melting point of ice, which will be 
discussed below. 

In Figs. 4, 5, and 6a circle signifies rolling friction, a triangle sliding friction, 
and a cross on top of a symbol water on the ice. Broken lines have been used 
in connecting points of one particular series of readings; when two or more 
series of readings have been combined in drawing a curve a full or unbroken 
line has been used. Since the load scale is one tenth of the friction force scale, 
coefficients can be read on the 10 kg. ordinate by merely placing a straight 
edge on the graph and by changing a decimal point. 
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Fic. 6. Friction force vs. load curves for 8 in. (diameter) synthetic standard tire showing 
strong curvature of rolling-friction curves toward load axis at 0° C. and approximate con- 
stancy of friction force at high loads. Circle, dry ice; cross on top of circle, wet ice. 
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Fic. 7. Coefficients of friction at —3° C. for 6 in. model tire and 8 in. (diameter) synthetic 
B standard tire. 
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In order to assist the reader in visualizing the type of value to expect for the 
friction coefficient of a rubber tire on ice which is well frozen but not extremely 
cold, Fig. 7 has been constructed from appropriate data in Figs. 4 and 5. The 
four lines in Fig. 7 indicate respectively rolling and sliding friction coefficients 
at —3° C. for the 6 in. (diameter) soft-pumped model tire and the 8 in. (dia- 
meter) hard-pumped synthetic B standard tire; the coefficients apply to very 
slow speeds as stated already. 

Table I gives some further details on the values of friction coefficients—also 
taken from the graphs. Some of these values have been read to the second 
significant figure but variations in ice conditions do not really warrant this 
accuracy. All values are approximate and the object of tabulating them is to 
demonstrate the wide variety of friction coefficients, each of which is a perfectly 
good value for a particular tire at a particular temperature on a particular ice 


surface. 
rABLE |! 


Friction coefficients 


Ice Sliding 
temperature, Load, (very slow 
Tire Ee kg. Rolling speed ) Remarks 
2 in. model Oto —1 0.8 0.35 0.12 
—3 0.8 0.45 0.31 
4 in. model Oto -1 1.0 0.37 0.11 
—3 0.8 0.4 0.36 
6 in. model Oto -1 0.8 0.35 0.14 to 0.37 — Including wet ice 
2.5 0.5 0.12 to 0.3 Including wet ice 
—3 0.8 0.33 0.2 to0.5 
2.5 0.5 0.3 to0.4 
—t 0.8 0.5 0.35 to 0.5 
8 in. standard 0 14 0.01 Without pouring water on ice 
(synthetic) Oto —1 44. 0.08 to0.2 0.04 to 0.09 Dry ice 
44 0.08 0.03 Wet ice 
—2to -—3 44 0.23 0.14 Dry ice 
—2 14 0.03 Wet ice 


—6 44 0.25 0.11 to0.16 Dry ice 
DISCUSSION OF RESULTS 

When a tire is deflected, the volume gets slightly smaller and the inflation 
pressure rises. If, however, it could be arranged by means of some device to 
keep the inflation pressure the same for all loads, and if the strength of the 
wall were negligible, the value of » in Thirion’s equation 1/u = a+bp would be 
constant, and since a and 6 are numerical constants depending on the sample, 
un would remain constant; the friction-force-versus-load relation would then be 
linear. The results on the model tires show, however, that this is only approxi- 
mately true: actually p increases with load and so the sliding-friction-force- 
versus-load curve should bend towards the load axis at high loadings—and 
this is what it does. One may, therefore, conclude that the sliding-friction 
results on the model tires do not conflict with Denny’s conclusion that Thirion’s 
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Law is widely applicable. Fig. 8 supports this conclusion: the straight line 
represents the equation 1/y = 1.85+2.7p; the points are reciprocals of 
sliding-friction values taken from the graph in Fig. 4 corresponding to the 


6 in. model tire at —3° C. 


FRICTION COEFFICIENT 
nm 


RECIPROCAL OF SLIDING 











0.20 0.25 0.30 
PRESSURE -Kg./sq. cm. 


Fig. 8. Verification of the applicability of Thirion’s Law to sliding friction by using values 
as given by the graph in Fig. 4 referring to the 6 in. inflated model tire at —3°C. 


However, what is of much greater importance is the fact that the rolling- 
friction results can never be expected to fit Thirion’s equation because the 
curvature—such as there is—is away from the load axis at high loadings: 
in other words as p increases, 1/u gets smaller, a condition which cannot be 
made to fit the equation if } is positive. 

The application of Thirion’s Law to the 8 in. wheels cannot be discussed 
because the accuracy of the contact-area measurements was not good enough. 

On the graph for the 6 in. model tire at —1° to 0° C. (see Fig. 4) the two 
‘“‘broken-line’”’ curves show widely differing results on sliding friction for dry 
ice: one of them is near the curve for wet ice and the other near the curve 
shown on an adjacent graph for dry ice at —3° C. A small change in tempera- 
ture is believed to have caused this, since evidence which will be discussed 
below in connection with Fig. 6 indicates that the friction properties of ice 
change markedly somewhere between 0° and —1° C. 

At —3° C., the sliding and rolling frictions of each of the model tires were 
comparable. A series of low points on the graph for the 6 in. model tire at 
—3° C. may have been due to tiny ice crystals on the surface, because that 
particular series of readings happened to have been taken on newly machined 
ice. Sliding friction can be higher than rolling friction under some circumstances. 
This could be verified visually if it were assumed that slow rolling friction is 
static friction (Taborex 1957), because the force necessary to balance the 
sliding friction when the table was turning slowly would draw the arm back as 
soon as the table stopped. This small detail suggests a peculiarity about rubber 
possibly connected with the fact that the braking coefficient of a tire on cement 
rises at first, quite sharply, with small percentage of slip before falling off at 
greater percentages (Férster 1956; Milwitzky et al. 1955). 

Fig. 5 shows that the friction of the standard tires was lower than that of 
the model tires. The hardness of the standard tires was partly responsible for 
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this; in addition they were made of synthetic rubber, while the model tires 
were made of natural rubber. The two synthetic tires had comparable friction 
values except near 0° C. but since small changes in temperature have a large 
effect near that temperature according to other evidence, this difference need 
not be stressed. For the hard synthetic tires, rolling friction obeyed very nearly 
a straight-line relation while the sliding friction did not, the curves bending 
toward the load axis. 

The results shown in Fig. 6 are entirely different from those in Fig. 5 and it 
is the difference in temperature which caused this. (The tire was so strong and 
hard that it deflected little even under a load of 44 kg. either at 0 or at 30 
p.s.i. inflation pressure.) The results in Fig. 6 were obtained before it was 
decided to work always at a temperature slightly below 0° C., and they refer 
to ice which was almost melting. 

It was only after the experimental work had been discontinued and the 
apparatus disconnected and removed from the cold room that the significance 
of these results was appreciated: they show that something happens between 
0° and —1° C. which is of far greater consequence as far as friction is concerned 
than anything which happens between —1° and —6° C. The explanation of 
this result is almost self-evident if one accepts the modern theory that the 
surface of a solid is covered with asperities and that friction is due to adhesion 
to these asperities. In this instance the ice is to be regarded as covered with 
asperities composed of water molecules held together by strong molecular 
forces, so long as the temperature is below 0° C. As soon as the temperature 
rises above 0° C. these forces let go and the rubber is then adhering to some- 
thing which has ceased to be solid; one would accordingly expect the friction 
to drop. In other words, since the freezing point of water is sharp, one would 
expect a discontinuity in the value of the friction at or very near —0° C. 

The results of the investigation under the controlled conditions of the cold 
room can be applied to confirm and explain driving experience on ice. The 
coefficient of rolling friction is either constant (Amontons’ Law) or it increases 
at high loading for temperatures definitely below 0° C. The coefficient of 
sliding friction on the other hand decreases with load particularly near 0° C. 
These facts explain why some drivers believe that heavy vehicles have more 
friction than light ones until sliding starts, because rolling friction force is 
either proportional to or better than proportional to the load on the tire, and 
so, if the load is unevenly distributed over the wheels, this force will appear 
small on the driven wheels when they carry less than their share of the load—a 
condition often found with light vehicles—and large, if they carry more than 
their share of the load—a condition often found with loaded trucks. There is at 
present insufficient evidence to say that the rolling-friction coefficient of an 
automobile tire on ice increases with load although the results on the model 
tires suggest this as a remote possibility. These tires were, however, deflected 
to a far greater extent than would be tolerated for the canvas-reinforced 
automobile tire. When the ice is exactly at 0° C. the work indicates that a 
heavily loaded hard tire should grip the road even less in proportion than a 
lightly loaded one, either rolling or sliding. 

Low-pressure natural-rubber tires seem to grip well on ice. The rolling friction 
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of rubber tires on ice—even on wet ice—at a temperature as mild as — 1° C. is 
by no means negligible. When friction is low, torque should be applied gently 
to avoid spinning as this causes a temporary reduction of friction; gentle 
braking is likewise desirable. The really difficult driving condition appears to 
arise when the ice is at exactly 0° C., and because it is critical to have the ice 
temperature at 0° C., impossible driving conditions occur relatively rarely. 
The overruling danger factor in driving on ice is not so much that the friction 
is so low as that it is exceedingly variable, changing unexpectedly, yet con- 
siderably, because of very slight change in temperature near 0° C. or because 
water or loose ice particles happen to be on the surface or because the surface 
is not of uniform smoothness. 

The important conclusion to draw from Table | is that 0.1 is not a sate 
coefficient to take for the rolling friction of a rubber tire on ice if the worst 
weather conditions are to be considered; 0.01 is nearer to the figure for safety 
in engineering design, at any rate for a hard tire. However, the hardness of the 
tire, the load, and the apparent contact area are such very important factors 
that until full-scale tests are made on an ice surface controlled at close to 
0° C.—much like the tests made on a concrete surface at the Langley 
impact basin (Milwitzky, Linguist, and Potter 1955)—a precise minimum 
cannot be placed on the rolling friction for an inflated rubber tire. The time is 
well Overdue to provide facilities for such tests in Canada because a small 
turntable cannot give all the answers; experience on the turntable has shown, 
however, that large-scale tests would be of little value unless the runway and 
ambient air temperature could be controlled to a temperature extremely close 
to 0°C.; at —1° C. the runway would no longer be at the danger temperature. 
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THE SOMMERFELD MODEL AND THE HEAT OF SOLUTION 
OF AN ALLOY! 


W. G. HENRY 


ABSTRACT 


An approximately self-consistent method, using a Sommerfeld model, and 
incorporating an approximation to the second-order perturbation energy 
(Lennard-Jones 1930), is used to calculate the heat of solution of both monovalent 
monovalent and monovalent—polyvalent alloy systems. The heat of solution, as 
given by the second-order approximation, varies directly with the square of 
the perturbing potential, and inversely with the mean total energy of the free 
electrons in the solvent. Varley’s (1954) result resembles this, except that his 
expression varies inversely with the mean kinetic energy of the free electrons in 
the pure solute. The first-order approximation to the energy of the alloy system 
is identical, except for a change of sign, with the expression suggested by Friedel 
(1952, 1954) for the heat of solution. The charge accumulation about a solute 
site is found to vary directly with the strength of the perturbation, and inversely 
with the mean energy of the free electrons in the solvent. Varley obtained a 
comparable expression; however, his function varies inversely with the mean 
kinetic energy of the free electrons in the pure solute. The calculated charge 
accumulations in the copper-silver system, and for gold dissolved in silver, agree 
qualitatively with those of Arafa (1949) and Huang (1948). 


1. INTRODUCTION 

Oriani (1956) made a general survey of the calculations of the heats of 
solution of alloys. The calculations which are relevant are those of Huang 
(1948), Arafa (1949), Friedel (1952, 1954a, 6), Varley (1954), and Shinohara 
(1955). 

Shinohara (1955) attempted a self-consistent calculation of the heat of 
solution of a polyvalent solute in a monovalent solvent. It was assumed that 
the perturbing potential due to a solute atom was of the screened Coulomb 
type, and the perturbation calculation was carried out to second order. The 
fact that the calculated heats of solution were 30 to 40 times too large was 
ascribed to a deficiency in the method. Varley (1954), who used a ‘“‘two-band”’ 
model, treated both monovalent-monovalent and monovalent—polyvalent 
systems. The form of the expression for the heat of solution obtained by 
Varley (1954) is similar to that of the present calculation. Friedel (1952, 
1954a, 6), on the basis of a cyclic process, derived an expression for the heat 
of solution of a polyvalent solute in a monovalent solvent which, on the 
basis of the present calculation, appears to be not related to the actual heat 
of solution. Arafa (1949) applied a cellular method to the copper-silver 
system, and obtained agreement with the heats of solution determined from 
the phase diagram. This method has not been applied to a monovalent- 
polyvalent system. 

Huang (1948) assumed a square-hole potential, and treated the case of gold 
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dissolved in silver by a self-consistent Thomas—Fermi method. After applying 
a large correction to take account of the kinetic energy associated with the 
rapid change of potential of the potential hole, he obtained good agreement 
with experiment. In the same paper he made an exact, non-self-consistent, 
quantum-mechanical calculation for the square-hole potential, and obtained 
—0.69 ev. per atom of gold for the heat of solution. The present approximately 
self-consistent calculation is an approximation to the exact square-hole 
potential problem, and is of more general applicability. It is a perturbation 
calculation which is carried out to second order, and in that respect it is 
similar to that of Shinohara. The method of carrying out the calculation, 
however, is quite different, in that an approximation to the second-order 
energy (Lennard-Jones 1930) is used in contrast to Shinohara’s attempt to 
evaluate exactly the second-order terms. 

The heats of solution of binary metal systems show the following general 
characteristics. The heats of solution of monovalent—monovalent systems and 
monovalent—polyvalent systems are of the same order of magnitude. There 
are no apparent systematic valence effects. The heat of solution is negative 
when either B is dissolved in A or vice versa. The present calculation attempts 
to relate the general features of the heat of solution to a simple model. The 
accumulation of charge at a solute site is also discussed. 

In Sections 2 and 3 the non-self-consistent calculation is given. In Section 4 
the self-consistent calculation is discussed. In Section 5 the problem of the 
polyvalent solute is discussed. 


2. THE FIRST-ORDER APPROXIMATION 
The integral heat of solution AH for N atoms of solution, where JN is 
Avogadro’s number, is defined as the difference in the heat content of the 
solution and of the components from which it is made. The integral heat of 
solution for small concentrations may be represented by the first term of a 
series 


(2.1) AH = AH'(0)n/N, 


where » is the number of solute atoms and AH’(0) is the slope of the integral 
heat of solution curve at zero concentration, which is equal to the heat of 
solution per gram-atom of solute at zero concentration. In the following 
sections, the difference between the heat content of a solution formed from 
N-—1 atoms of species A and one atom of B, and the heat content of the 
unalloyed components, is calculated. The difference according to (2.1) is 
equal to AH’(0)/N. 

In the Sommerfeld (1928) model of a one-dimensional pure metal the 
unperturbed wave functions ¥,;° are 


, , 1 . 
(2.2) v.° = (2/L)? sin(krx/L), 
and the energy associated with a level is given by the relation 


. h’k? 
9 PO, ARIE eee 
(2.3) W, J +207: 
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L is the length of the metal, 4°k?/8mL* is the kinetic energy of the kth 
level, and V° is the Sommerfeld potential which is independent of x. The 
total energy of the system of electrons W° is given by the expression 


Wek? 


9 \ PD ee To 5 
(2.4) Wi ii +253]: 


where the summation ts over all occupied levels. The zero level of energy is 
taken to be the energy of all ions and electrons at infinite separation. 

If it is assumed that one solute atom of species B, of the same valence and 
atomic length as the solvent A, may be represented by a perturbation (Fig. 1) 
of length 6 = L/N and of strength 1’,° equal to the difference in the Sommer- 
feld potentials of the free electrons in the pure components, then the change 
in energy of the system may be calculated by perturbation methods. The 
limit of the applicability of the method is discussed in Appendix 1. Since the 
energy levels of the one-dimensional unperturbed system are non-degenerate, 
the perturbation theory for non-degenerate systems may be applied. In three 
dimensions, it is assumed that the neglect of any near-degeneracy is not 


serious. 





Fic. 1. The perturbed potential in a one-dimensional metal of length ZL. 


The first-order approximation to W,, the energy of the &th level, is 
(Henry and Rogers 1956) 


Me Bie 
W,'” = W,° a f sin? #2 Vy dx 


= We+V,° tps sin( ne #28) |. 


By considering the levels in pairs, it may be shown (Appendix 2) that the 


to 
ur 


mean energy of a level is changed by an amount V’,°6/L. The result in three 
dimensions, which is discussed in Appendix 2, is V,°6,6,6,/L,L,L:. 
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For a binary system of monovalent metals with no volume change, the 
first approximation to the total energy of the perturbed system is 


ED. caer eS ryo 6 WR* 
(2.6) LW = NVSTNV SE T+ ga 
= (N—1)Va°+Va°+NG,°, 

since 6 = L/N and V,° = Vpg°—V,°, the difference in the Sommerfeld 


potentials of pure A and B. &,° is the mean kinetic energy of an electron in 
the solvent. The right-hand side of (2.6) is, however, identical with the total 
energy W, of the pure components. Thus, the difference in heat content 
between the solution and the components from which it is made is zero, 
according to this approximation. That is, the first-order perturbation gives 
only the energy of a phase mixture of the pure components, and indicates no 
term favoring the formation of a solution. 

For a binary system of monovalent metals with a volume change, there is 
an additional negative contribution to the heat of solution arising from the 
change in the kinetic energy. The total kinetic energy of the solution con- 
taining one solute atom is 


‘ hk? no 2ANG,° (3) ve 
2.47) 5 = VE,P-- = +31 —~] NG°+ ..., 
oo x 8m(La+A)~ Ca Lx " Es Gat 
where A = (Ly—L,)/N. The mean kinetic energy of an electron in the pure 
solute is 


' I h°k® ..¢ 26Nne,? (2°) : 
2.8 : = 5 =," — a § ee 
: \ » 8m(VA+L,) a Isp + Ls “= 


e Kinetic energy O 1e phase xture ©, 1S eretore 
Che kinet gy of th h mixture ©, is therefor 


. ; ; .6 2ANne,° nar 
(2.9) Oy = (N-1)G4°+ Eo anes +3(4 ) ee os 3c 
“A “A 


The difference between the kinetic energy of the solution and the phase 
mixture is approximately —3(L_y—L,)?G4°/L,a?. This is also valid for the 
three-dimensional model, where L, and Lz are now the cube edges of a gram- 
atom of solvent and solute respectively. If, for example, Lg and Ly differ 
by 10%, and &,° is 3.3 ev., then the contribution to the heat of solution is 
0.1 ev. If the length of the perturbation is assumed to be equal to 1/N of 
the length of the solution, then there is no additional contribution from this 
source. In the remainder of the paper it will be assumed that the volume of 


the solute equals the volume of the solvent. 
3. THE SECOND-ORDER APPROXIMATION 
The second-order approximation to the energy of a level W,° is given by 


the expression 


eat? W.? = W.4 S, 
l Wo-W, 
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where the dash indicates that the sum is over all values except / = k. H’ is 
the formal representation of the perturbation. The sum may be rearranged 
(Lennard-Jones 1930; Pauling and Wilson 1935), to permit an approximate 
value to be found, 
, , ( 2 ro , , 

(3.2) > An, _ (HT je Thy a x We HuHer 

oe W.-W? We W,.° W.°-W? 
The approximation consists in omitting the sum. The approximation is dis- 
cussed in Appendix 3. 

For a binary system of monovalent metals with no volume change, the 


J 


first term on the right-hand side of (3.2) becomes 


(Hx 


We we ae 


ne = Vy Ee ae = | 
Celt ek” ~~ UL 


A comparison with (2.5) shows that the mean energy of a level is lowered by 
an additional amount V,°°6/W,°L. The total energy change of the system 
resulting from this term is 
(3.4) 7, pet 

W.° "Lee 
where the sum is over all occupied states. N is the total number of electrons 
per gram-atom of alloy; V°+€° is the mean energy of an electron in the 


+ Eee V2? 


Sik Va ae 


solvent. 
The second term may be evaluated and vields 


(Hi,)° L|2 (. 2kax 1,6 -f 
Wwe = W,° rm sin i ve dx 


(3.5) - a(ey 
~~ WEF * 


This term is small in comparison with the term which contains 6/L to the 
first power. 

The total energy of the perturbed system is given by the second approxi- 
mation to be 


(2) ao * 6 
(3.6) > Ww,” = NW,°+NV,° “+e ra? 
where 4° = V4°+G,°. The first two terms on the right-hand side of (3.6) 


are, by (2.6), the energy of the phase mixture. Therefore the heat of solution 
per atom at zero concentration is 


(3.7) AH'(0)/N = V,° 


wn 


W,°. 
It may be concluded from the form of (3.7) that: (1) AHZ/’(0)/N is negative 
when either B is dissolved in A or vice versa; (2) the magnitude of A/7’(0) 
depends on whether B is dissolved in A or vice versa. 
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The result expressed in (3.7) may be compared with an exact calculation. 
Huang (1948) made an exact calculation of the square-hole potential problem 
using numerical methods for the single case of gold dissolved in silver. He 
obtained —0.69 ev. per atom for the heat of solution. For monovalent metals 
the Sommerfeld potential is approximately equal to the negative value of 
the sum of the absolute values of the sublimation energy A, the first ionization 
energy J“), and the mean kinetic energy ©. If the atomic volumes are equal, 
the difference in the Sommerfeld potentials is also equal to the difference in 
the mean energies of the electrons. The atomic volumes of silver and gold 
are approximately equal. The mean energies of the electrons in gold and 
silver, as determined from their heats of sublimation and first ionization 
energies, are —13.l ev. and —10.5 ev. respectively. Equation (3.7) yields 


(3.8) AH'(0)/N = (—2.6)2/(—10.5) = —0.64 ev., 


which is in reasonable agreement with the exact calculation. 

It is of interest to compare AH’(0)/N, as given by equation (3.7), with the 
result obtained by Varley (1954) using a ‘‘two-band”’ model. For simplicity, 
the comparison will be made for a binary system of monovalent metals with 
no volume change. The change in the mean interaction energy of the electrons 
and the change in the ground state energy, which result from the transfer of 
charge, will be neglected. In the notation of the present paper, Varley’s 
result for the difference in the energy between the solid solution and the 
phase mixture A/I(n/N) is 
y O2 

3 c(l—c) V, * 

(3.9) iG) = —-_ 
“x [cE,°+ (1—c)Ep*] ’ 

where c = n/N. The heat of solution per atom, at zero concentration, is 

given by the relation 

(é * 7 O29 re 

(3.10) AH'(0)/N = —i0 V, */Gp°. 

It is thus seen that the ‘‘two-band”’ model, to which it is somewhat difficult 

to give a physical interpretation, yields a result which is very similar in 

form to that of the present calculation. 


4, SELF-CONSISTENT CALCULATION 
Before the calculation may be made approximately self-consistent, it is 
necessary to calculate the amount of extra charge which accumulates in the 
volume occupied by the solute atom. Lennard-Jones (1930) has shown that 
the wave function which corresponds to the approximation (3.2) has the 
form 


Aix H' 
( ) Vic Vi WwW, W,° 
This is also the normalized wave function if terms in 1/W,°? are neglected. 

The extra charge which accumulates in the length 6 is given by the ex- 
pression 


»L—a L—a oO V. ° 
(4.2) = | f vr dx — vera ~ yo 


a 
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This result also holds for three dimensions. As may be seen from the form 
of the wave function, the accumulated charge comes uniformly from the 
whole lattice, and each wave function contributes to the localized charge. 
Since the denominator is negative, charge is accumulated about the solute 
when V,° is negative, and is attracted from the solute to the solvent when 
V,° is positive. Arafa (1949) also concluded that the charge comes uniformly 
from the lattice. 

To make the calculation self-consistent, an attempt is made to make q 
and V, compatible. The extra charge g which accumulates is given by the 
relation 
— 2(V,°+f(q)) _ 2V, 


(4.3) q —e =—— 
Ws Ws 

where V, is the self-consistent potential. If it is assumed that the change in 
potential is equivalent to that of the charge g in a spherical shell at a radius 
ry equal to the radius of a sphere of equivalent atomic volume, then f(g) is 
equal to q/ro. 

Since the accumulated charge comes uniformly from the whole lattice, then 
it is reasonable to assume that the potential of the N—1 solvent electrons 
is lowered by an amount equal to g/(N—1)ro. The total energy of the per- 
turbed lattice is now 


i 17° (yr ° , 6 
2, We”? = NW o+N(V,°+9/ro) . 


; N(V, °+¢/r9)° 6 Bs g 
t. Ae RS ic Be ees 
4) ree Ee 
er eee oe, 
= NWe+NV, SN apts - 


Thus the heat of solution per solute atom is, as before, 
(4.5) AH’'(0)/N = V,2/W 4°, 


where V, is the self-consistent potential. 

Gold in silver will be again used as an example. V,° is equal to —2.6 ev. 
The mean energy of the electrons in silver is —10.5 ev. The potential per 
unit of charge at 1.59 A, the equivalent atomic radius, is 9.0 ev. The charge 
accumulated is therefore given by (4.3), 


(4.6) gq = 2(—2.6)/[—10.5—2(9.0)] = 0.18 unit of charge. 
Therefore, the heat of solution is 
(4.7) AH'(0)/N = (—2.6+0.18X9.0)?/(—10.5) = —0.10 ev., 
which is of the correct order of magnitude. 
The field outside the region of the solute atom which results from the 
accumulated charge has been neglected thus far. This is difficult to handle 


by the methods developed in this paper. The magnitude of this effect may 
possibly be estimated by finding the polarization energy of the 12 nearest 
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neighbors in a field of ge/r’. If it is assumed that the polarizability of the 
neighboring atoms is of the order of the polarizability of argon, for example, 
then the polarization energy —ag’e?/2r+ amounts to approximately —2.0 ¢ 
ev. This term is always negative. For gold dissolved in silver the polarization 
energy is —0.03 ev., which is less than V,? W,°. 

It is not to be expected that the results of the calculation as a whole can 
be compared directly with experiment, since at least one important term 
which may be significant in many cases, namely the change in the energy 
of both the solute and solvent ions, cannot be treated by this model. 

In terms of the self-consistent potential, the charge accumulated is given 
by 
(4.8) q= V,/W,,*. 


Agreement in form is again obtained with the result calculated by Varley 
(1954). In the notation of the present paper, his result is 


(4.9) q = —pVn /Gz - 


Arafa (1949) found a charge accumulation of 0.32 unit round a copper 
ion in a silver lattice, and a diminution of 0.29 unit round a silver ion in a 
copper lattice. The present calculation yields an accumulation of 0.11 unit 
and‘a diminution of 0.10 unit. Huang (1948) derived, by using an approxi- 
mately self-consistent Thomas—Fermi method, a charge accumulation of 
about 0.5 unit round a gold ion in a silver lattice and a diminution of charge 
in the immediate neighborhood of equal magnitude. The present calculation 
gives an accumulation of 0.18 unit. 


5. POLYVALENT SOLUTES 

A polyvalent solute has Z+ 1 electrons outside a closed shell. The discussion 
will be restricted to monovalent-polyvalent systems. In general, the heats 
of solution of monovalent—polyvalent systems and monovalent-monovalent 
systems are of the same order of magnitude, and are negative, when either 
B is dissolved in A or vice versa. A special feature of the systems with a 
polyvalent solute is the absence of clearly defined valence effects, such as 
those inferred by Hume-Rothery, Mabbot, and Channel-Evans (1934) from 
the equilibrium diagrams of the alloy systems. Kleppa (1956), who investi- 
gated systematically the heats of formation of a number of polyvalent solutes 
in copper, silver, and gold, found no definite systematic valence effect. A 
study of the collected data of Kubaschewski and Catterall (1956) confirms 
these general conclusions. 

According to the present model, the heat of solution in monovalent—mono- 
valent systems results primarily because the energy levels of the conduction 
band of the solvent are perturbed by the presence of the solute atom (see 
(2.5) and (3.3)). A polyvalent solute with Z either even or odd will also 
cause a perturbing of the conduction band, and thus there will again be a heat 
of solution arising from the same source. 

When Z = 2, it will be assumed that the Sommerfeld model may be applied 
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to the lone electron. If the self-consistent calculation were carried out with a 
well corresponding to an ion with a charge of three, the charge accumulation 
would cause the perturbing potential to be small, and the perturbation 
method could be applied. It is also plausible that, on the average, two of the 
electrons would merely be bound, that is, have almost the same average 
potential and kinetic energy as in the pure solute. At the same time, only 
the energy of one electron of the solute, which replaces the electron of the 
solvent, would be altered appreciably. The trend in the sublimation energies, 
in the two series copper, zinc, gallium and silver, cadmium, indium, suggests 
that the pair of s electrons may contribute only a small amount to the binding 
of a trivalent element, and hence their energy should not depend markedly 
on their surroundings. The sublimation energies are 3.54, 1.87, 2.5 and 3.00, 
1.16, 2.56 ev. per atom respectively. 

If the change in the energy of the two s electrons is assumed to be small, 
and only the change in energy of the N electrons in the system is considered, 
it follows by analogy with (4.4) that the energy of the system is given by 
the relation 


(5.1) > Ww. = NW,°+NV,° oN os :. 
L Wao L 
where V,° = Vp°— Va4° and J’,° is the potential energy of the electron re- 
placing the solvent electron in the pure solute. 
If the volumes of the solute and solvent metals are equal, then the heat 
of solution per solute atom is 


(5.2) AH'(0)/N = V,2/W,°, 


which is the same as for the monovalent—monovalent systems. For a monovalent 
solute in a trivalent solvent, the same assumptions lead to the result that 
the heat of solution is equal to V,? li’ ,°. The heat of solution is thus negative 
if either B is dissolved in A or vice versa. Any change in the average potential 
and kinetic energy of the s electrons will make a direct contribution to (5.2). 
’,° may be evaluated approximately (see Section 3) by treating a trivalent 
metal as a monovalent metal. It is found, for example, that for the copper- 
gallium system V’,° is approximately 2.7 ev., which is almost the same as 
for the silver-gold system. In general this should be true, and it may account 
for the fact that the heats of solution of monovalent—monovalent systems 
are similar to those of monovalent-—trivalent systems. 

There is a relationship between the first two terms of (5.1) and the expression 
suggested by Friedel (1952, 1954, 6) for the heat of solution. The first two 
terms by themselves are the first-order non-self-consistent approximation. 
The first two terms of (5.1) represent the energy of a phase mixture con- 
taining one atom of solute. Hence the change in the energy of a phase mixture, 
when an atom of solvent is removed from the system and a solute atom is 
added, is equal to V,°. If the trivalent solute is treated as a monovalent 
metal of equal volume for the purposes of calculating V,°, then it is given 


> 


by (see Section 3) the relation 
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(5.3) Y = —([|Ap|+|/p®|— Ag|—|Z,]]. 


Friedel suggested the same expression as (5.3) for the heat of solution but 
with opposite sign. 

A cycle such as Friedel uses yields correctly AH’(0)/N, provided all the 
terms are included. The two terms which Friedel assumes to have approxi- 
mately equal absolute values, namely, the energy required to remove an At 
ion from neutral A, and the energy gained on the addition of a Bt ion to 
negatively charged A forming a neutral alloy, could, in fact, differ by several 
electron volts per atom, and should not be neglected. Friedel considered that 
the change in energy of the conduction band can be neglected, but, as was 
shown in Section 3, where the second-order approximation was discussed, the 
perturbing of the conduction band resulting from the substitution of a different 
species with a different associated potential field leads to a contribution to 
the heat of solution of the order of magnitude of the observed heats of solution. 

For Z = 1, a difficulty arises because of the equivalence of the valence 
electrons. The effective nuclear charge of neutral trivalent elements (Slater 
1930) is greater than that of monovalent elements. This suggests that: (1) the 
levels of the solvent lattice will be lowered, that is, V,° is negative; (2) the 
energy of the two electrons of the solute will be raised because of increased 
screening. If it is assumed that neither electron enters the band and _ that 
the volume of the solute atom may be excluded for the purposes of calculating 





&,°, then the energy of the system is given by 
Vy 


(5.4) > Ww, = (VN-1)W4°+2( p +2(9)) + Vi +55 5 
A 


where V,° cannot be well defined in terms of the present model and g(q) is 
the change in energy per electron of the pair due to the charge accumulation. 
The heat of solution per solute atom has the form 


» 2 
« 


PS + V,°+2g(q). 


D 
Va 


la 


(5.5) Al O ae 


= 


All that can be said is that the last two terms, which are of opposite sign, 
will tend to cancel. 

In order to avoid the occurrence of large positive terms which would result 
from the forcing of electrons to occupy levels at the top of the band, it is 
necessary to assume that not all of the valence electrons of trivalent and 
divalent solutes enter the band. 


6. SUMMARY 


An attempt is made to calculate the heat of solution at zero concentration 
for monovalent-monovalent and monovalent-polyvalent systems using a 
“‘Sommerfeld-like’’ model for an alloy. The calculated heats of solution are 
of the correct order of magnitude. The model consists of a small rectangular 
well representing the solute, inside a large rectangular well representing the 
solvent. The non-self-consistent first-order perturbation calculation gives the 
energy of a phase mixture. The first-order approximation is identical in form to 
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the expression suggested by Friedel (1952, 1954 a, b) for the heat of solution 
of a polyvalent solute except for a change of sign. The non-self-consistent 
second-order perturbation problem is solved by means of an approximation 
due to Lennard-Jones (1930). The second-order calculation gives, in addition 
to the energy of a phase mixture, a negative term which is identified as the 
heat of solution. The approximation was checked and found to be in good 
agreement with an exact calculation of Huang (1948). 

According to the present model the heat of solution arises primarily because 
the energy levels of the conduction band are perturbed by the potential 
associated with the solute atom. The heat of solution for monovalent—mono- 
valent systems varies directly as the square of the difference in the Sommerfeld 
potentials of the pure components, and inversely with the mean energy of 
the solvent electrons. For monovalent—polyvalent systems the heat of solution 
varies in a like manner, although the idea of the Sommerfeld potential must 
be modified. Varley (1954), who used a “‘two-band”’ model, found an expression 
for the heat of solution which varies directly as the square of the difference 
in the Sommerfeld potentials, but inversely with the mean kinetic energy 
of the solute. 

The calculation is refined by making it approximately self-consistent. The 
main reason for making the calculation self-consistent is to allow for the 
modification of the perturbing potential which results from the movement 
of charge. The movement of charge decreases the perturbing potential and 
hence the heat of solution. It is shown that the charge accumulation about 
a solute atom varies directly as the difference in the Sommerfeld potentials 
of the pure components, and inversely with the mean energy of the solvent 
electrons. For the copper-silver system and for gold dissolved in silver, the 
present calculations of the charge accumulation are in qualitative agreement 
with those of Arafa (1949) and Huang (1948) respectively. 

For polyvalent solutes it is argued that all of the valence electrons must 
not enter the conduction band, for the heats of solution would then be positive. 

Since the change in energy of both the solute and solvent ions cannot be 
treated by this model, it is not to be expected that the results of this calcula- 
tion may be compared quantitatively with experiment. The model does, how- 
ever, reproduce qualitatively the general properties observed experimentally. 

The author wishes to express his appreciation to Mr. R. Talman_ for 
discussions. 

APPENDIX 
1. The Applicability of the Perturbation Method 

Rojansky (1946) points out that a general rule for the convergence of the 
perturbation method is that the change in the energy of a level should be 
small in comparison with the separation between levels. The separation 
between the levels AW in the Sommerfeld model of a metal is equal to the 
inverse of the density of states N(@) which is given by Mott and Jones 
(1936), for example, and hence 


2 3/2 
(A.1) AW = =e ) a. 


Q \2m 








622 CANADIAN JOURNAL OF PHYSICS. VOL. 35, 1958 


where { is the volume and € is measured from the bottom of the band. For 
copper the maximum value of € is 7.1 ev. and hence, in a gram-atom of 
copper, the minimum separation of the levels is, by equation (A.1), 14 10~*4 
ev.; about two per cent of the levels have a separation twice as great as 
this. At the bottom of the band the spacing is of the order of 1X107'* ev. In 
three dimensions, the first-order perturbation change in the energy level is 

w ye, 


Vy4= iV, = 1.5X10°-“V,, 


») wt eee 
(A.2) 70 71X10 


where w» and Q are the volume per atom and volume per gram-atom of copper 
respectively. It is, therefore, to be expected that for perturbations of a few 
electron volts, reasonable results will be obtained for levels near the top of 
the band. For levels at the bottom of the band, quite large perturbations 
may be treated by this method. Since the ratio of the first-order perturbation 
to the separation between the levels is independent of the volume of the 
metal, these conclusions hold irrespective of the size of the metal specimen. 


2. Evaluation of W,\” and the Generalization to Three Dimensions 
Equation (2.5) shows that 


(A.3) W,? = Wiet+ AE +, sin( en —##8) | 


If the levels are considered in pairs, so that k = 2m and k = 2m-+1, then 


Wem+ Wome 


= Wt Weel 2 4 of amndy mts 


9 L*L on L *Gnti L 
Ms pasty _— fant) s | 
(A.4) ee L 3! (Qm+1)x zc grees 
1 yO yo ro 6 
= 3[Wont Wensilt Vp? 5 - 


Thus, if the levels are taken in pairs, the mean energy of a level is altered 
by V,°6/L. 
In two dimensions 


(A.5) Vy° = V,°(x) . Va°(y) 


and the normalized wave functions are 


° 4 \).. kav. lay 
(A.6) Ver = (2) sin C sin - , 


where L, and Ly are not in the ratio of integers, which ensures no degeneracy 
in the system. The first-order perturbation treatment involves products of 
integrals of the form 


Lehi 


lry 


(A.7) Hi, = f sin? *™* V,(x)dx J sin’ F> V,(y)dy. 








tS 
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If the integrals are arranged in a two-dimensional array 


1 ww. | [6 


(A.8) ; : : 

r fs “ home pe 
(N/2)?1] (N/2)?(N/2)? 

and the integrals then summed first horizontally and then vertically, it is 

seen, by equations (A.4) and (2.5), that the change in energy of the N/2 

levels is given by NV’,°6,6,/L,L,. The average change in energy of a level 

is thus V,°6,6,/L,L,. Similarly, in three dimensions, it is found that the 

average change in energy of a level is V,°6,6,6,/L,L,L:. 

3. Discussion of the 1 pproximation (3.2) 

For |W,°| < |W;,°|, that is, for levels above W,°, 


, Wid Hid] We (=) | 
AY ee a ss i ses Sica ae 
itt ewe 2 we Litwetlhwe) + 
For |,°| > |W,°|, that is, for levels below JW,°, 

_ Heh, j Hetil Wo (x) | 
A.10 —— = : ; : eee : 
(A10) pene = Te ae Litwetlive) ++ Ms 
The Lennard-Jones approximation is equivalent to neglecting all the terms 
in both equations (A.9) and (A.10) except 


sy ists (H").2 (His)’ 


i r We We We 


This is also equivalent to neglecting W,° in the denominator of 


(A.12) >.’ HA, / (W.-W). 
t 
Inspection shows that 


| , xl , , | 
(A.13) | 2’ MixlTe | < b> Ay: | 


oe W, | me We-WwW? 
and both terms are negative. Thus 
(A.14) >! AL, ,/Wee 
i>k 


is a lower limit for the complete expression. It is also true that, near IV;,°, 


: |<, HaHis| . | ~~ Helis 

(A.15) | oe < | ESS, 

I<k W, l<k W, —W l | 

where the term on the right-hand side is positive and the term on the left, 


which is to be added to the term on the left-hand side of the relation (A.15), 
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is negative. The question now arises as to whether or not this latter term 
can make 


(A.16) HH, /We 





> | YS Hy, H,,/ (We We) |. 


The result obtained by Huang (1948), and referred to in Section 3, indicates 
that for the Sommerfeld model of a monovalent metal it does not, and thus 
the approximation is reasonable. 
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NOTES 


SHOCK-INDUCED PHYSICAL AND CHEMICAL SURFACE CHANGES 
ON OXIDE POWDERS 


R. W. NIcHOLLS AND W. H. PARKINSON 


The shock tube is a very convenient source of transient high temperatures 
(~ 104° K.) with which to excite spectra of solids which are difficult to excite 
in more conventional ways. A recent study of the luminosities produced in 
the shock excitation of about 20 inorganic oxides (Nicholls and Parkinson 
1957) (using helium-driven shock waves in argon up to Mach 7.5) showed 
that ‘clean’, entirely molecular spectra of, for example, BO, AIO, VO, and 
WO were produced at the lowest Mach numbers (4), mixed molecular and 
atomic spectra occurred at medium Mach numbers, and entirely atomic 
spectra occurred at the highest attainable Mach numbers (7.5). Subsequent 
experiments, using time-resolved emission and time-resolved flash absorption 
spectroscopy, have established that the reflected shock wave (from the quartz 
window at the end of the tube) is largely responsible for the excitation. 

It is the purpose of this communication to comment briefly upon unusual 
physical and chemical changes produced upon the surfaces of the shock- 
treated powder particles of a-hematite (Fe,O ;), litharge (PbO), and metatung- 
stic acid (H2W,Q);). 

The most obvious gross change is in color. FeO; turned from red to black, 
the intensity of the apparent blackening increasing with repeated shock 
treatment. PbO turned from yellow to gray, and H2W,O,; turned from yellow 
to blue. Microscopic examination showed that the seat of the color change 
was on the grain surface, either a uniform change across the whole surface 
(PbO, H2W,O,3) or a localized change at sporadic points upon the surface 
(FesOs). Microscopic evidence of severe ablation of the surfaces is seen in the 


case of FesQs. 
FeO; 


ligures 1(a) and 1(6) are microphotographs of typical FeO; grains before 
and after shock treatment at about 5700° K. (calculated temperature of the 
reflected shock wave) for about 2 milliseconds. The black spherical inclusions 
of about 10 4 diameter upon the ablated surface of the shocked grain are 
responsible for the color change. Further, repeated moderate shock treatment 
(.Wincia = 5) was found to be more effective in increasing the surface density 
of the black inclusions than repeated treatment by stronger shock waves 
(\Vinea = 7.5). The luminosity arising from the mild shock excitation was 
almost entirely of FeO, 4 and B band systems (Pearse and Gaydon 1950), 
while the luminosity arising from the stronger shock treatment was almost 
entirely of Fel lines. The shocked powder was also magnetically quite active. 
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A number of the black spheres were detached from parent grains and found 
to be the seat of magnetic activity. X-Ray powder photographs of the black 
spheres showed them to be Fe;O, (magnetite or ferrosoferric oxide FeO.Fe.O;). 
This was confirmed by a wet chemical analysis, which established the presence 


of the ferrous ion in the spheres. 





Fic. 1. Microphotographs at equal magnification of a-Fe:O3 grains before and after shock 
treatment. (a) Unshocked a-FesO3 grain. (6) Shocked (at Mach 5) a-FesO3 grain showing 
ablation and the black spherical inclusions of Fe;Oy. 


It is tempting to speculate that the Fe;0, spheres were formed by a conden- 
sation of FeO radicals (which are certainly copiously present in the luminous 
vapor produced by the shock excitation as shown by the spectrum) upon the 
Fe.O; surface with the formation of FeO.Fe.O;. Whether the FeO radicals 
were ablated from the grains by mild shock waves or whether individual atoms 
(Fe, O) were ablated and later combined to form FeO is a matter for further 
investigation. Similar observations were made upon ilmenite (FeTiO;); the 
spectrum consisted of FeO (.1 and B systems) and TiO (alpha and beta 
systems) bands (Pearse and Gaydon 1950). Fe;0, was again present on the 
grain surface after shock treatment. These observations may have some 


application to paleomagnetism. 


PbO 


Shock excitation of PbO powder at about WW inca = 5 (reflected temperature 
5700° IX.) gave rise to luminosity which consisted almost entirely of the 
A'S+—X'E+ and B'S+—X'Z+ band systems of PbO (Pearse and Gaydon 
1950). Even under these relatively mild excitation conditions a few PbI lines 
were present and the observed band systems were weak. A color change from 
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yellow to gray was produced. Microscopic observation showed that the yellow 
grains of PbO (100 4 diameter) had been replaced by small spheres of 10 u 
diameter of two types: shiny gray and fused yellow. The latter were probably 
fused PbO. X-Ray powder photographs of the former showed that they were 
either Pb or Pb+PbO (Sidgwick 1950). In this case the transient heating 
caused by the reflected wave evidently caused significant dissociation in 
addition to excitation. 


H2W, 013 


Repeated mild shock treatment (Wineia = 5) of this compound caused a 
color change from yellow to blue. The spectrum of luminosity consisted 
almost entirely of a complex system of bands attributed to WO (Pearse and 
Gaydon 1950; Gatterer and Krishnamurthy 1952). The grains appear to be 
uniformly blue after shock treatment, and it seems likely that one of the higher 
tungsten oxides, e.g. W4,0,,, has been formed. 


The work is proceeding, particular attention being paid to the quantitative 
aspects of ablation. Particle size and shape measurements which may have 
application to meteor physics are in progress. 
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THERMOELECTRICITY AT LOW TEMPERATURES 
V. THE SUITABILITY OF LEAD AS A STANDARD REFERENCE MATERIAL! 


J.-P. Jan,” W. B. PEARSON, AND I. M. TEMPLETON 


INTRODUCTION 


In discussing the determination of an absolute scale of thermoelectric power 
at low temperatures, Pearson and Templeton (1955, III in this series) con- 
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cluded that lead was probably the most suitable conductor to adopt as a 
standard reference metal. A redetermination of this scale for lead will be 
published shortly (Christian et a/. 1958, VI in this series). However, before 
this redetermination of the absolute thermoelectric power was undertaken, 
it was necessary to establish that small amounts of common impurities in lead 
did not give rise to a large anomalous thermoelectric power at low temperatures. 
(The thermoelectric power of copper and gold at low temperatures is, for 
instance, very strongly affected by traces of certain impurities; as an example, 
the presence of only ~0.002 wt.% of iron changes the absolute thermoelectric 
power of copper from ~0 to —9 uv. per degree at 15° K.) It was also considered 
advisable, in view of the work reported by Preston-Thomas (1952), to establish 
that the presence of small amounts of impurities in the lead did not appreciably 
modify the thermoelectric properties in zero magnetic field in the temperature 
region immediately above the superconducting transition temperature. The 
present paper gives a brief account of these investigations. 




















EXPERIMENTAL MEASUREMENTS AND DISCUSSION 










Measurements of the thermoelectric force, as a function of temperature, of 
dilute lead alloys containing bismuth, tin, indium, or cadmium in solid solution, 
shown in Figs. 1 and 2, were made by an improved form of the method already 
described by MacDonald and Pearson (1953, 1954). The thermoelectric force 
was measured relative to ‘‘spec-pure’’ lead (Johnson, Matthey & Co.) from 
which also the alloys were prepared. Alloys containing ~0.1 at.% Pd, or Au 
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Fic. 1. Thermoelectric force of dilute Pb-Sn and Pb-Bi alloys relative to lead. Abscissa: 
temperature of “hot’’ junction; ‘‘cold” junction at 4.2° K. 

Fic. 2. Thermoelectric force of dilute Pb-In and Pb-Cd alloys relative to lead. Abscissa: 
temperature of “hot’’ junction; “‘cold’’ junction at 4.2° K. 
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showed practically no net thermoelectric force relative to pure lead, but 
‘residual’ resistance measurements (Ise x./Ras5°xK.), Which were made as a 
check on the nominal compositions of all specimens, indicated that very little 
of these solutes actually entered into solid solution in the lead. 

In order to search for any evidence of ‘‘foreshadowing”’ of the superconduct- 
ing transition point, samples containing 0.04 at.% Bi and 0.04 at.% Cd were 
selected so that alloys were examined in which the electron concentration of 
lead had been both increased and decreased. Careful thermoelectric and resisti- 
vity measurements were made throughout the transition region (Figs. 3, 4, 
and 5), the absolute thermoelectric measurements being made against the 
superconducting alloy Nb;Sn by the method described in Christian ef al. 


(1958). 


TEMP °K 
t.2 7.3 


i 2 @ 6M 6B © Fe 


Fic. 3. Observed absolute thermoelectric force of a Pb — 0.04 at.% Bi alloy (A) and a 
Pb — 0.04 at.0¢ Cd alloy (O) in zero field in the region of the superconducting transition 
temperature. Independent measurements of the transition temperatures of these alloys by 
an inductive method show that they actually transform at the same temperature as pure lead, 
i.e. 7.17;+0.010° K. 

Fic. 4. Absolute thermoelectric power of a Pb — 0.04 at.% Bi alloy (A) and a Pb- 0.04 
at.©; Cd. alloy (O) measured against the superconductor NbsSn. -—— Absolute thermoelectric 
power of Johnson, Matthey and Co. Ltd. ‘“spec-pure”’ lead from which the alloys were made. 
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The data shown in Figs. 1 and 2 indicate that small amounts of com- 
mon impurities dissolved in lead do not give rise to any notable anomalous 
variation of thermoelectric power as a function of composition or of tempera- 


TEMR °K 
7.0 7.2 7.4 7.6 7.8 8.0 8.2 8.4 8.6 8.8 


RESISTANCE 


INDUCTIVE 
MEASUREMENTS 





_ Fic. 5. Resistive and inductive measurements of zone-purified lead and dilute alloys 
in the Vicinity of the superconductive transition. 


O zone-purified lead \ 
A Pb- 0.04 at.% Bi alloy 
O Pb-—0.04 at.{( Cd alloy 


resistive measurements, 
current | or 2 ma. 


© zone-puritied lead inductive measurements 


ture, and only in the alloy containing 0.1 at.% Cd does the thermoelectric 
power relative to pure lead become fairly large, being about 0.6uv./° C. 
between 10° and 20° kK. It is thus clear that in the amounts in which these 
metallic impurities are present in pure lead (readily obtainable in a purity 
exceeding 99.99%), they should have no noticeable effect on its observed 
absolute thermoelectric power. Furthermore the impurities which we have 
chosen to examine are representative of elements lying in Groups I to V and 
VIII. The other common impurities of Groups II to V, if they enter solid 
solution in lead, would be expected to behave similarly to the elements we 
have studied, while the Group I elements and transition metals probably do 
not dissolve in lead. The data given in Figs. 3, 4, and 5 do not reveal any 
serious ‘‘foreshadowing”’ effects above the superconducting transition tempera- 
ture, although a small effect may exist in the Pb—Bi alloy. We can therefore 
conclude (see also Pearson and Templeton 1955) that lead seems to be a 
most suitable reference material against which measurements of the absolute 
thermoelectric power of metals and alloys can be made. 

We are grateful to Dr. D. kK. C. MacDonald for his interest, and for dis- 


cussions of this work. 
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THE VALENCE STATE OF METALLIC LEAD! 


W. B. PEARSON 


A connection between cohesion and the effective sizes of the atoms in 
elemental metals has long been recognized by metallurgists; it provides, for 
instance, the basis for the conclusion that lead is divalent rather than quadri- 
valent in the metallic state (cf. Hume-Rothery and Raynor 1954). The failure 
of metallic lead to ionize the two electrons in ‘6s orbitals’’ is assumed from 
the abnormally large increase in the diameter of the metallic atoms found at 
Pb (and TI) on proceeding along the Third Long Period of the Periodic 
Table. Such a conclusion is also consistent with the chemical evidence that 
the divalent state of lead is more stable than the quadrivalent state, and 
that lattice-spacing measurements of lead alloys (Tyzack and Raynor 1954) 
do not give any indication that the anomalous size of the metallic lead atom 
can, as in the case of aluminum, be interpreted in terms of “Brillouin zone 
effects” (Axon and Hume-Rothery 1948). 

It appears to us difficult, however, to reconcile this view with the electrical 
properties of some dilute lead alloys, containing Cd, In, Sn, or Bi as solutes, 
which were examined at low temperatures. If the 6s? electrons (of the free 
atom) remain bound in the solid state and are not concerned with cohesion 
(i.e. do not take part in the formation of resonating metallic bonds) we would 
expect lead to behave as if only two electrons were ionized when it is con- 
sidered as an electrical conductor. From our experience of the behavior of 
heterovalent solutes in copper, silver, and gold (cf. Linde 1932) we would in 
these circumstances expect a homovalent impurity for lead to be one which 
also ionized two electrons. The data shown in Fig. 1, which are taken from 
thermoelectric and resistive measurements of dilute lead alloys containing 
various solute proportions (Jan et a/. 1958), indicate that a homovalent impurity 
in lead is an atom lying in Group IV, not in Group II, since the electron 
scattering increases with valence difference on proceeding from Group IV to 
higher or lower groups. It seems worth while to draw attention to this fact, 

‘Issued as N.R.C. No. 4719. 


Can. J. Phys. Vol. 36 (1958) 





632 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


particularly as the data of Fig. 1 were derived from alloys which contained 
only 0.008 to 0.1 at.% of each solute, and so it is unlikely that the electronic 
structure of the lead will have been changed by alloying. 
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THE RATIO OF THE NUCLEAR g-FACTORS FOR B'" AND B®! 


W. Brouwer? AND H. E. PETCH 


The nuclear magnetic resonance (n.m.r.) and atomic hyperfine structure 
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methods for the measurement of a nuclear g-factor depend upon the inter- 
action of the nuclear magnetic moment with an applied external magnetic 
field or with an internal atomic field, respectively. In the n.m.r. method the 
ratio of the resonance frequencies for two isotopes gives the ratio of their 
nuclear g-factors directly since the applied external field cancels in the ratio. 
However, for some pairs of isotopes the terms containing the effects of the 
internal fields differ slightly. In such cases, the ratio of the A's, the internal 
interaction constants, obtained by the atomic hyperfine structure method is 
not exactly equal to the ratio of the g’s. Deviations between the values deter- 
mined by the two methods, which lead to the so-called h.f.s. anomalies, have 
been interpreted (Bohr and Weisskopf 1950) as due, at least in part, to the 
effect of the finite size of the nucleus. 

It was recently pointed out to us that, to establish the h.f.s. anomaly for 
boron, a more precise value is needed for the ratio of the nuclear g-values 
of B" and B'. The value for the ratio of the .1’s calculated from a hyperfine 
structure measurement by Title (1957) is 


A(B")/A(B"") = 2.9867+0.0002. 
This value is to be compared with Bitter’s (1949) n.m.r. value of 
g(B")/g(B') = 2.9861 +0.0009. 


The ratio of the -1’s agrees with that of the g’s to within the latter’s rather 
large limits of error so that there appears to be no h.f.s. anomaly. The object 
of this work was to obtain, using n.m.r. techniques, a more precise value of 
the ratio of the g's for B" and B'® with which to compare the value for the 
ratio of the .1’s obtained from h.f.s. measurements. 

The ratio of the nuclear g-factors for a pair of isotopes is given directly 
by the ratio of the resonance frequencies of the two nuclear species in an 
externally applied magnetic field. The g-factor ratio can often be determined 
to a much greater accuracy than either of the g-factors alone because the 
frequency ratio for isotopic pairs is not subject to several corrections which 
normally have to be applied. To obtain the greatest possible accuracy, the 
experimental arrangement was designed to permit simultaneous measurement 
of the two resonance frequencies. Two radio-frequency coils were wound, 
with their axes mutually perpendicular, on a spherical container which held 
2 cc. of sample. Each coil formed the inductive branch of the tank circuit 
of a separate marginal oscillator. Resonance conditions were made repetitive 
by modulating the steady magnetic field with an amplitude greater than the 
signal width and a frequency of 60c./sec. The recurrent resonance signals 
from the two marginal oscillators were, after detection and amplification, 
simultaneously displayed on a Dumont Type 322-A double-beam oscilloscope. 
The horizontal traces of the oscilloscope were driven by the same sinusoidal 
signal that modulated the steady magnetic field. The frequency of the mar- 
ginal oscillator used to detect the B' resonance was maintained constant at 
20 Mc./sec.; it was continuously checked against the 20 Mc./sec. output from 
a General Radio unit time/frequency calibrator Type No. 1213-C which in 
turn was periodically calibrated by means of the 10 Mc./sec. signal from 
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WWYV. The steady magnetic field was then adjusted (H») ~ 14,640 gauss) 
until the B"™ signal from the 20 Mc./sec. oscillator was at the mid-point of 
one of the oscilloscope traces. The frequency of the second marginal oscillator 
was then adjusted until the B!° resonance signal appeared at the mid-point 
of the second oscilloscope trace. When exact matching of the signals had 
been attained, the frequency (6.69728 Mc./sec.) of the B'® marginal oscillator 
was measured with a Hewlett-Packard frequency counter Model 524A. This 
instrument was also checked against WWV at 10 Mc./sec. 

In a preliminary investigation, the boron isotopes were studied in a super- 
saturated aqueous solution of NaBO»s. Matching of the two signals was 
difficult because of the large disparity in intensity caused by the low natural 
abundance and small g-factor of B!® relative to B'". A specimen was therefore 
prepared by allowing powdered boron, enriched in B!°, to react with fused 
NaOH. A supersaturated aqueous solution of the reaction product gave good 
B!’ and B" signals of nearly equal intensities and total widths of approximately 
1.9 gauss and 0.9 gauss respectively. At the same position the field inhomo- 
geneity limited the proton line width to about 0.3 gauss. 

The mean ratio of the resonance frequencies, obtained by three different 
observers in this laboratory, is 


f(BY) /f(B") = ¢g(B")/g(B™) = 2.98629+0.00002. 


The calculated probable error corresponds to a matching of the two signals 


to within about one-tenth of the B" line width which seems quite reasonable 
in view of the difficulty of exactly matching two signals of different widths. No 
evidence of systematic errors was found. This more precise value of the ratio 
of the g-factors of B" and B!° is in good agreement with the value previously 
determined by n.m.r. methods by Bitter (1949) but it differs from the ratio 
of Title’s (1957) 1’s by an amount which is probably significant. Confirmation 
of this point must await a more precise measurement of the hyperfine structure. 


The writers are grateful to Dr. Hin Lew of the National Research Council 
of Canada who suggested this problem and supplied the sample of B'®. We 
are also obliged to Drs. Duckworth, Graham, and Kleerekoper, all of McMaster 
University, whose contributions of advice or equipment materially aided the 
progress of this work. We are happy to have this opportunity to acknowledge 
the generous contribution made by the Ontario Research Foundation for the 
purchase of the large electromagnet. 
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COSMIC-RAY INTENSITY VARIATIONS AT SEA LEVEL DURING 
MAGNETIC-STORM PERIODS 






G. RAMASWAMY* AND S. D. CHATTERJEE 









In this paper we will discuss some of the geophysical aspects of cosmic 
radiation in the light of ion-chamber measurements taken at Ottawa during 
February—May, 1951 (Chatterjee and Bloom 1955). This period was chosen for 
two reasons. Firstly, though the year 1951 was relatively quiet magnetically, 
the magnetic activity was higher during the first half than in the second half. 
Secondly, during the same period, intensity measurements were made at 
Berkeley, California, by Chasson using wide-angle triple-coincidence-counter 
telescopes, and an analysis of magnetic-storm effects on these intensity 
variations has been reported (Chasson 1954). Since both data were obtained 
at sea level, some comparison could be validly attempted. 

It is well known that the time variations in cosmic radiation are due in 
major part to the changes in the distribution of atmospheric mass. By applica- 
tion of what is known as ‘barometric correction’ to the sea-level intensity 
much of the effects of these changes can be removed and the corrected intensity 
may then be examined for possible magnetic-storm effects. 

Figure 1 shows the variations of cosmic-ray intensity during February 
May, 1951, when corrected for barometric pressure. The intensity variations 
are plotted as daily percentage deviations from the prestorm average. This 
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Fic. 1. Variations in cosmic-ray intensity corrected for barometric pressure. 









type of continuous plotting of variations is preferable to short-interval charts 
used by some earlier investigators like Hess (1940). For, one cannot say from 
short-interval charts whether the intensity changes are but part of the general 
pattern of variations or are different from it. A most noticeable feature of this 
curve is the general irregular pattern in the variations. Whether this feature can 
be attributed to the pressure—altitude effect (due to decay of u-mesons) 
cannot be answered with definiteness. However, as Chasson (1954) has ob- 
served, a simple barometric correction would give a set of corrected intensities 
reliable enough to be used in examining for magnetic-storm effects. 
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Cosmic-ray intensity decreases, immediately after the commencement of 
magnetic storms, are observable from February 21, March 6, April 18, and 
May 9. It might be said that these decreases are not within the general pattern 
of variations and are therefore physically real. These decreases can therefore 
be called simultaneous decreases and might be explained as being due to the 
magnetic storms. The decrease from May 1 is not sufficiently clear to be 
ascribed to the storm commencing May 1. The storm of March 13 does not 
seem to have an effect on the intensity, except probably to maintain the 
decrease for some days to come. The variation during the storm period com- 
mencing April 2 is interesting, for the decrease in intensity begins on March 
31 itself and continues for several days up to April 7. To call this phenomenon 
an example of prestorm effect may not be wrong, if the marked irregularities of 
the curve during this period are considered as being due to other factors or to 
random causes. 

Of the four cosmic-ray effects associated with magnetic storms, viz., (1) 
prestorm decrease, (2) simultaneous decrease, (3) delayed decrease, and (4) 
null effect, the most noticeable in our data appears to be the simultaneous 
decrease. However, examples of the other effects are also present in our data. 
Because of these exceptions no precise relationship of magnetic-storm varia- 
tions can be attempted. Chasson’s data give examples of all these effects and 
his résults are more inconclusive than ours. Consideration of our results with 
those of Chasson shows that there is no good agreement between these results. 
It is of interest to note that correlating the time variations and cosmic-ray 
intensity (corrected for barometric effect) with the corresponding indices of 
geomagnetic activity (obtained from Agincourt Observatory) gives an inverse 
relationship in the sense that their correlation coefficient is highly negative 
(its value being —0.81). This negative correlation suggests that there is a 
tendency for relatively lower cosmic-ray intensity on magnetically disturbed 
days than during quiet periods. This same conclusion was obtained by Forbush 
(1938) in an earlier analysis of Huancayo data by comparing the cosmic-ray 
intensity for quiet and disturbed days, averaged for several months. 

As yet, no theory constructed has been successful in explaining the behavior 
of cosmic-ray intensity during magnetic-storm periods. However, from a 
comparison of our results and the studies of several others it can be said that in 
a majority of cases there is a decrease in cosmic-ray intensity associated in some 
manner with the magnetic storms and their sudden commencement. Whether 
this decrease is simultaneous or whether this is characterized by any lag or 
lead in time can not be adequately described with a knowledge of current 


studies. 


The measurements were taken while one of us (S. D. C.) was holding a 
National Research Council Postdoctorate Fellowship at Ottawa. The writers 
wish to acknowledge with thanks the generous help given by Dr. D. C. Rose, 
in whose laboratory the measurements were taken. 
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LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


Nuclear Orientation of Mn*‘ in Antiferromagnetic Single Crystals 


It has been suggested by various people that nuclei might be aligned in an antiferromagnetic 
single crystal at low temperatures (e.g. Daunt 1951; Gorter 1957). In all the cases so far 
reported, interactions in general, and antiferromagnetism in particular, appear to reduce 
nuclear orientation rather than to enhance it. For example, the anisotropy of y-rays from 
Mn* in cerium magnesium nitrate is much reduced on account of interactions, and is reduced 
still further below 0.003° K., the Néel temperature of the salt (Grace et al. 1954). We wish to 
report that we have used the interactions in an antiferromagnetic single crystal to produce 
nuclear alignment, and claim that this represents a significant extension of the class of sub- 
stances in which nuclear alignment can be produced. 

Three single crystals of MnCls.4H.O, with a total weight of about 250 mg., containing 
about 3 microcuries of Mn*4 were mounted in parallel crystallographic orientation in a demag- 
netization cryostat. They were sandwiched between cylinders of compressed potassium chrome 


50 100 
Minutes from demagnetization 
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alum, and the space between was filled with a sludge of alum powder and Apiezon oil B. The 
whole sample was cooled by adiabatic demagnetization, and observations were made as the 
specimen warmed up over a period of 214 hours. No external magnetic field was applied during 
this time. Manganese chloride crystallizes in the monoclinic system, and its external morphology 
has been described by Groth (1906), whose nomenclature we use. According to Gijsman 
(1957) the direction of easy magnetization is the crystallographic c-axis, and the next preferred 
axis is the b-axis. We placed scintillation counters along the crystallographic b and ¢ axes and 
counted the y-rays emitted in those directions. Mn** decays by K capture followed by the 
emission of an E2 gamma ray of 835 kev. energy (Grace et al. 1954). We can define an anisotropy 
of y-ray emission e, by the equation e = (J, — J-)/ZJ» where, for example, J, is the intensity of 
y-ray emission along the b-axis. Measurements were also made of 7*, the magnetic tem- 
perature of the chrome alum, during the warm-up period. Figure 1 shows e and 1/7™ as func- 
tions of time from demagnetization. 

It is noticed that the manganese nuclei were in fact aligned along the c-axis. The magnetic 
moment of Mn* is so far unknown, so we are unable to estimate the internal field acting on the 
manganese nuclei. The time required for the anisotropy to develop, about half an hour, is not 
inconsistent with estimates of the spin-lattice relaxation time, made according to the theory of 
Moriya (1956). 

A similar effect is seen in MnBry.4H,0, and the results of an exactly similar experiment are 
shown in Fig. 2. The axis of preferred orientation in this salt is also the c-axis (Bolger 1955; 
Gijsman 1955). 


Hours from demagnetization 


Fic. 2. The salt was warmed up to the temperature of the helium bath after 10.6 hours, to take a normalization 
count. The normalization count is included on this figure to indicate the stability of the counters. The relaxation 
time in MnBr2.4H20 is longer than in MnCle.4H,O, probably because the Néel temperature of MnBrz.4H20 is 


higher than that of MnClh.4H,0. 


This work was carried out under contract from Atomic Energy of Canada Ltd., to whom 
grateful acknowledgment is made. One of us (M.A.R.L.) acknowledges the award of a National 


Research Council Studentship. 
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Long 7. Lifetimes from Positronium Annihilations in Liquid Argon and Liquid 
Nitrogen* 


The time spectra of positron decays in liquid argon and liquid nitrogen have been inves- 
tigated using a standard coincidence circuit (Bell, Graham, and Petch 1952). The source 
was a finely deposited layer of Na®Cl on 1.6 mg./cm.? aluminum foil and was placed in a monax 
cryostat of } in. external diameter. Two-inch-thick NalI(Tl) phosphors in conjunction with 
RCA 6342 photomultipliers were used as detectors. The resolving time, 270, of the circuit 
was 3.4X107-° sec. at 60% coincidence efficiency, and 5.2107 sec. at 90% coincidence effi- 
ciency, with mean drop-off rates of 8X10~" sec. In this way over 1000 true coincidences per 
minute were detectable without using a source so strong as to overwork the photomultipliers 
at the very high operating voltages necessary for such high resolution (~2400 volts). 

Measurements revealed 72 components, such as are normally associated with the forma- 
tion of triplet positronium, in both liquids; see Fig. 1 and Table I. The values of J: in Table I 
were calculated making an allowance for the absorption of positrons in the source materials: 


POSITRON ANNIHILATIONS 
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Fic. 1, Time spectra for positron decays in argon and nitrogen liquids of 99.99% purity. 


*This work has been assisted in part by the Defence Research Board of Canada under Grant No. 9510-11. 
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this allowance was measured experimentally. A calculated correction for backse< attering into 
the source materials is also included: the backscattering causes less than 10% of all positrons 
to annihilate in the source. 


TABLE I 
TRIPLET POSITRONIUM ANNIHILATION DATA 











i Lifetime re, Temperature, 

Substance 10~° sec. units Intensity J2, % *s 
Liquid argon 99.99% purity 6.9+0.7 7 7+1 87.55+0.15 
Liquid nitrogen 99.99% purity ll. 340. 5 UX 542 77.5 +0.1 








Errors in 72 are 2Xstandard deviation+an allowance for calibration errors. 
Errors in J: are merely estimates. 


No observable temperature effect was found in liquid argon, but this substance is only 
liquid over a range of 4 centigrade degrees. 

These 72 components are interesting on account of their long mean-life which is far in excess 
of values reported for any other liquids except liquid helium (see Table IT). De Benedetti and 
Siegel (1954) state that the number of three-quantum annihilations in liquid nitrogen is ‘‘of 
the same order as that expected from an equal number of positrons stopping in a solid metal”’ 
This observation is in agreement with previous observations of R. L. Graham and myself 
(1955, 1956, unpublished) who observed no 72 component in liquid nitrogen, indicating that 
either positronium was not formed or the triplet lifetime was highly quenched. These earlier 
results can now be understood if one assumes a small percentage of oxvgen (a very common 
impurity in commercial nitrogen) to have been present in both de Benedetti’s and our samples 
of liquid nitrogen. The effect of the oxygen is to reduce the 72 lifetime by converting triplet 
positronium to singlet; this arises as a result of the magnetic moment of the oxygen molecule. 
Recent experiments at this laboratory (1957) with commercial grade liquid nitrogen again 
revealed no 72 component. 


TABLE II 
TRIPLET POSITRONIUM ANNIHILATION IN PURE LIQUIDS 





Temp., 72, 

L iquid =k. 10~° sec. units Ie, % Reference 
Isopropy! alcohol 293 2.3+40.4 ~30 Bell and Graham (1953) 
Water 293 1.7+0.2 ~30 Bell and Graham (1953) 
Water 293 1.8+0.1* 2145* Green and Bell (1957) 
Benzene 293 2.67 +0.08 3143 Berko and Zuchelli (1956) 
p-Dioxane 293 2.67 3647 Green and Bell (1957) 
Naphthalene 353-373 2.68 29 Landes, Berko, and Zuchelli (1956) 
Ethylenediamine 1.8+0.2 Hogg, Sutherland, Paul, and Hodgins (1956) 
Ammonia 1.7+0.2 Hogg, Sutherland, Paul, and Hodgins (1956) 
Helium I 120 +20 1642 Wackerle and Stump (1957) 
Helium I 91 +3t 15 5+2 2.é Paul and Graham (1957) 
Helium II 84.5 +4t mad 5 +2. Paul and Graham (unpublished) 





*Authors spe cifie d that thes se are standard errors. 
tErrors are 2 Xstandard deviation (statistical) +allowance for calibration errors. 
tThe value published by Paul and Graham was 13% but was in error owing to an arithmetical mistake. 


Deutsch (1952), working with argon gas, observed that oxygen would induce ortho-para 
conversion in positronium with a molecular cross-section ¢ = 4X107" cm.? This value will 
have been determined by assuming that positronium was in thermal equilibrium with the 
gas, and substituting the appropriate value of v, the positronium velocity, in 


ov = 3.1XK10-” cm.3/sec 


The product ov will be referred to as the “volume rate’’. A study of the effect of oxygen in 
liquid argon has been made using a new five-channel adaptation of the circuit of Bell, Graham, 
and Petch (1952). With this circuit the counting rates at five points on the 72 “tail” of a liquid 
argon resolution curve are simultaneously observed. Figs. 2 and 3 show some results. The 
quenching effect of oxygen is evidence that the 72 component is due to the formation of triplet 
positronium. The value for the ‘‘volume rate” due to oxygen in liquid argon is computed from 
Fig. 3 to be ov = 9.7107 cm.3/sec. per molecule, more than three times Deutsch’s value. 
While there is no a priori reason why these two observed volume rates should agree, they 
would have the same value if ¢ varied inversely with positronium velocity. If therefore a simple 

gas kinetic picture for positronium collisions in liquid argon is justifiable, it can be concluded 
that « does not follow a 1/v law. Volume rates have been compared here instead of cross sec- 
tions in order to avoid making any assumption about thermalization of positronium in liquid 


argon. 
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The wide difference in the two volume rates would be clarified by an investigation of the 
variation of o with v. Measurements on the quenching of triplet positronium by oxygen in 
ov 1,0 1” - A ac: i 
argon gas over a range of temperatures, say 100° K. to 270° K., would yield o as a function 
of v, at least over a limited range of positronium velocities. 
Measurements on liquid nitrogen are being continued at this laboratory. 
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